Strong QCD and Dyson- Schwinger Equations 



(N 

o 



Craig D. Roberts 

Physics Division, Argonne National Laboratory, 

Argonne, Illinois 60439, USA; 

and Department of Physics, Illinois Institute of Technology, 

Chicago, Illinois 60616-3793, USA. 









en 

rn 
O 
(N 



% 



Abstract. The real- world properties of quantum chromodynamics (QCD) - the strongly- 
interacting piece of the Standard Model - are dominated by two emergent phenomena: 
confinement; namely, the theory's elementary degrees-of-freedom - quarks and gluons - 
have never been detected in isolation; and dynamical chiral symmetry breaking (DCSB), 
which is a remarkably effective mass generating mechanism, responsible for the mass of 
more than 98% of visible matter in the Universe. These phenomena are not apparent 
in the formulae that define QCD, yet they play a principal role in determining Nature's 
observable characteristics. Much remains to be learnt before confinement can properly 
be understood. On the other hand, the last decade has seen important progress in the use 
of relativistic quantum field theory, so that we can now explain the origin of DCSB and 
are beginning to demonstrate its far-reaching consequences. Dyson-Schwinger equations 
have played a critical role in these advances. These lecture notes provide an introduction 
to Dyson-Schwinger equations (DSEs), QCD and hadron physics, and illustrate the use 
of DSEs to predict phenomena that are truly observable. 
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2 Craig D. Roberts 

1 Introduction 

The context for this contribution is provided by recent progress in mathematics, which 
followed a realisation that the process of renormalisation is naturally expressed via a Hopf 
algebra structure [1] . Indeed, the Hopf algebra approach allows for a comprehensive de- 
scription of the algebraic and combinatorial structures that underpin renormalisation, 
which is critical in four-dimensional quantum field theory. This enables a mathemati- 
cally sound approach to the problem of computing the /3-function, which is based on two 
elements: the existence of quantum equations of motion, namely, the Dyson-Schwinger 
equations (DSEs); and the consequences of the renormalisation group for local field the- 
ories. One needs the DSE umbrella in order to guarantee sufficient recursive structure in 
the theory such that a nonperturbative approach becomes feeisible. The Hopf algebraic 
foundations of these phenomena make the approach possible. Finally, the Hopf algebra 
description permits one to comprehend and explore basic ideas of renormalisation, driv- 
ing new applications of those ideas in the context of pure and applied mathematics. The 
progressive mathematical reformulation of the established physical procedure of renor- 
malisation has motivated people to imagine that this framework has the power to provide 
deeper insights into fundamental problems in quantum field theory. However, practical 
successes are rare [2]. 

This contribution is primarily an overview of the use of DSEs in hadron physics; i.e., 
in nonperturbative quantum chromodynamics (QCD). As such, the aims and problems 
are quite distinct from those discussed in connection with perturbation theory. Hadron 
physics practitioners are primarily concerned with finding a truncation of the infinite 
tower of DSEs that provides a realistic description of bound-states, rather than with 
bringing order to the recursive tower itself. A goal of these lecture notes is therefore to 
highlight the practical problems faced in the real-world use of DSEs and the physical 
connections between formalism and phenomena. 

A springboard to the widespread application of DSEs in hadron physics was provided 
by Ref. [3]. The approach has since been successful. Hence numerous other reviews have 
subsequently been written; e.g., Refs. [4-15]. Whilst those articles are a resource, it 
should not be necessary to read them before delving into this one. 

QCD is the strong-interaction part of the so-called Standard Model of Particle Physics. 
It concerns length-scales on the order of lfm= 10~^^m, or smaller; i.e., scales less than 
one femtometre (one "fermi"). This corresponds to an energy of Aqcd — 0.2 GeV or 
greater.^ Perturbation theory is a valid tool within QCD at energies > 10 Aqcd. It is 
widely hoped that phenomena outside the realm of the Standard Model will be encoun- 
tered at energies above 1 TeV = 1000 GeV. 

In order to understand the origin of the Standard Model, one should recall that in the 
early twentieth century the only matter particles known to exist were the proton, neutron 
and electron. However, with the advent of cosmic ray science and particle accelerators, 
numerous additional particles were discovered: muon (1937), pion (1947), kaon (1947), 
Roper resonance (1963), etc., so that by the mid-1960s it was apparent that not all the 
particles could be fundamental. A new paradigm was necessary. The constituent-quark 
theory of Gell-Mann [16] and Zweig [17] was a critical step forward. This was recog- 



'^"eV" is the abbreviation for electron- volt, which is used as a standard unit of measure in 
physics because of its usefulness in particle accelerator science: a particle with charge q has 
an energy E = qV after passing through a potential V; if q is quoted in units of the positron's 
charge and the potential in volts, one obtains an energy in eV. 1 GeV = lO^eV. Hadron physicists 
typically work with a system of "natural units" in which h = 1 = c, so that length- and mass- 
energy-scales are related as follows: Ifmsi I/Aqcjj. 



Strong QCD and Dyson-Schwinger Equations 3 

nised with the Nobel Prize to Gell-Mann in 1969: "for his contributions and discoveries 
concerning the classification of elementary particles and their interactions". Over the 
more than forty intervening years, the theory now called the Standard Model of Particle 
Physics has passed almost all tests. 

There are three pieces to the Standard Model. The first is the theory of electro- 
magnetism; i.e., quantum electrodynamics (QED), which was developed in the period 
1946-1950 by Feynman, Schwinger and Tomanaga, who shared the 1965 Nobel Prize in 
physics: "for their fundamental work in quantum electrodynamics, with deep-ploughing 
consequences for the physics of elementary particles" [18]. The second, developed in the 
period 1963-1973, consumed QED through the unification of a description of electromag- 
netism and weak interactions, such as radioactive decays, parity-violating decays and 
electron-neutrino scattering. This body of work, which produced the electroweak theory, 
was recognised with award of the 1979 Nobel Prize in physics to Glashow, Salam and 
Weinberg: "for their contributions to the theory of the unified weak and electromagnetic 
interaction between elementary particles, including, inter alia, the prediction of the weak 
neutral current" [19]. The final piece is QCD, which is supposed to describe: the existence 
and composition of the vast bulk of visible matter in the Universe - the proton, neutron, 
pion, etc.; the forces that form them; and the structure of nuclei. The development of 
QCD is far from complete but the theory is understood within the domain for which 
perturbation theory is a valid tool, following work by Politzer [20,21], Gross and Wilczek 
[22] in 1973-1974 that brought the 2004 Nobel Prize in physics: "for the discovery of 
asymptotic freedom in the theory of the strong interaction". As we shall see, however, 
hadron physics requires an understanding of nonperturbative QCD. 

Constituted from these elements, the Standard Model is a local quantum gauge field 
theory, which can be expressed in a very compact form. The Lagrangian possesses the 
following gauge symmetry: 

SUc{3) X SUl{2) X Uy{1) ^ ^. 

QCD electroweak 

and contains 19 parameters, which must be determined through comparison with ex- 
periment. Physics is, after all, an experimental science. The electroweak portion of the 
Standard Model involves 17 of the parameters. Most of these are tied to the Higgs boson, 
the model's only fundamental scalar, which has never been seen. The electroweak sector 
is essentially perturbative, so the parameters are readily determined. Only two of the 
Standard Model's parameters are intrinsic to QCD. However, as QCD is Nature's only 
example of a truly and essentially nonperturbative fundamental theory, the impact of 
these parameters is not yet fully known. 

The known particle content of the Standard Model is depicted in Fig. 1.1. The Higgs 
boson is obvious by its omission. Discovery of the Higgs is a primary mission for the 
Large Hadron Collider (LHC). Built at a cost of $7-billion, this machine accelerates 
two antiparallel beams of particles to almost the speed of light in a 27km tunnel that 
is 175m underground, before colliding them at four interaction points. During a ten 
hour experiment, each beam will travel 10-billion km; i.e., almost 100-times the earth- 
sun distance. The energy of each collision will reach 14TeV. It is the hope of many that 
this energy is high enough to expose the new particles or structures that are predicted 
by existing extensions of the Standard Model. There is also a fear that it is not. 

It is noteworthy that, at the time of writing, there is a growing acceptance that the 
Higgs-boson might not exist: the ATLAS and CMS detectors have excluded the existence 
of a Higgs over most of the mass region 145 — 466 GeV with 95% certainty. Optimists 
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Figure 1.1. Known particle content of the Standard Model. The Higgs boson is conspic- 
uously absent. It is possible that within nonperturbative QCD the distinction between 
matter- and force-fields disappears because bound-states formed solely from the force 
carrying gluons (gauge-bosons) are conceivable. 



view this as an opportunity, arguing that if the Higgs does not exist, its absence will 
point the way to new physics.^ 

There is certainly physics beyond the Standard Model; i.e., observational facts that are 
inexplicable within the model. For example, neutrinos have mass, which is not the case 
in the Standard Model. This was first detected in the Homestake mine experiment of the 
late 1960s [23,24], which observed a deficit in the flux of solar neutrinos Ue; i.e., neutrinos 
associated with the electron. Controversial from the first moment, this observation was 
verified conclusively at the Sudbury neutrino observatory [25] . It is explained by neutrino 
oscillations [26-28]; viz., so long as neutrinos have mass, a neutrino created in association 
with a specific lepton flavour (electron, muon or tau) can transform into that associated 
with another flavour: this is mixing between neutrino flavour- and mass-eigenstates, a 
phenomenon also seen in the weak interactions of quarks. 



^This was the case until Tues., 13 Dec. 2011. On that day, amidst great fanfare, CERN an- 
nounced that they believe they have seen a hint of the Higgs boson. The results are inconclusive. 
In that announcement, the following observations was made: "The Higgs field is often said to 
give mass to everything. That is wrong. The Higgs field only gives mass to some very simple 
particles. The field accounts for only one or two percent of the mass of more complex things like 
atoms, molecules and everyday objects, from your mobile phone to your pet llama. The vast 
majority of mass comes from the energy needed to hold quarks together inside atoms." This 
recognises the fundamental role of confinement and dynamical chiral symmetry breaking, key 
emergent phenomena within the strongly interacting part of the Standard Model, upon which 
much of this discourse will focus. 
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Neutrino masses are one concrete example of physics beyond the Standard Model. Our 
existence is another because we represent a macroscopic excess of matter over antimatter. 
Sakharov first observed [29] that in order to produce an excess of baryons, there must 
be processes that change baryon number. No such processes have yet been observed. 
Moreover, those processes must occur out of thermal equilibrium, otherwise they would 
merely balance matter and antimatter. (Actually, the electroweak component of the 
Standard Model is capable of satisfying Sakharov's conditions in principle but it is too 
weak to explain the observed matter-antimatter asymmetry.) The very likely existence 
of dark matter and dark energy are two more indications that the Standard Model is 
incomplete [30]. With this admission of incompleteness, one has opened Pandora's box: 
the known problems are few but the number of people working on physics beyond the 
Standard Model is practically innumerable. 

One means by which to gauge the foci of contemporary research in physics is to 
consider the top-problems lists that various communities have compiled: one can readily 
find a top-10 or top-24, etc. (Such lists have a long history, for who can forget Hilbert's 
1900 list [31] of the most important problems in mathematics for the twentieth century?) 
Here I choose to identify just five. 

What is dark matter? There appears to be a halo of invisible material engulfing 
galaxies, which is commonly referred to as dark matter. Existence of dark (=invis- 
ible) matter is inferred from the observation of its gravitational pull, which causes 
the stars in the outer regions of a galaxy to orbit faster than they would if only 
visible matter were present. Another indication is that we see galaxies in our own 
local cluster moving toward each other. 

What is dark energy? The discovery of dark energy goes back to 1998. A group of 
scientists had recorded several dozen supernovae, including some so distant that 
their light had started to travel toward Earth when the universe was only a fraction 
of its present age. Contrary to expectations, the scientists found that the expansion 
of the universe is not slowing, but accelerating [32,33]. This year, the leaders of 
those teams shared the Nobel Prize in Physics. 

What is the lifetime of the proton and how do \ve understand it? It was once 
supposed that protons, unlike, neutrons, are stable, never decaying into smaller 
pieces. Then, in the 1970's, theorists realised that their candidates for a grand 
unified theory, merging all the forces except gravity, implied that protons must be 
unstable. Wait long enough and, very occasionally, one should break apart. Is this 
truly a necessary consequence of Grand Unification? 

What underlies the enormous disparity between the gravitational scale and 
the typical mass scale of the elementary particles? Expressed differently, 
why is gravity so much weaker than the other forces? A small magnet can raise a 
paper clip even though the gravity of the whole earth is resisting. 

Can we quantitatively understand quark and gluon confinement in quantum 
chromodynamics and the existence of a mass gap? QCD is the theory de- 
scribing the strong nuclear force. Carried by gluons, it binds quarks into particles 
like protons and neutrons. Apparently, the quarks and gluons are permanently 
confined: one cannot remove a quark or a gluon from a proton because the re- 
straining force grows enormously with increasing distance and prevents them from 
escaping. 

This last brings me to the point of these lectures. 
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Figure 1.2. Left diagram - The leading-order contribution to light-by-light scattering is 
mediated by fermion-antifermion production and annihilation, which occurs at 0(a ); 
i.e., with probability ~ 10~®. Right diagrams - In QCD, gauge-boson self-interactions 
occur at tree-level, with these three-gluon [0((7)] and four-gluon vertices [0{g^)]. 



1.1 Quantum chromodynamics 

It is straightforward to express the local Lagrangian density that defines QCD: 

rqcD = g4i7^[^M].. -A«»j]*-^G^.Gr (1.2) 



= gi[n''9,i - m]qi - gC^qa^T^jqj - -01^0-^" (1.3) 



where G is the gluon gauge-field, q is the quark matter-field, g is the coupling constant, 
and the gluon field-strength tensor is 

G^. = d,Gl - d.Gl + gf^-^GlGl . (1.4) 

The underlined term in Eq. (1.4) generates gauge-field self-interactions, which have 
extraordinary consequences. This is readily elucidated through a comparison with QED, 
an Abelian gauge theory, which is perturbatively simple but, alone, undefined nonper- 
turbatively."^ A characteristic feature of QED is that gauge-boson self interactions do 
not occur at tree level. The leading contribution occurs at order a*. Fig. 1.2, and since 
OL ~ 1/137, such interactions take place with extremely small probability. In contrast, 
the non-Abelian character of QCD produces tree-level interactions between gauge bosons 
- the three- and four-giuon vertices generated by G^^^Ga" in Eq. (1.2); and one might 
guess this will have a big impact. Doing better than guessing led to the 2004 Nobel Prize 
in Physics [20-22]. 

Quantum gauge-field theories are typified by the feature that nothing is constant. 
The distribution of charge and mass, the number of particles, etc.; indeed, all the things 
that quantum mechanics holds fixed, depend upon the wavelength of the tool being used 
to measure them. Couplings and masses are renormalised via processes involving virtual- 
particles. Such effects make these quantities depend on the energy scale at which one 
observes them, and produce the running couplings illustrated in Fig. 1.3. At one-loop 
order, the QED coupling is 

oqed(Q) = 1^ Q-, (1.5) 

1 - -- In 

6ii me 

where a and tti^ are renormalised at the physical electron mass. The negative sign in 
the denominator signals that fermions screen electric charge. In QCD, the analogous 



^ Four-fermion operators become relevant in strong-coupling QED and must be included in 
order to obtain a well-defined continuum limit [34,35] . 
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Figure 1.3. Left panel - Running coupling in QED; and right panel - running coupling in 
QCD. N.B. On the same momentum domain, the QCD coupling changes 500 000-times 
more than the QED coupling and runs in the opposite direction. 



quantity is 



aQCD((5) 



(33-2iV/)ln- 



Q 



(1.6) 



Aqcd 

where Nf is the number of quark flavours that are active in the process under consid- 
eration and Aqcd is a dynamically generated mass-scale, whose value is experimentally 
determined; viz., Aqcd ~ 200 MeV. Plainly, quarks screen the colour charge. However, 
owing to their self-interactions, gluons antiscreen this charge. The enormous impact of 
that effect is evident in Fig. 1.3. 

The running coupling depicted in the right panel of Fig. 1.3 can be translated into 
an interaction between quarks and gluons that depends strongly on separation; i.e., an 
interaction which grows rapidly with separation between the sources. Taking the experi- 
mental value of Aqcd, the coupling is found to be huge at a separation r ~ 0.2 fm ~ jVp, 
where Vp is the proton's charge radius. This is a peculiar circumstance; viz., an interac- 
tion that becomes stronger as the participants try to separate. It leads one to explore 
some fascinating possibilities; If the coupling grows so strongly with separation, then 
perhaps it is unbounded; and perhaps it would require an infinite amount of energy in 
order to extract a quark or gluon from the interior of a hadron. Such thinking has led to 
the 

Confinement Hypothesis: 

Colour-charged particles cannot be isolated and therefore cannot be directly 

observed. They clump together in colour-neutral bound-states. 

This is hitherto an empirical fact. In the hope that it will become rather more, in 2000 
the Clay Mathematics offered a prize of $l-million for a proof that 5'L*'c(3) gauge theory 
is mathematically well-defined [36] , one necessary consequence of which will be an answer 
to the question of whether or not the confinement conjecture is correct. 

I would like to note here that pure-glue QCD is far simpler than real- world QCD, and 
that this is crucial in connection with the nature of confinement. Since bosons possess 
a classical analogue, it is relatively straightforward to define the pure-glue theory via 
a discrete probability measure that is amenable to direct numerical simulation using 
Monte-Carlo methods. In the absence of light-quarks, there is no perniciously nonlocal 
fermion determinant to interfere with the simulation. The quenched theory has long been 
subject to analysis and the "confinement of quarks" addressed in that context [37], with 
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the result that a Unearly-rising potential is known to exist between static sources. It 
is unfortunate, therefore, that this potential has no known connection with the physics 
of real- world QCD, wherein light-quarks are omnipresent and between which, owing to 
their low mass relative to the natural scale of QCD, it is impossible to define a potential. 
QCD is supposedly the theory underlying the strong interaction. It is asymptotically 
free, so that a chiral limit is defined and perturbation theory is a valid tool for Q^ 2> 
Aqqq. On the other hand, it is associated with an interaction that is strong over 98% 
of the interior of a typical hadron, such as the proton, and therefore impossible to treat 
using perturbation theory. QCD is Nature's only example of a truly nonperturbative, 
fundamental theory and, a priori, no one has any idea as to what such a theory is capable 
of producing. The study of nonperturbative QCD is the purview of Hadron Physics. 



1.2 Hadron physics 

The term "hadron," first used in 1962 [38], refers to a class of subatomic particles that are 
composed of quarks and/or gluons and take part in the strong interaction. Well known 
examples are the proton, neutron, and pion. The word is a derivative of the Greek word 
hadros {a5p6i;) , which variously means "large" or "massive," and is intended to contrast 
with leptos (AeTTTOi,), which means "small" or "light." The world of hadrons is subdivided 
into baryons, which are hadrons with half-integer spin; and mesons, hadrons with integer 
spin. 

Hadron physics is unique at the cutting edge of modern science because Nature has 
provided us with just one instance of a fundamental strongly-interacting theory; i.e., 
QCD. The community of science has never before confronted such a challenge as solv- 
ing this theory. In 2007, the USA's Nuclear Science Advisory Council delivered a long 
range plan to the Department of Energy's Office of Nuclear Physics, which included the 
statement: "A central goal of (the DOE Office of) Nuclear Physics is to understand 
the structure and properties of protons and neutrons, and ultimately atomic nuclei, in 
terms of the quarks and gluons of QCD." Internationally, this is an approximately $1- 
billion/year effort in experiment and theory, with roughly $375-million/year spent in the 
USA. On the order of 90% of these funds are spent on experiment. To provide a context, 
$l-billion/year is the operating budget of CERN. 

The hadron physics community now has a range of major facilities that are accumu- 
lating data, of unprecedented accuracy and precision, which pose important challenges 
for theory. The opportunities for researchers in hadron physics promise to expand with 
the use of extant accelerators, and upgraded and new machines and detectors that will 
appear on a five-to-ten-year time-scale, in China, Germany, Japan, Switzerland and the 
USA. A short list of facilities may readily be compiled: Beijing's electron-positron col- 
lider; in Germany - COSY (Jiilich Cooler Synchrotron), ELS A (Bonn Electron Stretcher 
and Accelerator), MAMI (Mainz Microtron), and FAIR (Facility for Antiproton and Ion 
Research) under construction near Darmstadt, with new generation experiments in 2015 
(or later, as now seems more probable); in Japan - J-PARC (Japan Proton Accelerator 
Research Complex) under construction in Tokai-Mura, 150km NE of Tokyo, expected to 
begin operation toward the end of 2012, and KEK, Tsukuba; in Switzerland, the ALICE 
and COMPASS detectors at CERN; and in the USA, both RHIC (Relativistic Heavy 
Ion Collider) at Brookhaven National Laboratory, which focuses on the strong phase 
transition, physics just 10/iS after the Big Bang, and the Thomas Jefferson National Ac- 
celerator Facility (JLab), exploring the nature of cold hadronic matter, currently being 
upgraded at a cost of $310-million, with new generation experiments expected in 2016. 
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The interpretation and prediction of phenomena at JLab provide the themes for much 
of the material presented herein. 

A theoretical understanding of the phenomena of hadron physics requires use of the 
full machinery of relativistic quantum field theory, which is the only way to reconcile 
quantum mechanics with special relativity. Relativistic quantum field theory is based on 
the relativistic quantum mechanics of Dirac and its unification with quantum mechanics 
took some time. Questions still remain as to a practical implementation of a Hamiltonian 
formulation of the relativistic quantum mechanics of interacting systems. 

The Poincare group has ten generators: six associated with the Lorentz transforma- 
tions (rotations and boosts); and four associated with translations. Quantum mechanics 
describes the time evolution of a system with interactions and that evolution is gener- 
ated by the Hamiltonian, or some generalisation thereof. However, the Hamiltonian is 
one of the generators of the Poincare group, and it is apparent from the Poincare algebra 
that boosts do not commute with the Hamiltonian. This is partly because the existence 
of antiparticles is concomitant with relativistic quantum mechanics; i.e., the equations 
of relativistic quantum mechanics admit negative energy solutions. However, once one 
allows for negative energy, then particle number conservation is lost: 

-Esystom = -Esystom + (-Epj + Ep-^) + . . . ad infinitum, (1.7) 

where Ej. = —E^- This poses a fundamental problem for relativistic quantum mechanics: 
few particle systems can be studied, but the study of (infinitely) many bodies is difficult 
and no general theory currently exists. Following from these observations, the state 
vector calculated in one momentum frame will not be kinematically related to the state 
in another frame, a fact that makes a new calculation necessary in every momentum 
frame. The discussion of scattering, which takes a state of momentum p to a different 
state with momentum p' is therefore problematic [39,40]. 

Relativistic quantum field theory provides a way forward. In this framework the fun- 
damental entities are fields, which can simultaneously represent infinitely many particles. 
The neutral scalar field, 4'{x), provides an example. One may write 



(27r)3 2^ 



[a(fc)e-*'=-"-ha^(/fc)e''=-"] , (1.8) 



where: uu ~ y \k\'^ + m? is the relativistic dispersion relation for a massive particle; 

the four- vector (k'^) — {ctJk,k); a{k) is an annihilation (creation) operator for a particle 
(antiparticle) with four-momentum k (— fc); and a' (fc) is a creation (annihilation) operator 
for a particle (antiparticle) with four-momentum k {—k). With this plane- wave expansion 
of the field one may proceed to develop a framework in which the nonconservation of 
particle number is not a problem. That is crucial because key observable phenomena in 
hadron physics are essentially connected with the existence of virtual particles. 

Relativistic quantum field theory has its own problems, however. For example, the 
question of whether a given quantum field theory is rigorously well defined is an un- 
solved mathematical problem. All relativistic quantum field theories admit analysis via 
perturbation theory, and perturbative renormalisation is a well-defined procedure that 
has long been used in quantum electrodynamics (QED) and quantum chromodynamics 
(QCD), as emphasised by the 2004 Nobel Prize in Physics [20-22]. However, a rigorous 
definition of a theory means proving that the theory makes sense nonperturbatively. This 
is equivalent to proving that all the theory's renormalisation constants are nonperturba- 
tively well-behaved. I have written that understanding the properties of hadrons requires 
solving QCD. However, QCD is not known to be a rigorously well-defined theory and 
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hence it might not have a solution. Thus, more properly, the task of hadron physics is 
to determine whether QCD is truly the theory of the strong interaction. 

Physics is an experimental science and experiment advances quickly. Hence, de- 
velopment of an understanding of observable phenomena cannot wait on mathematical 
rigour. Assumptions must be made and their consequences explored. Hadron physics 
practitioners therefore assume that QCD is (somehow) well-defined and follow where 
that idea may lead. That means exploring and mapping the hadron physics landscape 
with well-understood probes, such as the electron at JLab; and employing established 
mathematical tools, refining and inventing others, in order to use the Lagrangian density 
of QCD to predict what should be observable real-world phenomena. 

A primary aim of the world's current hadron physics programmes in experiment and 
theory is to determine whether there are any contradictions with what can actually 
be proved in QCD. Today, there are none which are uncontroversial. Indeed, in my 
view there are none and those who claim otherwise are misguided; see, e.g., Refs. [41- 
44]. Another is to identify what QCD will probably predict. In this field the interplay 
between experiment and theory is the engine of discovery and progress, and the discovery 
potential of both is high. Much has been learnt in the last five years - these notes will 
provide a perspective on the meaning of some of those discoveries - and one can safely 
expect that many surprises remain in hadron physics. 



1.3 Euclidean metric 

In the first formulation of (quenched-)QCD as a numerical problem, a Euclidean metric 
was employed [37]. The Euclidean-QCD action defines a probability measure, for which 
many numerical simulation algorithms are available. Working in Euclidean space is, how- 
ever, more than simply pragmatic: Euclidean lattice field theory is currently a primary 
candidate for the rigorous definition of an interacting quantum field theory. This relies 
on the plausibility of it being possible to define the generating functional via a proper 
limiting procedure. The moments of the Euclidean measure; i.e., "vacuum expectation 
values" of the Euclidean fields, are the n-point Schwinger functions; and the quantum 
field theory is completely determined once all its Schwinger functions are known. The 
time-ordered Green functions of the associated Minkowski space theory can be obtained 
in a formally well-defined fashion from the Schwinger functions [45-48]. This is always at 
the back of the minds of some theorists who are interested in essentially nonperturbative 
QCD. 

However, there is another very important reason to work in Euclidean space; viz., 
owing to asymptotic freedom, all results of perturbation theory are strictly valid only 
at spacelike-momenta. The set of spacelike momenta correspond to a Euclidean vec- 
tor space; and the continuation to Minkowski space rests on many assumptions about 
Schwinger functions that are demonstrably correct only in perturbation theory. It is 
assumed, e.g., that a Wick rotation is valid. Namely, that QCD dynamics do not non- 
perturbatively generate any singularity structure that is not encountered in perturbative 
diagrams. This is a brave assumption, which appears to be completely false in the case 
of coloured states. If the appearance is reality, then QCD must be defined in Euclidean 
space; and the properties of the real-world are then determined only from a continuation 
of colour-singlet quantities. This will pose problems because nothing elementary in QCD 
is a colour singlet. One must somehow solve real-world problems, such as determining 
the spectrum and interactions of complicated two- and three-body bound-states, before 



Strong QCD and Dyson-Schwinger Equations 11 

returning to the real world. That will require imagination and a good toolbox. The DSEs 
provide the latter. 

The Euclidean conventions used herein are easily made plain. For 4-vectors a, b: 

4 

a-b := a^^hv S^^ := ^J Q-i &i , (1-9) 

where S^^ is the Kronecker delta and the metric tensor. Hence, a spacelike vector, Q^, 
has Q^ > 0. The Dirac matrices are Hermitian and defined by the algebra 

{7^.,74 = 2 5^.. (1.10) 

I use 

hij-jliy] and 75 := -71727374 , (l-H) 



t 



a 



/jf 



2' 
so that 

tr [757m7i'7p7ct] = -^e^^pcr , £1234 = 1 • (1-12) 

One may obtain the Euclidean version of any Minkowski space expression by using the 
following transcription rules: 

Configuration Space Momentum Space 



(1)/ « 


'^x^' 


- 




(2) p^ 


^J^■ 


a^ 




(3)4^ 


-i7^ ■ A^ 




(4) A^B> 


^ — )> 


-A^- 


■B^ 


(5) x^d^ 


— )> X' 


b.qe 





(1) / . 


'k" 


-r 


d^k^ 


(2) ^^ 


-i-Y^ ■ k^ 




(3)4^ 


-i-y^ ■ A^ 




(4) fc^g" 


-> - 


-k^-q^ 




(5) k^x" 


-> - 


-k^ ■ x" 





These rules are legitimate in perturbation theory; i.e., the correct Euclidean space integral 
for a given diagram will be obtained by applying them to the Minkowski integral. The 
rules take account of the change of variables and Wick rotation of the contour. When one 
begins with Euclidean space, as I do, the reverse is also true. However, for diagrams that 
represent DSEs which involve dressed n-point functions, whose analytic structure is not 
known a priori, the Minkowski space equation obtained using this prescription will have 
the right appearance but it's solutions may bear no relation to the analytic continuation 
of the solution of the Euclidean equation. Any such differences will be nonperturbative 
in origin. As mentioned above, it is this fact which makes a choice of metric crucial at 
the outset. 



1.4 Quarks and contemporary nuclear physics 

The Standard Model contains six quark fiavours but for the bulk of observable material 
in the universe it is only the u-, d— and s-quarks which are relevant. Whilst there are 
nonetheless good reasons to study systems involving c- and b-quarks, e.g., exploring just 
which states the Standard Model can support [49] and searching for physics beyond the 
Standard Model, I will focus herein on the lighter quarks. 

There is arguably no hadron more interesting than the pion. It was discovered by Cecil 
Powell [50], who was awarded the 1950 Nobel Prize in Physics: "for his development of the 
photographic method of studying nuclear processes and his discoveries regarding mesons 
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Figure 1.4. Left panel - the proton (P) and p-meson each have a mass that is typical of 
hadrons but the pion (vr-meson) is remarkably light, with a mass more typical of leptons. 
Right panel - "eightfold way" pictures of the proton (left) and pion (right), wherein both 
are quantum mechanical objects comprised of constituent-quarks bound by a potential. 



made with this method." The Abstract of Ref. [50] remarks that: "Several features of 
these processes remain to be elucidated, but we present the following account of the 
experiments because the results appear to bear closely on the important problem of 
developing a satisfactory meson theory of nuclear forces." Indeed they did. 

The discovery of the pion was the beginning of modern particle physics. Following 
this came a disentanglement of confusion between the muon and the pion, which have 
similar masses but vastly different interaction patterns: the muon is a lepton and the 
pion, a meson. Then the discovery of strangeness [51]; and subsequently 28 new hadrons 
with masses less-than 1 GeV. Amongst this collection it was abundantly clear that the 
pion is too light: it interacts as a hadron but possesses the mass of a lepton, Fig. 1.4. 

A semblance of order was brought to the hadron spectrum via the "eightfold way" 
[52,53], a picture based on the Lie group SU{3) wherein the u-, d- and s-quarks are repre- 
sented by the basis vectors in the fundamental representation. In such a picture, the pion 
is a quantum mechanical system comprised from a constituent-quark and constituent- 
antiquark, interacting via a potential; and the proton is not vastly different: it is three 
constituent-quarks bound via a symmetry-related potential. In the usual case this would 
mean that the mass of the proton is roughly three-times the mass of the constituents; 
i.e., the constituent-quark mass is Mij = Md ~ Mp/?> ~ 320 MeV. However, the same 
notion produces 77i,r ~ {Mu + Md) ~ 640 MeV, a value that is five-times too large. The 
notion cannot quickly be discarded, however, because the p-meson, which appears to 
be identical to a pion except that the constituent-quark spins are aligned within the p, 
instead of anti-aligned as in the pion, has a mass 770 MeV ~ (Mu + Md)- Furthermore, 
the properties of the vast majority of ground state hadrons appear to be explained by 
such simple counting laws. 

The correct question therefore appears to be: what is wrong with the pion; or how 
does one build an almost massless particle from two massive constituent-quarks. Within 
quantum mechanics one could naturally tune a potential so that the ground-state is 
massless but in that case mbound-stato rises linearly with the sum of the constituents' 
masses. This is not acceptable for the pion because it has long been known that [54] 

m^ (X riniass — term) l-*^-^*^/ 

where Hmass-term is whatever mass-term should be associated with those quarks that ac- 
tually comprise the pion in the Hamiltonian density which describes them. It will become 
clear that in order to resolve the dichotomy of the pion, as a nearly massless bound-state 
of apparently massive constituents, it is necessary to adopt a theory and an approach in 
which one can express: a well-defined and valid chiral limit; and an accurate realisation 
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of dynamical chiral symmetry breaking (DCSB). This is highly nontrivial, impossible in 
quantum mechanics, and plausible only in asymptotically-free gauge theories. 

The flavour index is suppressed in Eq. (1.2), but here I note that M represents a 
current-quark mass matrix: 



M 



rriu 











rrid 











m 



mt 



mt 



(1.14) 



The current-quark masses are external parameters in QCD. They are generated by cou- 
plings between quark currents and the Higgs boson within the Standard Model. With 
six different masses, there are six different couplings. This mechanism produces the mass 
values [55] but doesn't explain them; e.g., no one knows why mt ~ 40 000m„, although 
the number of suggested answers is inversely proportional to this level of certainty. 

Interacting gauge theories, in which it makes sense to speak of massless fermions, 
possess a nonperturbative chiral symmetry. This means that one can set some or all of 
the current-quark masses to zero, and define and preserve the associated set of Ward- 
Takahashi identities. (I will return to this in Sec. 5 because the chiral limit can be a 
complicated concept.) A related concept is helicity, which is the projection of a particle's 
spin, J, onto its direction of motion: 



A = J- 



(1.15) 



For a massless fermion, helicity is a Lorentz-invariant spin-observable: A = ±; i.e., it is 
parallel or antiparallel to the direction of motion. This is readily understood: massless 
particles travel at the speed of light, hence no observer can overtake the particle and 
thereby view its momentum as having changed sign. 

The generator of chiral transformations involves 75 and the simplest chiral rotation 
is 

q{x) ^- e'"'^ q{x) , q{x) ^r q{x) e'"^ . (1-16) 

In connection with a fermion-antifermion bound-state, a rotation through 6 = 7r/4 
changes a parity-positive state into a parity- negative state: J — >■ J^; i.e., a partic- 
ular chiral rotation is equivalent to the operation of parity conjugation (see Sec. 3.2 of 
Ref. [7]). Therefore, if a theory is chirally symmetric, then its spectrum should be ex- 
pected to exhibit degenerate parity partners. (At this point I have overlooked the feature 
of chiral anomalies in gauge field theories because they have no bearing on this argument. 
Section 4.2 canvasses aspects of the non-Abelian anomaly.) 

Two of the current-quarks in perturbative QCD are light [56]: 



« 0.5, rud 



4MeV. 



(1.17) 



Hence the spectrum of states constituted from u- and d-quarks should exhibit almost 
exact parity doubling. However, the splitting between parity partners is typically greater- 
than 100-times the current-quark mass-scale; e.g.. 



state (J^) 
mass (MeV) 



nucleon ( ^ 
940 



i+\ 



nucleon ( | ) 
1535 



(1.18) 
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Plainly, something is happening in QCD: some inherent, dynamical effect is dra- 
matically changing the pattern by which the Lagrangian-density's chiral symmetry is 
expressed. The effect is qualitatively different from spontaneous symmetry breaking, a 
la the Higgs mechanism. Nothing is added to QCD. One has only the original fermion 
and gauge-boson fields, yet the mass operator produces a spectrum with no sign of the 
chiral symmetry which characterised the light-quark sector of the Lagrangian density. 

The key challenges of real-world QCD are now plain. Quarks and gluons are ap- 
parently confined; viz., no matter how hard one strikes the proton, one cannot liberate 
an individual quark or gluon. There is a very unnatural pattern of bound-state masses: 
whilst the Lagrangian (or perturbative) bare masses are very small, the spectrum does 
not exhibit a degeneracy between parity partners. One can express this as a dynami- 
cal breakdown of chiral symmetry. Neither of these phenomena is apparent in QCD's 
Lagrangian density, yet they are the dominant characteristics of real-world QCD. The 
task of hadron physics is to explain the manner through which such complex behaviour 
emerges from rules so simple that they can be expressed in just a few lines; viz., Eqs. (1.3), 
(1.4). 



2 Emergent Phenomena 



Two of the most important hadrons are the neutron and proton, collectively termed the 
nucleons. They are characterised by two static properties: electric change and mag- 
netic moment, hm- The charges are obvious: proton, -fl, and neutron, 0; however, the 
magnetic moments are not quantised. The proton's magnetic moment was discovered 
by Stern in 1933 [57], who received the Nobel Prize for Physics in 1943 [58]: "for his 
contribution to the development of the molecular ray method and his discovery of the 
magnetic moment of the proton." The measured magnetic moment of the proton is 

Mp = (1 + 1.79)^, (2.1) 

where e is the proton's charge and M is its mass. At the time of discovery, this imme- 
diately posed a problem because in Dirac's relativistic quantum mechanics a pointlike 
fermion with charge q and mass m, interacting with an electromagnetic field, has a mag- 
netic moment fi = q/[2m\ [59]. The actual value is almost three-times Dirac's prediction. 

Numerous explanations were proposed but we now know this large discrepancy orig- 
inates from the finite size of the proton. This feature of the proton was not uncovered 
until 1955 [60], a discovery that brought Hofstadter a half-share of the 1961 Nobel Prize 
in physics [61]: "for his pioneering studies of electron scattering in atomic nuclei and for 
his thereby achieved discoveries concerning the structure of the nucleons." (The 1961 
prize also recognised Mossbauer: "for his researches concerning the resonance absorption 
of 7 radiation and his discovery in this connection of the effect which bears his name.") 

The proton is nonpointlike because it is a composite object. The baryon-number 
carrying constituents of the proton are quarks, whose existence was confirmed via high- 
energy electron-proton scattering experiments at the Stanford Linear Accelerator Center 
(SLAC) in 1968. This research resulted in Friedman, Kendall and Taylor being awarded 
the 1990 Nobel Prize in physics [62-65]: "for their pioneering investigations concerning 
deep inelastic scattering of electrons on protons and bound neutrons, which have been 
of essential importance for the development of the quark model in particle physics." 
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The elastic and inelastic scattering of electrons from an hadronic target are still key 
tools in hadron physics. The electron in particular, and leptons in general, are excellent 
probes because, on the hadron scale at least, they are structureless. Their interactions 
are described by a simple current 

j;iP\P) = -ieMP')A,{P, Q)u.{P) -PP-"F-«°" -ieue(P')7M^e(P) , (2.2) 

where the momentum transfer Q — P' — P. Consequently, structures revealed by such 
experiments can typically be interpreted as properties of the hadron target. 
The nucleon's electromagnetic current is 

j^{P',P) = ieuN{P')A^{P,Q)uN{P) (2.3a) 



ieuM{P ) 



1^F^n[Q') + j^jO^,Q,F2n{Q') 



un{P). (2.3b) 



Here, -Fijv((5^) and F2n{Q^) are, respectively, the nucleons' Dirac and Pauli form factors: 
Fi is called the Dirac form factor because of its association with the matrix structure of 
the current for a pointlike fermion; and F2 is the Pauli form factor because it is associated 
with an additional ~ or anomalous Pauli [66] - contribution to the particle's magnetic 
moment. The Sachs electric and magnetic form factors are expressed in terms of _Fi,2 as 
follows: 

Gen{Q^) = F^m{Q^)-^F2n{Q'), (2.4a) 

Gmn{Q^) = Fin{Q^) + F2n{Q^). (2.4b) 

If one is studying a system for which a nonrelativistic limit is a useful or valid approx- 
imation, then Fourier transforms of the Sachs form factors produce configuration-space 
densities that describe the distribution of charge and magnetisation within the composite 
object. A comparison between Eqs. (2.2) and (2.3b) reveals that for a structureless or 
simply-structured fermion F\{Q^) = 1 and F2{Q^) = 0, and hence Ge{Q^) = Gm{Q^)- 
Thus the distribution of charge and magnetisation is identical within composite but 
simply-structured fermions. 

As observed in Sec. 1.2, a central goal of nuclear physics is to understand the struc- 
ture and properties of protons and neutrons, and ultimately atomic nuclei, in terms of 
the quarks and gluons of QCD. Confinement and dynamical chiral symmetry breaking, 
in fact, a fundamental theory of unprecedented complexity means that this aim is far 
from being achieved. The struggle with QCD defines the difference between nuclear and 
particle physics: nuclear physicists try to solve this theory, whereas particle physicists 
escape to a place where perturbation theory, the last refuge of a scoundrel, is all that is 
necessary. 

Before one can make progress toward this goal, one must first grasp what is meant 
by "the quarks and gluons of QCD." Is there such a thing as a constituent-gluon or a 
constituent-quark? After all, these are the concepts for which Gell-Mann was awarded 
the Nobel Prize. Do they or can they be rigorously defined in QCD, perhaps as quasi- 
particle degrees-of- freedom? If not, then with what should they be replaced and how 
should the goal be rigorously interpreted? Finally, why is it important to achieve this 
goal? 

The last question is readily answered. Recall the dichotomy of the pion, discussed 
in connection with Eq. (1.13): how does one compose an almost massless object from 
two very massive constituents? The alternative is to have m-n ~ rUp. In that case the 
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attractive and repulsive forces in the nucleon-nucleon potential have the same range and 
there is no intermediate range attraction. Under these circumstances, it is probable 
that [67]: the deuteron would not be stable, big-bang nucleosynthesis would not have 
occurred; ^^C could not have been produced and wouldn't be stable if it were; etc. These 
questions matter because their answers explain how we came to exist and ask them. 

Neither confinement nor DCSB is apparent in QCD's Lagrangian and yet they play 
the dominant role in determining the observable characteristics of real-world QCD. The 
physics of hadrons is ruled by emergent phenomena such as these, which can only be 
elucidated through the use of nonperturbative methods in quantum field theory. This is 
both the greatest novelty and the greatest challenge within the Standard Model. We must 
find essentially new ways and means to explain precisely via mathematics the observable 
content of QCD. 



2.1 Dyson- Schwinger equations 

Beginning with the quark model, many insights into the nature of hadrons have been ob- 
tained by using simple models of their structure. A problem with such studies, however, 
is that they typically employ degrees-of-freedom, potentials, etc., that have no ascertain- 
able connection with QCD. At another extreme is the brute-force numerical simulation 
of lattice-regularised QCD. This is today a sizeable international endeavour, with the 
annual five-day lattice conference attracting approximately four-hundred participants. 
There are numerous reviews, recent amongst them being Refs. [68-71]. 

In between, however, there is a compromise. In 1994 the following observations were 
made [3]: "As computer technology continues to improve, lattice gauge theory [LGT] 
will become an increasingly useful means of studying hadronic physics through investiga- 
tions of discretised quantum chromodynamics [QCD] . . . However, it is equally important 
to develop other complementary nonperturbative methods based on continuum descrip- 
tions. In particular, with the advent of new accelerators such as CEBAF (VA) and 
RHIC (NY), there is a need for the development of approximation techniques and mod- 
els which bridge the gap between short-distance, perturbative QCD and the extensive 
amount of low- and intermediate-energy phenomenology in a single covariant framework 
. . . Cross-fertilisation between LGT studies and continuum techniques provides a par- 
ticularly useful means of developing a detailed understanding of nonperturbative QCD." 
These statements have proved true and a useful compromise is provided by QCD's Dyson- 
Schwinger equations (DSEs). 

It has long been known [72-74] that one can derive a system of coupled integral 
equations relating all the Green functions for a theory, one to another. The gap equation, 
illustrated in Fig. 2.1, is one such expression. The DSEs are nonperturbative equivalents 
in quantum field theory of the Lagrange equations of motion and, amongst their many 
uses, they are essential in simplifying the general proof of renormalisability of gauge field 
theories. 

The DSEs are a particularly well suited for use in QCD. At the simplest level they 
provide a generating tool for perturbation theory. It therefore follows from asymptotic 
freedom that in all applications any model-dependence is restricted to the infrared or, 
equivalently, the long-range domain. On the other hand, a nonperturbative solution of 
the DSEs enables the study of: hadrons as composites of dressed-quarks and -gluons; the 
phenomena of confinement and DCSB; and therefrom an articulation of any connection 
between them. The solutions of the DSEs are Schwinger functions and because all cross- 
sections can be constructed from such n-point functions, the DSEs can be used to make 
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Figure 2.1. In QCD this gap equation describes the manner by which a quark's acquires 
a momentum-dependent self-energy, E(p), via interaction with its own gluon field. This 
diagram represents a nonlinear integral equation because S{p) — l/[i'y-p+m + T,{p)]. The 
kernel of that integral equation is composed from the dressed-gluon propagator, D^i/(fc), 
and the dressed-quark-gluon vertex, r^{£,p), both of which satisfy integral equations of 
their own. (The quark is represented by the solid line, the gluon by the spring-like line 
and the quark-gluon vertex by the open circle.) 



predictions for real-world experiments. One of the merits in this is that assumptions 
about the infrared form of n-point functions can be tested, verified and improved, or 
rejected in favour of more promising alternatives. In this mode the DSEs provide a bridge 
between experiment and theory, and thereby a means by which to use nonperturbative 
phenomena to chart or at least constrain the infrared behaviour of QCD's /3-function. 

Let's return to the dressed-quark propagator, which is given by the solution of QCD's 
gap equation. Fig. 2.1: 

S{py^ = i7-p -f m-f S(p), (2.5) 

^iP) - / (gl ^^ ^^^iP - ^) ^-T ., . lAiP) + BiP) ^"(^'^)- (2.6) 

The most general Poincare covariant solution of this gap equation involves two scalar 
functions. There are three common, equivalent expressions: 

^(^) ^ i^-pA{p^) + B{p^) ^ ^J■p + Mip^) ^ -'^ • ^"^'(^') + "^(^'^ ■ (2.7) 

In the second form, Z{p^) is called the wave-function renormalisation and M(p'^) is the 
dressed-quark mass function. 

Given that a weak coupling expansion of the DSEs produces every diagram in pertur- 
bation theory, one can use Eq. (2.6), and its analogues for the dressed-gluon propagator 
and dressed-quark gluon vertex, to obtain 



Bo 



.(/)^^[l + ^ln[^] +...], (2.8) 



where the ellipsis denotes terms of higher order in a that involve (ln[rn,^/p^])^ and 
i\n\a[m? /p^]), etc. However, at arbitrarily large finite order in perturbation theory it is 
always true that 

lim Sport (p^) = 0. (2.9) 

(N.B. In contrast to QED, a chiral limit is nonperturbatively defined in QCD owing to 
asymptotic freedom.) Equation (2.9) means that if one starts with a chirally symmetric 
theory, then in perturbation theory one also finishes with a chirally symmetric theory: 
the fermion DSE cannot generate a mass-gap if there is no bare-mass seed in the first 
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place. Thus DCSB is impossible in perturbation theory. The question is whether this 
conclusion can ever be avoided; namely, are there circumstances under which it is possible 
to obtain a nonzero dressed-quark mass function in the chiral limit? 



2.2 Exploring dynamical mass generation 

To begin the search for an answer, recall the gap equation, Eq. (2.6), and consider the 
expression obtained with the following forms for the dressed-gluon propagator and quark- 
gluon vertex 

5^D^.(p-^) = 5^,^, r;:(fc,p) = 7,^, (2.10) 

IJIq 1 

wherein mc ~ 0.8 GeV is a mass-scale, typical of gluons, as we shall see, and aia is 
a parameter that specifies the interaction strength. (One may view the interaction in 
Eq. (2.10) as being within a class inspired by the Nambu-Jona-Lasinio model [75]. The 
application to which it is here put most certainly relates to Nambu's half-share in the 
2008 Nobel Prize in physics [76]: "for the discovery of the mechanism of spontaneous 
broken symmetry in subatomic physics.") Using Eqs. (2.10), the gap equation becomes 

S \p)=i^.p + m+—--^ j j^-i^S{q)^^, (2.11) 

an expression in which the integral possesses a quadratic divergence, even in the chiral 
limit. When the divergence is regularised in a Poincare covariant manner, the solution 

is [77] 

^(p) = -■ — \rM ' (2.12) 

i.e., A{p'^) = 1 and M{p^) = M, a momentum-independent quantity determined by the 
following one-dimensional nonlinear integral equation: 

M = m + A/i^ fds s \— . (2.13) 

3nm% J s + AP ^ ' 

(One arrives at these results via a little matrix algebra involving the Dirac matrices and 
capitalising on 0(4) invariance.) 

A heat-kernel-like regularisation procedure is useful in hadron physics applications; 
viz., one writes [78]: 



roo fT 



dre 



-T{a + M^) 



-(s + M^)t-^ 



s + M^ Jo J ^2 s + M^ 

(2.14) 
where rir.uv are, respectively, infrared and ultraviolet regulators. Since Eq. (2.10) does 
not define a renormalisable theory, then Auv :— I/t^v cannot be removed but instead 
plays a dynamical role, setting the scale of all dimensioned quantities. Using Eq. (2.14), 
the gap equation becomes 






M = m-HM:j--fi-C'"(M"), (2.15) 



where C''{l\/P)/M^ = r{-l,M^T^^) - r(-l, AfVj?), with r{a,y) being the incomplete 
gamma- function. 
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In the current application, one may set Air = 0. However, when employed in hadron 
physics phenomenology, a finite value of rir =: l/Air implements confinement by ensuring 
the absence of quark production thresholds in ^-matrix amplitudes [79] . This is apparent 
from Eq. (2.14). Subsequently, Auv defines the only mass-scale in a nonrenormalisable 
model. Hence one can set Auv = 1 and hereafter merely interpret all other mass-scales 
as being expressed in units of Auv, in which case C°^{M^) — M'^r(— 1, M^). 

Now consider Eq. (2.15) in the chiral limit; i.e., mo — 0: 

M = Mi^ C°'(Af'), (2.16) 

37rm^ 

where rhc = mo/ Auv One solution of Eq. (2.16) is obviously M = 0. This is the 
result that connects smoothly with perturbation theory: one begins with no mass, and 
no mass is generated. In this instance the theory is said to realise chiral symmetry in the 
Wigner-Weyl mode. 

Suppose, on the other hand, that M 7^ in Eq. (2.16). The expression then makes 
sense if, and only if, the following equation has a solution: 



4qir ^01/ 
Svrm? 



1 = TZrFC"'(M^). (2.17) 



C^^{M'^) is a monotonically decreasing function with maximum value one at M = 0. 
Consequently, 3Af 7^ solution if, and only if, 

>-mG. (2.18) 

TV 4 

It is thus apparent that for any gluon mass-scale, mo, there is always a domain of values 
for QfiR upon which a nontrivial solution of the gap equation can be found. If one supposes 
Auv ~ rnc , then 3M 7^ solution for 

^>\. (2.19) 

■K 4 

This result, derived in a straightforward manner, is remarkable enough to have de- 
livered a Nobel Prize to Nambu. It reveals the power of a nonperturbative solution to 
nonlinear equations. Although one began with a model of massless fermions, the interac- 
tion alone has provided those fermions with mass. This is the phenomenon of dynamical 
chiral symmetry breaking; namely, the generation of mass from nothing; and when it hap- 
pens, chiral symmetry is said to be realised in the Nambu-Goldstone mode. It is clear 
from Eqs. (2.18), (2.19) that DCSB is guaranteed to be possible so long as the interaction 
exceeds a particular minimal strength. This behaviour is typical of gap equations. 

Figure 2.2 depicts the QiR-dependence of the self-consistent solution to Eq. (2.15) 
obtained with a current-quark mass m — 0.01 (measured in units of Auv). The vertical 
line marks the critical coupling in Eq. (2.19); namely, the value Oj^ — 3/4 below which 
a DCSB solution of Eq. (2.16) is impossible. Evidently, for qir < Qjr, the dressed mass 
M « mo; i.e., the self-consistent solution is well approximated by the perturbative result. 
However, a transition takes place for a ~ a" , and for oib. > am the dynamical mass 
is much greater than the bare mass, with M increasing rapidly as the coupling qir is 
increased. 



2.3 Dynamical mass and confinement 

One aspect of quark confinement is the absence from the strong interaction spectrum 
of free-particle-like quarks. What does the model of Sec. 2.2 have to contribute in this 
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Figure 2.2. Dynamical mass obtained from Eq. (2.15) as a function of the coupling, q/tt, 
with bare mass 0.01, marked by the horizontal dashed line, Air = and Auv = ma- The 
vertical dotted line marks the critical coupling in Eq. (2.19). 



connection? Well, whether one works within the domain of the model on which DCSB 
takes place, or not, the quark propagator always has the form, Eq. (2.12), 



5^-'^(p) 



i-y- p[A{p^) ^ 1] + [B{p^) ^ M] 



-i'y ■ p -\- M 



(2.20) 



where M is constant. This expression has a pole at p^ + M^ = and thus is always 
effectively the propagator for a noninteracting quasiparticle fermion with mass M: the 
fermion propagates as a plane wave. Hence, while it is true that models of the Nambu- 
Jona-Lasinio type support DCSB, they do not generally exhibit confinement. (This is 
true so long as Air = 0.) 

Consider an alternative theory associated with a gap equation defined via the Munczek- 
Nemirovsky interaction [80] 



g^D^,{k) = {2i.fG5\k)[5^,. _ :^] , r,(Ap) = 7. 



(2.21) 



Here G defines the model's mass-scale and the interaction is a i5-function in momentum 
space, which is the antithesis to models of the Nambu-Jona-Lasinio type wherein the 
interaction is instead a 5- function in configuration space. Using Eq. (2.21), the gap 
equation is 



i-i-pA{p) + B{p ) = »7-p + mo + G7M ^242(^2)^^2(^2) T^ , 

which yields the following coupled nonlinear algebraic equations: 

A{p^) 



A{p') = 1 + 2 



B{p^) = mo + 4 



p242(p2)+B2(p2) 

B{p^) 



p2yl2(p2)+B2(p2) 



(2.22) 

(2.23) 
(2.24) 
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Equation (2.22) yields an ultraviolet finite model and hence there is no regularisation 
mass-scale. In this instance one can therefore refer all dimensioned quantities to the 
model's mass-scale and set G = 1. 

Consider the chiral limit of Eq. (2.24); 

^(^^) = %M^(A%(pr ^'-''^ 

Obviously, like Eq. (2.16), this equation admits a trivial solution B{p^) = that is 
smoothly connected to the perturbative result, but is there another? The existence of 
a, B ^ solution; i.e., a solution that dynamically breaks chiral symmetry, requires (in 
units of G) 

p'A"(p') + B'(p') = 4. (2.26) 

Suppose this identity to be satisfied, then its substitution into Eq. (2.23) gives 

A(/) - 1 = i Aip^) => Aip^) = 2 , (2.27) 

which in turn entails 



B(p') = 2\/l-p2. (2.28) 

A complete chiral-limit solution is composed subject to the physical requirement that 
the quark self-energy is real at spacelike momenta, and hence 

^^^'^ - {j^(i + v/TTW); P^>i ^'-''^ 

B{P') = 1^^^^^' <z\ (2.30) 



; p-' > 1 . 

In this case both scalar functions characterising the dressed-quark propagator differ sig- 
nificantly from their free-particle forms and are momentum dependent. (We will see that 
this is also true in QCD.) It is noteworthy that the magnitude of the model's mass-scale 
plays no role in the appearance of this DCSB solution of the gap equation; viz., rein- 
stating G, in the chiral limit one has, VG 7^ 0, -B(O) = 2G 7^ 0. Thus, in models of 
the Munczek-Nemirovsky type, the interaction is always strong enough to support the 
generation of mass from nothing. 

The DCSB solution of Eqs. (2.29), (2.30) is defined and continuous for all values of 
p^, including timelike momenta, p^ < 0. It gives a dressed-quark propagator whose 
denominator 

p'^"(p") + B'(p')>0, Vp". (2.31) 

This is an extraordinary result. It means that the propagator does not exhibit any free- 
particle-like poles. It has long been argued that this feature can be interpreted as a 
realisation of quark confinement; e.g., Refs. [79-83]. 

The interaction in Eq. (2.21) has taken a massless quark and turned it into something 
which at timelike momenta bears little resemblance to the perturbative quark. It does 
that for all nonzero values of the model's mass-scale. This model exemplifies an intriguing 
possibility; namely, that all models with quark confinement necessarily exhibit DCSB. It 
is obvious from Sec. 2.2 that the converse is certainly not true. 

In the chirally asymmetric case the gap equation yields 
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M(5) Munczek-Nemirovsky 
M(s) = 0.015 
M(s) = |s|'", s<0 



Figure 2.3. With a bare mass m — 0.015, the dressed- quark mass function M{s = p^) 
obtained using the interaction specified by Eq. (2.21) - solid curve. A free-particle-like 
mass-pole would occur at that value of s for which the solid and green curves intersect. 
As the figure suggests, this never happens in models of the Munczek-Nemirovsky type. 
(All dimensioned quantities measured in units of G, the model's mass-scale.) 



The second of these is a quartic equation for B(p^). It can be solved algebraically. There 
are four solutions, obtainable in closed form, only one of which possesses the physically 
sensible ultraviolet spacelike behaviour: B{p^) — ^ m as p^ — > cxd. The physical solution is 
depicted in Fig. 2.3. At large spacelike momenta, M{s = p^) — >■ m^ and a perturbative 
analysis is reliable. That is never the case for s < 1 (in units of G), on which domain 
M[s) >• m, a defining characteristic of DCSB; and the difference {M{s) — \s\) is always 
nonzero, a feature that is consistent with confinement. 

I have illustrated two variations on the theme of dynamical mass generation via the 
gap equation. A general class of models for asymptotically free theories may be discussed 
in terms of an effective interaction 



g'D{Q^)^47T 



a(Q' 



(2.34) 



where a{Q'^) is such that g^D{Q^) evolves according to QCD's renormalisation group in 
the ultraviolet. This situation is depicted in Fig. 2.4. The types of effective interaction 
fall into two classes. In those for which a(0) < 1 [recall Eq. (2.19)], the only solution of 
the gap equation in the chiral limit is M{s) = 0. Whereas, when a(0) > 1, M{s) 7^ is 
possible in the chiral limit and, indeed, corresponds to the energetically favoured ground 
state [3,84-86]. In this case M(0) 7^ is a new, dynamically generated mass-scale. If 
large enough, it can explain how a theory that is perturbatively massless nevertheless 
possesses the spectrum of a massive theory. 

It should now be plain that confinement and DCSB are key emergent phenomena 
in QCD, and that their understanding requires the nonperturbative solution of a fully- 
fledged relativistic quantum field theory. Both mathematics and physics are still far 
from being able to accomplish that. On the other hand, I have illustrated that confine- 
ment and DCSB are expressed in QCD's propagators. These phenomena are also evident 
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Figure 2.4. Two classes of effective interaction in the gap equation. A theory will exhibit 
dynamical chiral symmetry breaking if, and only if, a{Q = 0) > 1. Here the ultraviolet 
behaviour of the interactions is constrained to be that of QCD. 



in dressed vertices. Such nonperturbative modifications should have observable conse- 
quences. The DSEs are a useful analytical and numerical tool for the nonperturbative 
study of relativistic quantum field theory. They can be used to show that simple models 
exhibit DCSB but that DCSB 7^ confinement; and conversely, that simple models can 
exhibit confinement and that confinement => DCSB. However, the question of what can 
be established for QCD remains. 



2.4 Confinement 



Gluon and quark confinement is an empirical fact: no matter how hard one has thus 
far struck the proton, or any other hadron, no individual gluon or quark has yet been 
liberated. However, there is no agreed, theoretical definition of light-quark confinement. 
One thing is nevertheless certain; namely, true confinement entails quark-hadron duality, 
which is the statement that all observable consequences of a confining QCD could, in 
principle, be computed using an hadronic basis. 

Confinement is often discussed in terms of Wilson's area law, exposed in an influential 
article entitled "Confinement of quarks," which introduced the concept of lattice gauge 
theory [37]. The area law identifies confinement with a linearly rising potential between 
coloured charges. Phenomenologists with an attachment to quantum mechanics and 
string theorists, too, have long been nourished by this result. Unfortunately, the only 
quarks in the article are static; i.e., infinitely heavy, so the concept doesn't bear on real- 
world QCD. It is therefore worth stating plainly that the potential between infinitely- 
heavy quarks measured in numerical simulations of quenched lattice-regularised QCD - 
the so-called static potential - is simply irrelevant to the question of confinement in the 
real world, in which light quarks are ubiquitous. In fact, it is a basic feature of QCD 



24 Craig D. Roberts 

Normal particle Confined particle 

complex-P^ complex-P^ 



Figure 2.5. Left panel - An observable particle is associated with a pole at timelike-P^. 
This becomes a branch point if, e.g., the particle is dressed by photons. Right panel - 
When the dressing interaction is confining, the real-axis mass-pole splits, moving into 
pairs of complex conjugate poles or branch points. No mass-shell can be associated with 
a particle whose propagator exhibits such singularity structure. 



that light-particle creation and annihilation effects are essentially nonperturbative and 
therefore it is impossible in principle to compute a potential between two light quarks 
[87]. 

Some consider that a belief in linear potentials and their relevance is justified by the 
appearance of Regge trajectories in the hadron spectrum; i.e., that associated with a 
given ground state composite, there are higher mass composite states with larger total 
angular momentum (spin J), all of which lie on a straight line in the plane {J,M'^). 
The appearance does not necessarily match reality, however, as observed in Ref. [88] : "In 
time the Regge trajectories thus became the cradle of string theory. Nowadays the Regge 
trajectories have largely disappeared, not in the least because these higher spin bound 
states are hard to find experimentally. At the peak of the Regge fashion (around 1970) 
theoretical physics produced many papers containing families of Regge trajectories, with 
the various (hypothetically straight) lines based on one or two points only!" A detailed 
analysis supports this view; viz., Ref. [89]: "Early Chew-Frautschi plots show that meson 
and baryon Regge trajectories are approximately linear and non-intersecting. In this 
paper, we reconstruct all Regge trajectories from the most recent data. Our plots show 
that meson trajectories are non-linear and intersecting. We also show that all current 
meson Regge trajectories models are ruled out by data." 

Drawing on a long list of sources; e.g., Refs. [80-83], a perspective on confinement 
was laid out in Ref. [79]. As illustrated above and in Fig. 2.5, confinement can be related 
to the analytic properties of QCD's Schwinger functions. For example, it can be read 
from the reconstruction theorem [46,47] that the only Schwinger functions which can be 
associated with expectation values in the Hilbert space of observables - namely, the set of 
measurable expectation values - are those that satisfy the axiom of refiection positivity. 
This is an extremely tight constraint. It can be shown to require as a necessary condition 
that the Fourier transform of the momentum-space Schwinger function is a positive- 
definite function of its arguiuents. This condition suggests a practical confinement test, 
which can be used with numerical solutions of the DSEs (see, e.g.. Sec. III.C of Ref. [90] 
and Sec. IV of Ref. [91]). The implications and use of refiection positivity are discussed 
and illustrated in Sec. 2 of Ref. [10]. 

It is notable that any 2-point Schwinger function with an infiexion point at p^ > 
must breach the axiom of refiection positivity, so that a violation of positivity can 
be determined by inspection of the pointwise behaviour of the Schwinger function in 
momentum space (Sec. IV. B of Ref. [92]). 



in 
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Figure 2.6. Dyson-Schwinger equation for the gluon propagator. The gluon vacuum 
polarisation tensor n,i^ = (Sf^i, — k^ky /k )n(fc ) is the gluon's self energy, from which 
the propagator is obtained: in Landau gauge, D^^ik) = {S^^ — k^k^ /k^)\l/ {k'^[l+Yi{k^)\)\ 
is the propagator obtained (The broken line represents the propagator for the ghost field.) 



The gluon two-point function satisfies the gap equation depicted in Fig. 2.6. Consider 
then A(fc ) — l/(fc [1 -|- n(fc )]), which is the single scalar function that describes the 
dressing of a Landau-gauge gluon propagator. A large body of work has focused on 
exposing the behaviour of A(fe^) in the pure Yang- Mills sector of QCD. These studies 
are reviewed in Ref. [15]. A connection with the expression and nature of confinement 
in the Yang-Mills sector is indicated in Fig. 2.7. The appearance of an inflexion point in 
the two-point function generated by the gluon's momentum-dependent mass-function is 
impossible to overlook. Hence this gluon cannot appear in the Hilbert space of observable 
states. At present, DSE and lattice-QCD agree that the gluon is truly described by a 
momentum-dependent mass function but the infrared scale is only known to within a 
factor of two; i.e., 

mc — lim fc^n(fc^) ^2-4Aqcd. (2.35) 

From the perspective that confinement can be related to the analytic properties of 
QCD's Schwinger functions, the question of light-quark confinement can be translated 
into the challenge of charting the infrared behavior of QCD's universal /3-function. (Al- 
though this function may depend on the scheme chosen to renormalise the theory, it is 
unique within a given scheme [93]. Of course, the behaviour of the /?- function on the 
perturbative domain is well known.) This is a well-posed problem whose solution is an 
elemental goal of modern hadron physics and which can be addressed in any framework 
enabling the nonperturbative evaluation of renormalisation constants. It is the /3- function 
that is responsible for the behaviour evident in Fig. 2.7, and Fig. 2.9 below; and one of 
the more interesting of contemporary questions is whether it is possible to reconstruct 
the /3-function, or at least constrain it tightly, given empirical information on the gluon 
and quark mass functions. 

In this connection I note that the DSEs connect the /3-function to experimental ob- 
servables. Hence, comparison between DSE computations and observations of the hadron 
mass spectrum, and elastic and transition form factors, can be used to chart /3-function's 
long-range behaviour. Extant studies show that the properties of hadron excited states 
are a great deal more sensitive to the long-range behaviour of the /3-function than those 
of ground states. This is illustrated in Refs. [94,95], which through a study of ground- 
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Figure 2.7. A(fc^), the function that describes dressing of a Landau-gauge gluon propa- 
gator, plotted for three distinct cases. A bare gluon is described by A{k^) = l/fc^ (the 
dashed line), which is plainly convex on k^ £ (0, oo). Such a propagator has a representa- 
tion in terms of a non-negative spectral density. In some theories, interactions generate a 
mass in the transverse part of the gauge-boson propagator, so that A(fc^) = l/(fc^ -t-rrig), 
which can also be represented in terms of a non-negative spectral density. In QCD, how- 
ever, self-interactions generate a momentum-dependent mass for the gluon, which is large 
at infrared momenta but vanishes in the ultraviolet [15]. This is illustrated by the curve 
labelled "IR-massive but UV-massless." With the generation of a mass-function, A{k^) 
exhibits an inflexion point and hence cannot be expressed in terms of a non-negative 
spectral density. 



state, radially-excited and exotic scalar-, vector- and flavoured-pseudoscalar-mesons in 
rainbow-ladder truncation, which is leading order in the most widely used nonperturba- 
tive scheme [96,97], produced the effective coupling and running gluon mass depicted in 
Fig. 2.8. 



2.5 Apergu 

Whilst confinement is empirically true. Fig. 2.9 shows that dynamical chiral symme- 
try breaking is a fact in QCD. Contemporary DSE- and lattice-QCD studies agree on 
the basic features of the mass function for a dressed light-quark. The mass function 
is momentum dependent. For light quarks with nonzero Lagrangian mass it exhibits 
logarithmic growth as one proceeds toward the infrared from large momenta: this is 
the perturbative domain upon which asymptotic freedom is evident. At a momentum 
scale of approximately 5Aqcd, however, the rate of increase accelerates dramatically as 
power-law growth sets in: this is the beginning of the nonperturbative domain, where- 
upon dynamical chiral symmetry breaking is key to the evolution of a dressed- quark's 
mass. As the m = curve in Fig. 2.9 shows, even a quark that is perturbatively massless 
would acquire a mass in QCD. This is kindred to the behaviour of the running gluon 
mass in Fig. 2.8. The inflexion point at p ~ Aqcd, which is predicted by DSE studies and 
appears to be supported by the lattice results, indicates that beyond this point one has 
very probably entered the confinement domain. The inflexion point is connected with 
the behaviour illustrated in the right panel of Fig. 2.5. 
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0.6 0.8 



Figure 2.8. Left panel - Rainbow-ladder gluon running-mass; and nght panel - rainbow- 
ladder effective running-coupling. Matching related curves in each panel (solid to solid 
or dashed to dashed) gives an equivalent description of observables within the rainbow- 
ladder truncation. 



The reality of DCSB means the Higgs mechanism is largely irrelevant to the bulk of 
normal matter in the Universe. Instead the most important mass generating mechanism 
for light-quark hadrons is the strong interaction effect of DCSB; e.g., one may identify it 
as being responsible for 98% of a proton's mass [102,103]. There is a caveat; namely, as so 
often, the pion is exceptional. Its mass is given by the simple product of two terms, one 
of which is the ratio of two order parameters for DCSB whilst the other is determined by 
the current-quark mass (Sec. 4.1). Hence the pion would be massless in the absence of a 
mechanism that can generate a current-mass for at least one light-quark. The impact of a 
massless, strongly-interacting particle on the physics of the Universe would be dramatic. 

It is QCD's /^-function that is responsible for the behavior evident in Figs. 2.7 and 
2.9, and thereby the scale-dependence of the structure and interactions of dressed-giuons 
and -quarks. Elucidation of an unambiguous connection between observables and the 
behaviour evident in Figs. 2.7 and 2.9 is therefore a promising method for empirically 
charting the long range behaviour of the strong interaction, something which will not 
be achieved by theory alone in the foreseeable future. As described and illustrated in 
Sees. 8, 9 below, this can potentially lead to an explanation and understanding of the 
origin of the interaction strength at infrared momenta that guarantees DCSB through 
the gap equation and ties confinement to DCSB. 



3 Gap and Bethe-Salpeter Equations 



In order to proceed it is useful to recapitulate on the best known and simplest DSE. The 
Dyson or gap equation describes how quark propagation is influenced by interactions; 
viz., for a quark of flavour /, 



Sf{p)-' = Z2ii^-p + m)n + Z^j g^D^,.{p - q)^j^Sf{q)^rUq, 



(3.1) 



where: D^^ is the gluon propagator; F^, the quark-gluon vertex; J , a symbol that 
represents a Poincare invariant regularization of the four-dimensional Euclidean inte- 
gral, with A the regularization mass-scale; 7n^™(A), the current-quark bare mass; and 
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Figure 2.9. Dressed-quark mass function, M(p): solid curves - DSE results, [98-100], 
"data" - lattice-QCD simulations [101]. (N.B. m = 70MeV is the uppermost curve. 
Current-quark mass decreases from top to bottom.) The constituent mass arises from a 
cloud of low-momentum gluons attaching themselves to the current-quark: DCSB is a 
truly nonperturbative effect that generates a quark mass from nothing; namely, it occurs 
even in the chiral limit, as evidenced by the tti = curve. 



^i,2(C'^) A'^), respectively, the vertex and quark wave-function renormalisation constants, 
with C, the renormalisation point - dependence upon which is not usually made explicit. 
The gap equation will be familiar from Eq. (2.5) but here I've been more careful with the 
definitions. 

The gap equation's solution is the dressed-quark propagator, 

1 _ z{p\e) 



S{P) 



i-fpA(p'^,C)+B{p'^,C2) i-fp + M{p'^) ' 



(3.2) 



which is obtained from Eq. (3.1) augmented by a renormalisation condition. A mass- 
independent scheme is a useful choice and can be implemented by fixing all renormalisa- 
tion constants in the chiral limit. (See, e.g., Ref. [104] and references therein; or Ref. [105] 
for a detailed discussion of renormalisation.) 

The mass function, M(p^) = _B(p^, ("^)/j4(p^, (^^), is independent of the renormalisa- 
tion point, ("; and the renormalised current-quark mass. 



.(C,A)m'^-(A) 



^r'^: 



bm 



(3.3) 



wherein Z4 is the renormalisation constant associated with the Lagrangian's mass-term. 
Like the running coupling constant, this "running mass" is familiar from textbooks. 
However, it is not commonly appreciated that m'' is simply the dressed-quark mass 
function evaluated at one particular deep spacelike point; viz. 



mj=Mf{e). 
The renormalisation-group invariant current-quark mass may be inferred via 



rhf — lim 

p^ — ^oc 



— In . „ 

2 A2 



P 



QCD 



Mfip' 



(3.4) 



(3.5) 
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where 7^ = 12/(33 — '2Nf) (see Fig. 1.3). The chiral hmit is expressed by 

m/=0. (3.6) 



Moreover, 



777- 

VC»Aqcd, -^='^^. (3.7) 



However, I would like to emphasise that in the presence of DCSB the ratio 



^C=P= Af^,(p2 



(3.8) 



is not independent of p^: in the infrared; i.e., Vp^ < Aqqq, it then expresses a ratio of 
constituent-like quark masses, which, for light quarks, are two orders-of-magnitude larger 
than their current-masses and nonlinearly related to them [102,103]. 

The gap equation illustrates the features and fiaws of each DSE. It is a nonlinear 
integral equation for the dressed-quark propagator and hence can yield much-needed 
nonperturbative information. However, the kernel involves the two-point function D^^ 
and the three-point function Vl. The gap equation is therefore coupled to the DSEs 
satisfied by these functions, which in turn involve higher n-point functions. Hence the 
DSEs are a tower of coupled integral equations, with a tractable problem obtained only 
once a truncation scheme is specified. It is unsurprising that the best known truncation 
scheme is the weak coupling expansion, which reproduces every diagram in perturbation 
theory. This scheme is systematic and valuable in the analysis of large momentum trans- 
fer phenomena because QCD is asymptotically free but it precludes any possibility of 
obtaining nonperturbative information. 

Given the importance of DCSB in QCD, it is significant that the dressed-quark prop- 
agator features in the axial-vector Ward-Takahashi identity, which expresses chiral sym- 
metry and its breaking pattern: 

P^ri^^{k;P)+i[mf{0+m,{0]ri'{k;P)^Sj\k+)iro+i-/sS;\k.), (3.9) 

where P = pi + p2 is the total-momentum entering the vertex and k is the relative- 
momentum between the amputated quark legs. To be explicit, fc = (1 — 77)pi + rjp2, 
with 77 G [0, 1], and hence k+ = p\ = k -\- r/P, fc_ = p2 = fe — (1 — rj)P. In a Poincare 
covariant approach, such as presented by a proper use of DSEs, no observable can depend 
on 77; i.e., the definition of the relative momentum. N.B. Sec. 4.2 discusses the important 
differences encountered in treating fiavourless pseudoscalar mesons. 

In Eq. (3.9), Fj^ and F5 are, respectively, the amputated axial-vector and pseu- 
doscalar vertices. They are both obtained from an inhomogeneous Bethe-Salpeter equa- 
tion (BSE), which is exemplified here using a textbook expression [106]: 

[r,^ik■P)]tu = Z2h5-(^.]tu+ f [S{q+)rs^iq;P)Siq-)]srKi:{q,k-P), (3.10) 

Jq 

in which K is the fully-amputated quark-antiquark scattering kernel, and the colour-, 
Dirac- and flavour-matrix structure of the elements in the equation is denoted by the 
indices r, s, t, u. N.B. By definition, K does not contain quark-antiquark to single gauge- 
boson annihilation diagrams, nor diagrams that become disconnected by cutting one 
quark and one antiquark line. 
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Figure 3.1. Illustration of the procedure used to produce a nonperturbative systematic 
symmetry-preserving truncation. It produces a modified skeleton expansion, in which 
the propagators are fully-dressed but the vertices are constructed term-by-term. 



The Ward-Takahashi identity, Eq. (3.9), entails that an intimate relation exists be- 
tween the kernel in the gap equation and that in the BSE. (This is another example 
of the coupling between DSEs.) Therefore an understanding of chiral symmetry and 
its dynamical breaking can only be obtained with a truncation scheme that preserves 
this relation, and hence guarantees Eq. (3.9) without a fine-tuning of model-dependent 
parameters. Until 1995-1996 no one had a good idea about how to do this. Equations 
were truncated, sometimes with good phenomenological results and sometimes with poor 
results. Neither the successes nor the failures could be explained. 

That changed with Refs. [96,97], which described a procedure, illustrated in Fig. 3.1, 
that generates a Bethe-Salpeter equation from the kernel of any gap equation whose 
diagrammatic content is known. That this is possible and achievable systematically is 
necessary and sufficient to prove numerous exact results in QCD. The procedure also en- 
ables the formulation of practical phenomenological models that can be used to illustrate 
the exact results and provide predictions for experiment with readily quantifiable errors. 



4 Exact Results in Hadron Physics 



Before detailing some of the exact results, it is worth explicating the range of applica- 
tions; e.g., veracious statements about the pion cr-term [102]; radially-excited and hybrid 
pseudoscalar mesons [107,108]; heavy- light [109] and heavy- heavy mesons [110]; novel 
results for the pion susceptibility obtained via analysis of the isovector-pseudoscalar vac- 
uum polarisation [111], which bear upon the essential content of the so-called "Mexican 
hat" potential that is used in building models for QCD; and a derivation [112] of the 
Weinberg sum rule [113]. 
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4.1 Pseudoscalar meson mass formula 

The first of tlie results was introduced in Ref. [114]; namely, a mass formula that is exact 
for flavour non-diagonal pseudoscalar mesons: 

where: mi are the current-masses of the quarks constituting the meson; and 

fH^^P^.^{0\qf2l5J^.qn\Ho-) = ^2 trcD / iTo'y^.Sf,{q+)^„^_ {q; P)Sf,{q-) , (4.2) 

J q 

ip„^^_ = -(0|q/,i75g/i 1^0- > = Zi trcD / 755/, (g+)rH^_ (g; P)Sf, (q-) , (4.3) 

Jq 

where Th _ is the pseudoscalar meson's bound-state Bethe-Salpeter amplitude: 

rH„_(fc;P) = T.[^EH^,-{k;P) + i-PFH„_{k-P) 

+7 • k Gh^_ (k; P) - a^,^k^,P,HH^_ {k; P) ] , (4.4) 

which is determined from the associated homogeneous BSE; viz., Eq. (3.10) without the 
inhomogeneity. 

Equation (4.2) is the pseudovector projection of the pseudoscalar meson's Bethe- 
Salpeter wave-function onto the origin in configuration space; namely, the meson's lep- 
tonic decay constant, which is the strong interaction contribution to the strength of 
the meson's weak interaction. Similarly, Eq. (4.3) is the pseudoscalar projection of the 
wave- function onto the origin in configuration space; viz., a pseudoscalar analogue of 
the meson's leptonic decay constant. These observations will subsequently become quite 
important herein. 

It is worth emphasising that the quark wave- function and Lagrangian mass renormal- 
isation constants, Z2,4(C, A), respectively, depend on the gauge parameter in precisely the 
manner needed to ensure that the right-hand sides of Eqs. (4.2), (4.3) are gauge- invariant. 
Moreover, Z2{(,A) ensures that the right-hand side of Eq. (4.2) is independent of both 
(" and A, so that fn _ is truly an observable; and Z4{(^,A) ensures that p)^ is inde- 
pendent of A and evolves with ( in just the way necessary to guarantee that the product 
rr6 p%[ is renormalisation-point-independent. 

In addition, it should be noted that Eq. (4.1) is valid for every pseudoscalar meson 
and for any value of the current-quark masses; viz., rhj^ G [0, cxj), i = 1, 2. This includes 
the chiral limit, in whose neighbourhood Eq. (4.1) can be shown [114] to reproduce the 
familiar Gell-Mann-Oakes-Renner relation [54], Eq. (1.13). It also includes arbitrarily 
large values of one or both current-quark masses. In the former case; viz., one quark mass 
fixed and the other becoming very large, so that rfiqjrfiQ <^ 1, then /h _ oc 1/y/mH _ , 
pH _ oc \/mH _, and one arrives at niH-^ = rng [109]. The mass formula thereby 
provides the path to the QCD proof of a potential model result. 

The axial- vector Ward-Takahashi identity, Eq. (3.9), is a crucial bridge to Eqs. (4.1) - 
(4.3); and on the way one can also prove the following Goldberger-Treiman-like relations 
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[114]: 

f\^_E„^^_{k-Q) = B\e), (4.5) 

Fn{k-0) + 2f%^_F„^^_{k-0) = A°(k^), (4.6) 

GR(k;0) + 2f^^_^_G„^^_{k;0) = rf^^''(^')' (4-7) 

fffl(fc;0) + 2/]^^_i/H^_(fc;0) = 0, (4.8) 

wherein the superscript indicates that the associated quantity is evaluated in the chiral 
limit, and Fr, Gr, Hr are analogues in the inhomogeneous axial- vector vertex of the 
scalar functions in the Hg- -meson's Bethe-Salpeter amplitude. 

These identities are of critical importance in QCD. The first, Eq. (4.5), can be used 
to prove that a massless pseudoscalar meson appears in the chiral-limit spectrum if, and 
only if, chiral symmetry is dynamically broken. Moreover, it exposes the fascinating con- 
sequence that the solution of the two-body pseudoscalar bound-state problem is almost 
completely known once the one-body problem is solved for the dressed-quark propaga- 
tor, with the relative momentum within the bound-state identified unambiguously with 
the momentum of the dressed-quark. This latter emphasises that Goldstone's theorem 
has a pointwise expression in QCD. The remaining three identities are also important 
because they show that a pseudoscalar meson must contain components of pseudovec- 
tor origin. This result overturned a misapprehension of twenty-years standing; namely, 
that only Eh _ (fc; 0) is nonzero [115]. These pseudovector components materially in- 
fiuence the observable properties of pseudoscalar mesons [43,44,77,116,117], something I 
will illustrate in Sec. 7, as do their analogues in other mesons [118-120]. 

I have explained that pseudoscalar mesons hold a special place in QCD, Sec. 1.4. 
There are three such states, composed of it,d-quarks, in the hadron spectrum with masses 
below 2GeV [55]: 7r(140); 7r(1300); and 7r(1800). Of these, the pion [7r(140)] is naturally 
well known and much studied. The other two are observed, e.g., as resonances in the 
coherent production of three-pion final states via pion- nucleus collisions [121]. In the 
context of a quantum mechanical constituent-quark-model Hamiltonian, these mesons 
are often viewed as the first three members of a QQ n ^Sq trajectory, where n is the 
principal quantum number; i.e., the 7r(140) is viewed as the 5-wave ground state and the 
others are its first two radial excitations. By this reasoning the properties of the 7r(1300) 
and 7r(1800) are likely to be sensitive to details of the long-range part of the quark-quark 
interaction because the constituent-quark wave-functions will possess material support at 
large interquark separation. There is a caveat here; namely, in comparison with 7r(1300), 
the 7r(1800) is narrow, with a width of 210MeV [55], and has a decay pattern that may be 
consistent with its interpretation as a hybrid meson in constituent-quark models [122].'' 
This picture has the constituent-quarks' spins aligned to produce Sqq — 1 with J — 
obtained by coupling Sqq to a spin-1 excitation of the confinement potential. Neverthe- 
less, from this perspective, too, the properties of 7r(1800) are sensitive to the interaction 
responsible for light-quark confinement. Hence the development of an understanding of 
the properties of these three pseudoscalars might provide information about light-quark 
confinement; and, as has been emphasised, developing an understanding of confinement 
is one of the most important goals in mathematics and physics. 



*It is worth providing some definitions. Exotic mesons - states whoso quantum numbers 
cannot be supported by quantum mechanical quark-antlquark systems; and hybrid mesons - 
states with quark-model quantum numbers but a non-quark-model decay pattern. Both systems 
are suspected to possess "constituent gluon" content, which translates into a statement that they 
are expected to have a large overlap with interpolating fields that explicitly contain gluon fields. 
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Figure 4.1. Lattice-QCD results for the ratio of the decay constants for the first-excited- 
and ground-state pseudoscalar mesons as a function of the pion mass squared. (Lattice 
parameters: volume = 16 x 32; fi — 5.2, spacing a ~ 0.1 fm, two flavours of degener- 
ate sea quarks; Wilson gauge action and clover fermions.) The "not improved" results 
were obtained from a fermion action with poor chiral symmetry properties. In this case 
|/,ri//iro| ~ 0.4, consistent with expectations based on quantum mechanics. The "im- 
proved" results were obtained through implementation of the full ALPHA method for 
the nonperturbative improvement of the fermion action, which greatly improves the sim- 
ulation's chiral symmetry properties. In this case, l/Tri/Jirol ~ 0.01. (N.B. The sign of 
the ratio was not determined in the simulation but is discussed in Ref. [95].) 



The puissance of Eq. (4.1) is further emphasised by the fact that it is also applicable 
here [107,108]. The formula holds at each pole common to the pseudoscalar and axial- 
vector vertices and therefore it also impacts upon the properties of non ground state 
pseudoscalar mesons. I will work with a label n > for the pseudoscalar mesons: 7r„, 
with n = denoting the ground state, n = 1 the state with the next lowest mass, and 
so on. Plainly, by assumption, m^^ ,(, > tti^q , and hence ro,^^ ,q > in the chiral limit. 
Moreover, the ultraviolet behaviour of the quark-antiquark scattering kernel in QCD 
guarantees that Eq. (4.3) is cutoff independent. Thus 



A ^ 

< Pvr, 



(C) := lim P^ 

m— ^0 



.(0 <oo, Vn, 



(4.9) 



where m is defined in Eq. (3.5). Hence, it is a necessary consequence of chiral symmetry 
and its dynamical breaking in QCD; viz., Eq. (4.1), that 



/." 



;0,Vn> 1. 



(4.10) 



This is the statement that Goldstone modes are the only pseudoscalar mesons to possess 
a nonzero leptonic decay constant in the chiral limit when chiral symmetry is dynamically 
broken. The decay constants of all other pseudoscalar mesons on this trajectory, e.g., 
radial excitations, vanish: DCSB impacts upon every pseudoscalar meson. On the other 
hand, in the absence of DCSB the leptonic decay constant of each such pseudoscalar 
meson vanishes in the chiral limit; i.e, Eq. (4.10) is true Vn > 0: absent DCSB, then all 
pseudoscalar mesons decouple from the weak interaction. 
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Figure 4.2. This simple flavour-mixing diagram is immaterial to the resolution of the 77- ry' 
conundrum, as is any collection of processes for which the figure may serve as a skeleton 
diagram. (Straight lines denote quarks and springs denote gluons.) 



From the perspective of quantum mechanics, Eq. (4.f0) is a surprising fact. The 
leptonic decay constant for S-wave states is typically proportional to the wave-function 
at the origin. Compared with the ground state, this is smaller for an excited state because 
the wave-function is broader in configuration space and wave-functions are normalised. 
However, it is a modest effect. For example, consider the e"'"e~ decay of vector mesons. A 
calculation in relativistic quantum mechanics based on light-front dynamics [123] yields 
l/pi / fpo 1 = 0-5, consistent with the value inferred from experiment and DSEs in rainbow- 
ladder truncation [94]: |/pi/,/po| = 0.45. Thus, it is not uncommon for Eq. (4.10) to be 
perceived as "remarkable" or "unbelievable." Notwithstanding this, in connection with 
the pion's first radial excitation, the value of /,ri ~ — 2MeV predicted in Ref. [107] is 
consistent with experiment [124] and simulations of lattice-QCD [125], as illustrated in 
Fig. 4.1. It is now recognised that the suppression of /,ri is a useful benchmark, which can 
be used to tune and validate lattice QCD techniques that try to determine the properties 
of excited states mesons. 



4.2 Flavourless pseudoscalar mesons 

In connection with electric-charge-neutral pseudoscalar mesons, Eq. (4.1) is strongly mod- 
ified owing to the non-Abelian anomaly. This entails that whilst the classical action 
associated with QCD is invariant under UA{Nf) (non-Abelian ax;ial transformations gen- 
erated by A075, where Ao oc diag[l, . . . , Ijv,]), the quantum field theory is not. The 
modification is particularly important to properties of rj and rj' mesons. The latter is 
obviously a peculiar pseudoscalar meson because its mass is far greater than that of any 
other light-quark pseudoscalar meson; e.g., m^i = 1.75 jn,,. N.B. The diagram depicted 
in Fig. 4.2 is often cited as central to a solution of the rj-rj' puzzle. However, as will be- 
come clear below, whilst it does contribute to flavour-mixing, the process is immaterial 
in resolving the rj-rj' conundrum. 

The correct mass formula for flavourless pseudoscalars follows from consideration of 
the complete UA{Nf) Ward-Takahashi identity: 

P^rt^ik; P) = S~\k+)i'y5T'' + i75.F"5-'(fc_) - 2iM'"'rl{k- P) - X(fc; P) , (4.11) 

which generalises Eq. (3.9). In Eq. (4.11), {J-'^\ a = 0, . . . , N^ — 1} are the generators of 
U{Nf) in the fundamental representation, orthonormalised according to txT'^T^ = \5°'^; 
the dressed-quark propagator S = diag[5'u, 5*^, Sa, Sc, St, ■ ■ ■] is matrix- valued; and 

M"" = trj. [{.;^^M'■"}J^''] , (4.12) 
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where TVC' is a matrix of renormalised current-quark masses and the trace is over flavour 
indices. 

The final term in the last line of Eq. (4.11) expresses the non-Abelian axial anomaly. 
It can be written 

X(fc;P) = S-\k+)5''"Au{k;P)S-'\k^), (4.13) 

Au{k;P) = fd^xd%e''''^+'^-''''y'>Nf{T°q{x)Q{0)q{y)), (4.14) 



and since A"^ "(fc; P) is a pseudoscalar, it has the general form 

A\k; P) = J^-ys [i£A{k; P)+'y PTA{k; P)+j-kk- PGAik; P) + a^.k^P,HA{k; P)] . 

(4.15) 
The matrix element in Eq. (4.14) represents an operator expectation value in full QCD; 
the operation in Eq. (4.13) amputates the external quark lines; and 

Q{x) = i^trc [e^.p.F^,,Fp^{x)] = d^,K^{x) (4.16) 

47r 

is the topological charge density operator, where the trace is over colour indices and 
F^y ~ ^^""Fl^y is the matrix- valued gluon fleld strength tensor. It is plain and important 
that only ^"^^ is nonzero. N.B. While 2(2;) is gauge invariant, the associated Chern- 
Simons current, A"^, is not. Thus in QCD no physical state can couple to K^ and hence 
no state which appears in the observable spectrum can contribute to a resolution of the 
so-called f/yi(l)-problem; i.e., physical states cannot play any role in ensuring that the 
ri' is not a Goldstone mode. 

As described in Sec. 4.1, if one imagines there are Nf massless quarks, then DCSB is a 
necessary and sufficient condition for the a 7^ components of Eq. (4.11) to guarantee the 
existence of N^ — 1 massless bound-states of a dressed-quark and -antiquark. However, 
owing to Eq. (4.13), a = in Eq. (4.11) requires special consideration. One case is easily 
covered; viz., it is clear that if A^ = 0, then the a = component of Eq. (4.11) is no 
different to the others and there is an additional massless bound-state in the chiral limit. 

On the other hand, the large disparity between the mass of the ly'-meson and the 
octet pseudoscalars suggests that A'^ 7^ in real- world QCD. If one carefully considers 
that possibility, then the Goldberger-Treiman relations in Eqs. (4.5) - (4.8) become [126] 

2f„^_EBs{k;0) = 2B'\k^)^£A{k;0), (4.17) 

FUk;0)+2fH^_FBs{k;0) = A^ik^) ~ TA{k;0), (4.18) 

G%{k-0) + 2,fH^^_GBsik;0) = 2^A\k^) - gA{k;0), (4.19) 

HUk;0) + 2fH^^_HBs{k;0) = -H^(fc;0), (4.20) 

It follows that the relationship 

£A{k;0)^2B'\k^), (4.21) 

is necessary and sufficient to guarantee that r5^(fc; P), the ffavourless pseudoscalar vertex, 
does not possess a massless pole in the chiral limit; i.e., that there are only Nf — 1 
massless Goldstone bosons. Now, in the chiral limit, B'^ik^) 7^ if, and only if, chiral 
symmetry is dynamically broken. Hence, the absence of an additional massless bound- 
state is only assured through the existence of an intimate connection between DCSB and 
an expectation value involving the topological charge density. 
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This critical connection is further highlighted by the following result, obtained through 
a few straightforward manipulations of Eqs. (4.11), (4.13) and (4.14): 

{gg)^J = -lim4 (m) = - lim Z4(C', A^) trcD / S°(g, C) (4.22) 

J q 

= ^ j d\{q{x)i^,q{x)Qm'-\ (4.23) 

The absence of a Goldstone boson in the a = channel is only guaranteed if this explicit 
identity between the chiral-limit in-meson condensate and a mixed vacuum polarisation 
involving the topological charge density is satisfied. (I will elucidate the concept of in- 
meson condensates [127-130] in Sec. 5.) 

Mass formulae valid for all pseudoscalar mesons have also been obtained [126] 



where 



fn^^ m\^ = 2 M"V^„_ + 5"° ««„_ , (4.24) 



f^^_P^ = Z2tv f T''t.1^xh^,-{<1;P), (4.25) 



1 

A 



iPH^SO = ^4tr A .F"75Xffo-('Z;^) 

J q 




(4.26) 



n^o- = V 9 ^«o-' ^«o- =(0121^0-)- (4.27) 



For charged pseudoscalar mesons, Eq. (4.24) is equivalent to Eq. (4.1), but the novelty of 
Eq. (4.24) is what it expresses for flavourless pseudoscalars. To illustrate, consider the 
case of a U{Nf = 3)-symmetric mass matrix, in which all Nf = 3 current-quark masses 
assume the single value m^, then this formula yields: 

m^,/°, =n„, -f2m^p°f. (4.28) 

Plainly, the rj' is split from the Goldstone modes so long as n^' 7^ 0. Numerical simula- 
tions of lattice-QCD have confirmed this identity [131,132]. 

It is important to elucidate the physical content of n^/ . Returning to the definition: 

^{0\Q\V'), (4.29) 

it is readily seen to be another type of in-meson condensate. It is analogous to those that 
will be discussed in Sec. 5 but in this case the hadron-to-vacuum transition amplitude 
measures the topological content of the rj' . One may therefore state that the rj' is split 
from the Goldstone modes so long as its wave-function possesses nonzero topological 
content. This is plainly very different to requiring that the QCD vacuum is topologically 
nontrivial. 

Within QCD the properties of the rj' can be computed via the BSE, just like other 
mesons. A nonzero value of n^i can be achieved with a Bethe-Salpeter kernel that 
contains the contribution depicted in Fig. 4.3 because one may argue from Eqs. (4.14) 
and (4.16) that an anomaly-related contribution to a meson's Bethe-Salpeter kernel can- 
not contain external quark or antiquark lines that are connected to the incoming lines: 
purely gluonic configurations must mediate, as illustrated in Fig. 4.3. Furthermore, it 
is straightforward to see that no finite sum of gluon exchanges can serve this purpose. 
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Figure 4.3. An illustration of the nature of the contribution to the Bethe-Salpeter kernel 
associated with the non-Abelian anomaly. All terms have the "hairpin" structure illus- 
trated in the lower panel. No finite sum of such intermediate states is sufficient. (Straight 
lines denote quarks, with /i and /2 independent, and springs denote gluons.) 



Indeed, consider any one such single contribution in the chiral limit. It will be propor- 
tional to the total momentum and hence vanish for P = 0, in conflict with Eq. (4.28). 
This lies behind the need for something like the Kogut-Susskind ghost; i.e., the coupling 
of a massless axial- vector gauge- like field to the Chern-Simons current, which does not 
appear in the particle spectrum of QCD because the current is not gauge invariant. (See 
Ref. [133] and Sec. 5.1 of Ref. [134].) 

It is argued [135,136] that in QCD with Nc colours. 



/Nc 



(4.30) 



and it can be seen to follow from the gap equation, the homogeneous BSE and Eqs. (4.25), 
(4.26) that 



f°,^VK^p'^:{0. 



One thus obtains 



= -^+2m(C)^ 

J ll' 



■°(C) 



•I n 



(4.31) 



(4.32) 



The first term vanishes in the limit Nc -^ oo while the second remains finite. Subsequently 
taking the chiral limit, the rj' mass approaches zero in the manner characteristic of 
all Goldstone modes. (N.B. One must take the limit Nc -^ oo before the chiral limit 
because the procedures do not commute [137].) These results are realised in the effective 
Lagrangian of Ref. [138] in a fashion that is consistent with all the constraints of the 
anomalous Ward identity. It is noteworthy that this is not true of the so-called 't Hooft 
determinant [134,139,140]. 

Implications of the mass formula in Eq. (4.24) were exemplified in Ref. [126] using 
an elementary dynamical model that includes a one-parameter Ansatz for that part of 
the Bethe-Salpeter kernel related to the non-Abelian anomaly, an illustration of which is 
provided in Fig. 4.3. The study compares ground-state pseudoscalar- and vector-mesons 
constituted from all known quarks, excluding the f-quark. Amongst the notable results 
is a prediction for the mixing angles between neutral mesons; e.g.. 



-15.4° , 9r,' = -15.7° 



(4.33) 
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N.B. There are necessarily two mixing angles, with each determined at the appropriate 
pole position in the inhomogeneous vertex. It is interesting that the angles are approxi- 
mately equal and compare well with the value inferred from a single mixing angle analysis 
[141] e = -13.3° ±1.0°. 

It is worth explicating the nature of the flavour-induced difference between the tv 
and TT masses. If one ignores mixing with mesons containing other than u, d-quarks; 
viz., works solely within SU{Nf = 2), then 771.^0 — rn^+ — — 0.04MeV. On the other 
hand, the full calculation yields m^o — m^+ — —0.4 MeV, a factor of ten greater, and one 
obtains a -K^-rj mixing angle, whose value at the neutral pion mass shell is 

e^^imlo) = 1.2°. (4.34) 

For comparison, Ref. [142] infers a mixing angle of 0.6° ± 0.3° from a i('-matrix analysis 
of the process pd — >■ ''HeTr". Plainly, mixing with the 77-meson is the dominant non- 
electromagnetic contribution to the iv -it" mass splitting. The analogous angle at the 77 
mass-shell is 

6i^^(m^) = 1.3°. (4.35) 

The angles in Eq. (4.33) correspond to 

Itt") ~ 0.72 uu- 0.69 dd- 0.013 ss, (4.36) 

\ri) ~ 0.53 MM -h 0.57 dd- 0.63 ss, (4.37) 

\ri') ~ O.U uu + 0A5 dd + 0.78 ss. (4.38) 

Evidently, in the presence of a sensible amount of isospin breaking, the tt" is still predom- 
inantly characterised by J^^ but there is a small admixture of ss. It is found in Ref. [126] 
that mixing with the tv° has a similarly modest impact on the flavour content of the r; 
and rj' . It's effect on their masses is far less. 



5 Condensates are Confined within Hadrons 

Dynamical chiral symmetry breaking and its connection with the generation of hadron 
masses was first considered in Ref. [75]. The effect was represented as a vacuum phe- 
nomenon. Two essentially inequivalent classes of ground-state were identified in the 
mean-field treatment of a meson-nucleon field theory: symmetry preserving (Wigner 
phase); and symmetry breaking (Nambu phase). Notably, within the symmetry breaking 
class, each of an uncountable infinity of distinct configurations is related to every other 
by a chiral rotation. This is arguably the origin of the concept that strongly-interacting 
quantum field theories possess a nontrivial vacuum. 

With the introduction of the parton model for the description of deep inelastic scat- 
tering (DIS), this notion was challenged via an argument [143] that DCSB can be realised 
as an intrinsic property of hadrons, instead of through a nontrivial vacuum exterior to 
the observable degrees of freedom. Such a perspective is tenable because the essential 
ingredient required for dynamical symmetry breaking in a composite system is the exis- 
tence of a divergent number of constituents and DIS provided evidence for the existence 
within every hadron of a sea of low-momentum partons. This view has, however, received 
scant attention. On the contrary, the introduction of QCD sum rules as a theoretical arti- 
fice to estimate nonperturbative strong-interaction matrix elements entrenched the belief 




Strong QCD and Dyson-Schwinger Equations 
Luminiferous Ether 



39 



Since the Earth is in motion, 
the flow of aether across 
the Earth should producea 
detectable "aether wind" 
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Figure 5.1. Physics theories of the late 19th century postulated that, just as water 
waves must have a medium to move across (water), and audible sound waves require a 
medium to move through (such as air or water), so also light waves require a medium, 
the "luminiferous aether." This was apparently unassailable logic until, of course, one of 
the most famous failed experiments in the history of science to date [144] . 



that the QCD vacuum is characterised by numerous distinct, spacetime-independent 
condensates. Faith in empirical vacuum condensates might be compared with an earlier 
misguided conviction that the universe was filled with a luminiferous aether. Fig. 5.1. 

Notwithstanding the prevalence of the belief in empirical vacuum condensates, it 
does lead to problems; e.g., entailing, as explained below, a cosmological constant that is 
10*®-times greater than that which is observed [128,145]. This unwelcome consequence 
is partly responsible for reconsideration of the possibility that the so-called vacuum con- 
densates are in fact an intrinsic property of hadrons. Namely, in a confining theory - and 
confinement is essential to this view - condensates are not constant, physical mass-scales 
that fill all spacetime; instead, they are merely mass-dimensioned parameters that serve 
a practical purpose in some theoretical truncation schemes but otherwise do not have an 
existence independent of hadrons [128-130,146-148]. 

To account for the gross overestimate, recall that in Sec. 1 I remarked that the ex- 
pansion of the universe is accelerating [145]. One piece of evidence is provided by the 
observations of type la supernovae reported in Refs. [32,33]; and another by measure- 
ments of the composition of the universe, which point to a missing energy component 
with negative pressure. To explain the latter, the curvature of the universe may be 
characterised by a mass-energy density, pu. There is a critical value of this density for 
which the universe is flat: pF- Observations of the cosmic microwave background's power 
spectrum, e.g., Ref. [149], indicate that fio = pu / pF = 1 ± 0.04. In a fiat universe, the 
matter density, Q.m, and the energy density must sum to the critical density. However, 
matter only contributes about one-third of the critical density, Q.m = 0.33 ± 0.04. Thus 
two-thirds of the critical density is missing. 

In order to have escaped detection, the missing energy must be smoothly distributed. 
In order not to interfere with the formation of structure (by inhibiting the growth of den- 
sity perturbations) the energy density in this component must change more slowly than 
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matter (so that it was subdominant in the past). The universe's accelerated expansion 
can be accommodated in general relativity through the cosmological constant, A, and 
observations determine an associated density, Pa'". (Recall that Einstein introduced the 
repulsive effect of the cosmological constant in order to balance the attractive gravity of 
matter so that a static universe was possible. However, he promptly discarded it after 
the discovery of the expansion of the Universe.) 

It has been suggested that the advent of quantum field theory makes consideration 
of the cosmological constant obligatory not optional [145]. Indeed, the only possible 
covariant form for the energy of the (quantum) vacuum; viz., 

T^i' ^ pvAcS^u (5.1) 

is mathematically equivalent to the cosmological constant. The vacuum is [145] "... a 
perfect fluid and precisely spatially uniform ..." so that "Vacuum energy is almost the 
perfect candidate for dark energy." Now, if the ground state of QCD is truly expressed 
in a nonzero spacetime- independent expectation value (qq), then the energy difference 
between the symmetric and broken phases is of order A/qcd ~ 0.3 GeV, as indicated by 
Fig. 2.9. One obtains therefrom: 

^QCD ^ ^o*V?'- (5.2) 

In fact, the discrepancy is far greater if the Higgs vacuum expectation value is treated in 
a similar manner. 

This mismatch has been called "the greatest embarrassment in theoretical physics." 
However, it vanishes if one discards the notion that condensates have a physical exis- 
tence, which is independent of the hadrons that express QCD's asymptotically realisable 
degrees of freedom [128]; namely, if one accepts that such condensates are merely mass- 
dimensioned parameters in one or another theoretical truncation scheme. This appears 
mandatory in a confining theory [129,130], a perspective one may embed in a broader con- 
text by considering just what is observable in quantum field theory [150]: ". . . although 
individual quantum field theories have of course a good deal of content, quantum field 
theory itself has no content beyond analyticity, unitarity, cluster decomposition and sym- 
metry." If QCD is a confining theory, then the principle of cluster decomposition is only 
realised for colour singlet states [79] and all observable consequences of the theory, in- 
cluding its ground state, can be expressed via an hadronic basis. This is quark-hadron 
duality. 

It is worthwhile to recapitulate upon the arguments in Refs. [129,130]. To begin, 
note that Eq. (4.2) is the exact expression in QCD for the leptonic decay constant of a 
pseudoscalar meson. It is a property of the pion and, as consideration of the integral 
expression reveals, it can be described as the pseudovector projection of the pion's Bethe- 
Salpeter wavefunction onto the origin in configuration space. Note that the product 
tp — STS is called the Bethe-Salpeter wavefunction because, when a nonrelativistic limit 
can validly be performed, the quantity ^ at fixed time becomes the quantum mechanical 
wavefunction for the system under consideration. (N.B. In the neighborhood of the 
chiral limit, a value for fn _ can be estimated via either of two approximation formulae 
[84,151,152]. These formulae both illustrate and emphasize the role of /h _ as an order 
parameter for DCSB.) 

If chiral symmetry were not dynamically broken, then in the neighborhood of the 
chiral limit /h _ oc rh [107] . Of course, chiral symmetry is dynamically broken in QCD 
[98-101] and so for the ground-state pseudoscalar 

hm /«^_ (m) = /]^ / . (5.3) 

771— )-0 ^ 
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Taken together, these last two observations express the fact that fn _ , which is an 
intrinsic property of the pseudoscalar meson, is a bona fide order parameter for DCSB. 
An analysis within chiral perturbation theory [153] suggests that the chiral limit value, 
/)^ _ , is ~ 5% below the measured value of 92.4 MeV; and efficacious DSE studies give 
a 3% chiral- limit reduction [127]. 

Now, Eq. (4.3) is kindred to Eq. (4.2); it is the expression in quantum field theory 
which describes the pseudoscalar projection of the pseudoscalar meson's Bethe-Salpeter 
wavefunction onto the origin in configuration space. It is thus truly just another type of 
pseudoscalar meson decay constant. 

In this connection it is therefore notable that one may rigorously define an "in-meson" 
condensate; viz. [114,127]: 

-{Qf2<lfi)H^_ = -fH„- (0|g/2 759/1 1^0-) = fH^.p-H^_ =: i4i^_ {rn) ■ (5.4) 

Now, using Eq. (4.5), one finds [114] 

lim 4 (m) = Z, trcD r^^.S^q; () = -{gqf, . (5.5) 

Hence the so-called vacuum quark condensate is, in fact, the chiral-limit value of the 
in-meson condensate; i.e., it describes a property of the chiral-limit pseudoscalar meson. 
One can therefore argue that this condensate is no more a property of the "vacuum" 
than the pseudoscalar meson's chiral-limit leptonic decay constant. Moreover, Ref. [154] 
establishes the equivalence of all three definitions of the so-called vacuum quark con- 
densate: a constant in the operator product expansion [155,156]; via the Banks-Casher 
formula [157]; and the trace of the chiral-limit dressed-quark propagator. Hence, they 
are all related to the in-meson condensate via Eq. (5.5) and none is defined essentially in 
connection with the vacuum. 

It is worth remarking that in the presence of confinement it is impossible to write 
a valid nonperturbative definition of a single quark or gluon annihilation operator; and 
therefore impossible to rigorously define a second quantised vacuum (ground state) for 
QCD upon a foundation of gluon and quark (quasiparticle) operators. To do so would 
be to answer the question: What is the state that is annihilated by an operator which 
is - as appears at present - unknowable? However, with the assumptions that confine- 
ment is absolute and that it entails quark-hadron duality, the question changes com- 
pletely. In this case, the nonperturbative Hamiltonian of observable phenomena in QCD 
is diagonalised by colour-singlet states alone. The ground state of this nonperturbative 
strong-interaction Hamiltonian is the state with zero hadrons. One may picture the cre- 
ation and annihilation operators for such states as rigorously defined via smeared sources 
on a spacetime lattice. The ground-state is defined with reference to such operators, 
employing, e.g., the Cell-Mann - Low theorem [158], which is applicable in this case be- 
cause there are well-defined asymptotic states and associated annihilation and creation 
operators. 

In learning that the so-called vacuum quark condensate is actually the chiral-limit 
value of an in-pion property, some respond as follows. The electromagnetic radius of 
any hadron which couples to pseudoscalar mesons must diverge in the chiral limit. This 
long-known effect arises because the propagation of massless on-shell colour-singlet pseu- 
doscalar mesons is undamped [159-162]. Therefore, does not each pion grow to fill 
the universe; so that, in this limit, the in-pion condensate reproduces the conventional 
paradigm? 

Confinement, again, vitiates this objection. Both DSE- and lattice-QCD studies 
indicate that confinement entails dynamical mass generation for both gluons and quarks. 
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see Sees. 2.3, 2.4. The dynamical gluon and quark masses remain large in the limit of 
vanishing current-quark mass. In fact, the dynamical masses are almost independent of 
the current-quark mass in the neighbourhood of the chiral limit. It follows that for any 
hadron the quark-gluon containment-radius does not diverge in the chiral limit. Instead, 
it is almost insensitive to the magnitude of the current-quark mass because the dynamical 
masses of the hadron's constituents are frozen at large values; viz., 2 — 3 Aqcd- These 
considerations show that the divergence of the electromagnetic radius does not correspond 
to expansion of a condensate from within the pion but rather to the copious production 
and subsequent propagation of composite pions, each of which contains a condensate 
whose value is essentially unchanged from its nonzero current-quark mass value within a 
containment-domain whose size is similarly unaffected. 

There is more to be said in connection with the definition and consequences of a chiral 
limit. Plainly, the existence of strongly-interacting massless composites would have an 
enormous impact on the evolution of the universe; and it is naive to imagine that one 
can simply set rhu.d = and consider a circumscribed range of manageable consequences 
whilst ignoring the wider implications for hadrons, the Standard Model and beyond. 
For example, with all else held constant, Big Bang Nucleosynthesis is very sensitive to 
the value of the pion-mass [67] . We are fortunate that the absence of quarks with zero 
current-quark mass has produced a universe in which we exist so that we may carefully 
ponder the alternative. 

The discussion of Ref. [129] was restricted to pseudoscalar mesons. It is expanded in 
Ref. [130] via a demonstration that the in-pseudoscalar-meson condensate can be repre- 
sented through the pseudoscalar-meson's scalar form factor at zero momentum transfer. 
With the aid of a mass formula for scalar mesons, revealed therein, the analogue was 
shown to be true for in-scalar-meson condensates. As argued, the concept is readily 
extended to all hadrons so that, via the zero momentum transfer value of any hadron's 
scalar form factor, one can readily extract the value for a quark condensate in that hadron 
which is a measure of dynamical chiral symmetry breaking. 

Given that quark condensates are an intrinsic property of hadrons, one arrives at 
a new paradigm, as observed in the popular science press [163]: "EMPTY space may 
really be empty. Though quantum theory suggests that a vacuum should be fizzing with 
particle activity, it turns out that this paradoxical picture of nothingness may not be 
needed. A calmer view of the vacuum would also help resolve a nagging inconsistency 
with dark energy, the elusive force thought to be speeding up the expansion of the 
universe." In connection with the cosmological constant, putting QCD condensates back 
into hadrons reduces the mismatch between experiment and theory by a factor of 10*^ . If 
technicolour-like theories are the correct scheme for extending the Standard Model [164], 
then the impact of the notion of in-hadron condensates is far greater still. 



6 Many Facets of DCSB 



The importance and interconnection of confinement and DCSB are summarised in Sec. 2.5; 
and some of the profound implications of DCSB for pseudoscalar mesons are detailed in 
Sec. 4. The latter could be proved owing to the existence of at least one systematic 
nonperturbative symmetry-preserving DSE truncation scheme (Fig. 3.1 and the associ- 
ated discussion). On the other hand, the practical application of this particular scheme 
has numerous shortcomings. For example, at leading-order (rainbow-ladder) the trun- 
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cation is accurate for ground-state vector- and electrically-charged pseudoscalar-mesons 
because corrections in these channels largely cancel, owing to parameter-free preservation 
of the Ward-Takahashi identities. However, they do not cancel in other channels [165- 
168]. Hence studies based on the rainbow-ladder truncation, or low-order improvements 
thereof, have usually provided poor results for scalar- and axial- vector-mesons [169-174], 
produced masses for exotic states that are too low in comparison with other estimates 
[94,95,169,172,174], and exhibit gross sensitivity to model parameters for tensor-mesons 
[175] and excited states [94,95,107,176]. In these circumstances one must conclude that 
physics important to these states is omitted. One anticipates therefore that signifi- 
cant qualitative advances in understanding the essence of QCD could be made with 
symmetry-preserving kernel Ansatze that express important additional nonperturbative 
effects, which are impossible to capture in any finite sum of contributions. Such an ap- 
proach has recently become available [152] and will be summarised below. It is hoped that 
the mathematics of Faa di Bruno Hopf algebras will assist hadron physics practitioners 
to improve further upon this. 



6.1 DCSB in the Bethe-Salpeter kernel 

In order to illustrate the decisive importance of DCSB in the Bethe-Salpeter kernel, 
consider, e.g., flavoured pseudoscalar and axial-vector mesons, which appear as poles in 
the inhomogeneous BSE for the axial- vector vertex, TA, where f,g are flavour labels. 
An exact form of that equation is (fc, q are relative momenta, P is the total momentum 
flowing into the vertex, and q± = qjz P/2, etc.) 



r^«(fc;P) = Z2JS1,.- f g'DMk^q)^7.Sfiq+)Ti'^iq;P)S,{q-)^rl{q.,k. 



A 



+ I g'D^,B{k-q)^Jc.Sf{q+)^Aif^^{k,q-P), (6.1) 



1 



where Ag „ is a 4-point Schwinger function. [The pseudoscalar vertex satisfies an ana- 
logue of Eq. (6.1).] This form of the BSE was first written in Ref. [167] and is illustrated 
in the lower-panel of Fig. 6.1. The diagrammatic content of the right-hand-side is com- 
pletely equivalent to that of Eq. (3.10), which is depicted in the upper-panel of the figure. 
However, in striking qualitative opposition to that textbook equation, Eq. (6.1) partly 
embeds the solution vertex in the four-point function. A, whilst simultaneously expli- 
cating a part of the effect of the dressed-quark-gluon vertex. This has the invaluable 
consequence of enabling the derivation of both an integral equation for the new Bethe- 
Salpeter kernel. A, in which the driving term is the dressed-quark-gluon vertex [167], and 
a Ward-Takahashi identity relating A to that vertex [152]. No similar equations have yet 
been found for K and hence the textbook form of the BSE, whilst tidy, is very limited 
in its capacity to expose the effects of DCSB in bound-state physics. 

As emphasised above, no study of light-quark hadrons is dependable if it fails to 
comply with the axial- vector Ward-Takahashi identity, Eq. (3.9). The condition 

P^A^'«^(fc,q;P)+i[m^(C)-fm,(C)]A^^(fc,q;P)=r^(Q+,fc+)i75 + n5r«(g_,fc_), (6.2) 

where Agg is the analogue of Ag^ o in the pseudoscalar equation, is necessary and sufficient 
to ensure the Ward-Takahashi identity is satisfied by the solution of Eqs. (3.1) and (6.1) 

[152]. 
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Figure 6.1. Omitting the inhoinogeneity, the upper panel illustrates the textbook form 
of the Bethe-Salpeter equation, Eq. (3.10), whereas the lower panel depicts the form ex- 
pressed in Eq. (6.1). The reversal of the total-momentum's flow is immaterial here. N.B. 
In any symmetry-preserving truncation, beyond the leading-order identified in Ref. [97]; 
i.e., rainbow-ladder, the Bethe-Salpeter kernel is nonplanar even if the vertex in the gap 
equation is planar [167]. This is illustrated in Fig. 3.1. 



Consider Eq. (6.2). Rainbow-ladder is the leading-order term in the systematic DSE 
truncation scheme of Refs. [96,97]. It corresponds to Vl =7^, in which case Eq. (6.2) 
is solved by Ag^g = = Ag|. This is the solution that indeed provides the rainbow- 



ladder forms of Eq. (6.1). Such consistency will be apparent in any valid systematic 
term-by-term improvement of the rainbow-ladder truncation. 

However, since the two-point functions of elementary excitations are strongly modified 
in the infrared, one must accept that the same is generally true for three-point functions; 
i.e., the vertices. Hence the bare vertex will be a poor approximation to the complete 
result unless there are extenuating circumstances. This is readily made apparent, for 
with a dressed-quark propagator of the form in Eq. (2.7), one finds immediately that the 
Ward-Takahashi identity is breached; viz., 



Pm^m / S~\k + P/2) - S-\k - P/2) , 



(6.3) 



and the violation is significant whenever and wherever the mass function in Fig. 2.9 
is large. This was actually realised early on, with studies of the fermion-gauge-boson 
vertex in Abelian gauge theories [177] that have inspired numerous ensuing analyses. The 
importance of this dressing to the reliable computation of hadron physics observables 
was exposed in Refs. [178,179], insights from which have subsequently been exploited 
effectively; e.g., Refs. [14,114,120,180-185]. 

The most important feature of the perturbative or bare vertex is that it cannot 
cause spin-flip transitions; namely, it is an helicity conserving interaction. However, one 
must expect that DCSB introduces nonperturbatively generated structures that very 
strongly break helicity conservation. These contributions will be large when the dressed- 
quark mass- function is large. Conversely, they will vanish in the ultraviolet; i.e., on the 
perturbative domain. The exact form of the vertex contributions is still the subject of 
study but their existence is model-independent. 

Critical now is a realisation that Eq. (6.2) is far more than just a device for checking 
a truncation's consistency. For, just as the vector Ward-Takahashi identity has long 
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been used to build Ansdtze for the dressed-quark-photon vertex [3,177,186,187], Eq. (6.2) 
provides a tool for constructing a symmetry-preserving kernel of the BSE that is matched 
to any reasonable form for the dressed-quark-gluon vertex which appears in the gap 
equation. With this powerful capacity, Eq. (6.2) achieves a goal that has been sought 
ever since the Bethe-Salpeter equation was introduced [106]. As will become apparent, it 
produces a symmetry-preserving kernel that promises to enable the first reliable Poincare 
invariant calculation of the spectrum of mesons with masses larger than 1 GeV. 

The utility of Eq. (6.2) was illustrated in Ref. [152] through an application to ground 
state pseudoscalar and scalar mesons composed of equal-mass u- and d-quarks. To this 
end, it was supposed that in Eq. (3.1) one employs an Ansatz for the quark-gluon vertex 
which satisfies 

P^iri(k+,k.)=B{P^) [Sj\k+)^Sj'{k^)] , (6.4) 

with B flavour-independent. (N.B. While the true quark-gluon vertex does not satisfy 
this identity, owing to the form of the Slavnov- Taylor identity which it does satisfy, it is 
plausible that a solution of Eq. (6.4) can provide a reasonable pointwise approximation 
to the true vertex [188].) Given Eq. (6.4), then Eq. (6.2) entails {I = q — k) 

ilp4${k, q- P) = B{f) \ci%q; P) - Vi%k- P)] , (6.5) 

with an analogous equation for P^lpiKl^' g{k, q\ P). This identity can be solved to obtain 

Ail{k,q-P)~B{{k-qf)^,,Kl'{k.,q-P), (6.6) 

with, using an obvious analogue of Eq. (4.4) , 

Tii'ik, q- P) = 2lp [iAe^ (g, fc; P) + 7 • PA^., (g, k; P)] + 7^ Eg^ (g, k- P) 

+ 2£p^-£AG,{<l,k;P) + [7,3,7 • P]^HAq,k;P) + 2^^ [7 • £, 7 • P]AH,(g,fc;P) , (6.7) 

where £ ^ {q + k)/2, E^{q, k; P) = [$(g; P) + $(fc; P)]/2 and As(g, k; P) = [$(g; P) - 
<E.(fc;P)]/[g2-fc2j^ 

Now, given any Ansatz for the quark-gluon vertex that satisfies Eq. (6.4), then the 
pseudoscalar analogue of Eq. (6.1), and Eqs. (3.1), (6.6), (6.7) provide a symmetry- 
preserving closed system whose solution predicts the properties of pseudoscalar mesons. 
The relevant scalar meson equations are readily derived. With these systems one can 
anticipate, elucidate and understand the influence on hadron properties of the rich non- 
perturbative structure expected of the fully-dressed quark-gluon vertex in QCD: in par- 
ticular, that of the dynamically generated dressed-quark mass function, whose impact is 
quashed at any finite order in the truncation scheme of Ref. [97] , or any kindred scheme. 

To proceed one need only specify the gap equation's kernel because, as noted above, 
the BSEs are completely defined therefrom. To complete the illustration [152] a simplified 
form of the effective interaction in Ref. [127] was employed and two vertex Ansdtze were 
compared; viz., the bare vertex E^ = 7^, which defines the rainbow-ladder truncation 
of the DSEs and omits vertex dressing; and the Ball-Chiu (BC) vertex [177], which 
nonperturbatively incorporates some of the vertex dressing associated with DCSB: 

ir» (g, k) = iEAs (g', fc") 7m + 2^M [n • ^ Aa9 (g', A:') + Ass (g^ fc')]. (6.8) 

A particular novelty of the study is that one can calculate the current-quark-mass- 
dependence of meson masses using a symmetry-preserving DSE truncation whose dia- 
grammatic content is unknown. That dependence is depicted in Fig. 6.2 and compared 
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Figure 6.2. Dependence of pseudoscalar (left panel) and scalar (right) meson masses 
on the current-quark mass, m. The Ball-Chiu vertex (BC) result is compared with the 
rainbow-ladder (RL) result. (Figure adapted from Ref. [152].) 



with the rainbow-ladder result. The m-dependence of the pseudoscalar meson's mass pro- 
vides numerical confirmation of the algebraic fact that the axial- vector Ward-Takahashi 
identity is preserved by both the rainbow-ladder truncation and the BC-consistent Ansatz 
for the Bethe-Salpeter kernel. The figure also shows that the axial- vector Ward-Takahashi 
identity and DCSB conspire to shield the pion's mass from material variation in response 
to dressing the quark-gluon vertex [9,167,168]. 

As noted in Ref. [152], a rainbow-ladder kernel with realistic interaction strength 
yields 



2M(0) 



2A/(0) 



(0.3 ±0.1), 



(6.9) 



which can be contrasted with the value obtained using the BC-consistent Bethe-Salpeter 
kernel; viz.. 



e^'' < 0.1 . 



(6.10) 



Plainly, significant additional repulsion is present in the BC-consistent truncation of the 
scalar BSE. 

Scalar mesons are commonly identified as ^Po states, see Fig. 6.3. This assignment 
expresses a constituent-quark-model perspective, from which a J^ = 0^^ fermion- 
antifermion bound-state must have the constituents' spins aligned and one unit of con- 
stituent orbital angular momentum. Hence a scalar is a spin and orbital excitation of 
a pseudoscalar meson. Of course, no constituent-quark-model can be connected sys- 
tematically with QCD. Nevertheless, the presence of orbital angular momentum in a 
hadron's rest frame is a necessary consequence of Poincare covariance and the momentum- 
dependent vector-boson-exchange character of QCD [10,110,189], so there is a realisation 
in QCD of the quark-model anticipation. 

Extant studies of realistic corrections to the rainbow-ladder truncation show that 
they reduce hyperfine splitting [168]. Hence, with the comparison between Eqs. (6.9) and 
(6.10) one has a clear indication that in a Poincare covariant treatment the BC-consistent 
truncation magnifies spin-orbit splitting, an effect which can be attributed to the influence 
of the quark's dynamically-enhanced scalar self-energy [10] in the Bethe-Salpeter kernel. 
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Pseudoscalar meson: L=0 & S=0 



Scalar meson: L=l & S=l, & L+S=0 



Figure 6.3. In constituent-quark-like models, pseudoscalar mesons are ^So states - con- 
stituent spins antiparallel and zero orbital angular momentum; and scalar mesons are 
^Po states - constituent spins parallel and one unit of orbital angular momentum. Hence 
a scalar is a spin and orbital excitation of a pseudoscalar meson. 



6.2 Quark anomalous magnetic moments 

It was conjectured in Ref. [152] that a full realisation of DCSB in the Bethe-Salpeter 
kernel will have a big impact on mesons with mass greater than 1 GeV. Moreover, that 
it can overcome a longstanding failure of theoretical hadron physics. Namely, no extant 
continuum hadron spectrum calculation is believable because all symmetry preserving 
studies produce a splitting between vector and axial- vector mesons that is far too small: 
just one-quarter of the experimental value (see, e.g., Refs. [170-173]). Significant devel- 
opments have followed that conjecture [182,184] and will now be related. 

As described in Sec. 2, when using Dirac's relativistic quantum mechanics, a fermion 
with charge q and mass m, interacting with an electromagnetic field, has a magnetic 
moment fi — q/[2m\ [59]. For the electron, this prediction held true for twenty years, 
until improvements in experimental techniques enabled the discovery of a small deviation 
[190], with the moment increased by a multiplicative factor: 1.00119 ± 0.00005. This 
correction was explained by the first systematic computation using renormalised quantum 
electrodynamics (QED) [191]: 



J_ 

2m 



1 + 



2tvJ 2m 



(6.11) 



where a is QED's fine structure constant. The agreement with experiment established 
quantum electrodynamics as a valid tool. The correction defines the electron's anomalous 
m,agnetic m,om,ent, which is now known with extraordinary precision and agrees with 
theory at 0(0^^) [192]. 

The fermion-photon coupling in QED is described by: 



d X iq ■ip{x)^^^{x) A^{x) 



(6.12) 



where ^^{x), ip(x) describe the fermion field and A^{x) describes the photon. As I ex- 
plained in Sec. 2, this interaction generates the following electromagnetic current for an 
on-shell Dirac fermion {k — pf — pi), 



iqu{pf) 



2m 



u{pi) , 



(6.13) 
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where: Fi{k^), F2{k^) are form factors; and u{p), u{p) are electron spinors. Using their 
Euchdean space definition, one can derive a Gordon-identity; viz., with 21 = pf -\- pi, 

2mu{pf)i-y^u{j)i) = u{j)f) [2£^ + ia^^k^,] u{pi) . (6.14) 

With this rearrangement one sees that for massive fermions the interaction can be de- 
composed into two terms: the first describes the spin-independent part of the fermion- 
photon interaction, and is common to spin-zero and spin-half particles, whilst the second 
expresses the spin-dependent, helicity flipping part. Moreover, one reads from Eqs. (6.13) 
and (6.14) that a point-particle in the absence of radiative corrections has Fi = 1 and 
F2 = 0, and hence Dirac's value for the magnetic moment. The anomalous magnetic 
moment in Eq. (6.11) corresponds to ^2(0) = a/2n. 

One infers from Eq. (6.14) that an anomalous contribution to the magnetic moment 
can be associated with an additional interaction term: 

d x^qip{x)a)j,ui>{x)Ff_iv{x) , (6.15) 

where Ffj,i,{x) is the gauge-boson field strength tensor. This term is invariant under local 
[/(I) gauge transformations but is not generated by minimal substitution in the action 
for a free Dirac field. 

Consider the effect of the global chiral transformation ^{x) — >■ exp{i9"/5)ip{x). The 
term in Eq. (6.12) is invariant. However, the interaction of Eq. (6.15) is not. These 
observations facilitate the understanding of a general result: F2 = for a massless 
fermion in a quantum field theory with chiral symmetry realized in the Wigner mode; i.e., 
when the symmetry is not dynamically broken. A firmer conclusion can be drawn. For 
m = it follows from Eq. (6.14) that Eq. (6.12) does not mix with the helicity-fiipping 
interaction of Eq. (6.15) and hence a massless fermion does not possess a measurable 
magnetic moment. 

A reconsideration of Ref. [191] reveals no manifest conflict with these facts. The 
perturbative expression for ^2(0) contains a multiplicative numerator factor of m and 
the usual analysis of the denominator involves steps that are only valid for m ^ 0. 
Fundamentally, there is no conundrum because QED is not an asymptotically free theory 
and hence, as mentioned in Sec. 1.1, does not have a well-defined nonperturbative chiral 
limit. 

On the other hand, in QCD the chiral limit is rigorously defined nonperturbatively 
[127]. (It remains to be seen whether the theory thus obtained is meaningful, as indi- 
cated in the antepenultimate paragraph of Sec. 5.) The analogue of Schwinger's one-loop 
calculation can then be carried out to find an anomalous c/iromo-magnetic moment for 
the quark. There are two diagrams in this case: one similar in form to that in QED; 
and another owing to the gluon self-interaction. One reads from Ref. [193] that the per- 
turbative result vanishes in the chiral limit. However, Fig. 2.9 demonstrates that chiral 
symmetry is dynamically broken in QCD and one must therefore ask whether this affects 
the chromomagnetic moment. 

Of course, it does; and it is now known that this is signalled by the appearance of 
Ass in Eq. (6.8). If one writes the quark-gluon vertex as 

iVfj.{pf,pi\k) = Xi{pf,pi;k)i'yfj, + 21^ 17 • £ X2{pf,Pi;k) + Xj{pj,p^;k) \ + .. . , (6.16) 

then contemporary simulations of lattice-regularised QCD [194] and DSE studies [188] 
agree that 

A3(p,P;0)«-^B(p^C) (6.17) 
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Figure 6.4. Direct anomalous chromomagnetic moment contribution to the dressed- 
quark-gluon vertex computed in quenched-QCD with current-quark mass m ~ 100 MeV 
[194]. The one-loop perturbative result is shown for comparison. Plainly, the nonpertur- 
bative result is two orders-of-magnitude larger than the perturbative computation. This 
level of magnification is typical of DCSB; e.g., with a current-quark mass of 4 MeV, one 
obtains M(p^ = 0) ~ 400 MeV. (Figure adapted from Ref. [194].) 



and also on the form of Ai(p,p;0), which is functionally similar to A{p ,(^). However, 
owing to non- orthogonality of the tensors accompanying Ai and A2, it is difficult to obtain 
a lattice signal for A2. One must therefore consider the DSE prediction for A2 in Ref. [188] 
more reliable. 

As pointed out above, perturbative massless-QCD conserves helicity so the quark- 
gluon vertex cannot perturbatively have a term with the helicity-flipping characteristics 
of A3. Equation (6.17) is thus remarkable, showing that the dressed-quark-gluon vertex 
contains at least one chirally-asymmetric component whose origin and size owe solely to 
DCSB; and Sec. 6.1 illustrates that A3 has a material impact on the hadron spectrum. 

This reasoning is extended in Ref. [182]: massless fermions in gauge field theories can- 
not possess an anomalous chromo/electro-magnetic moment because the term that de- 
scribes it couples left- and right-handed fermions; however, if chiral symmetry is strongly 
broken dynamically, then the fermions should also posses large anomalous magnetic mo- 
ments. Such an effect is expressed in the dressed-quark-gluon vertex via a term 



"=(p/,Pi;A:) = G^^k^Ts{pf,Pi,k) . 



(6.18) 



That QCD generates a strongly momentum-dependent chromomagnetic form factor 
in the quark-gluon vertex, T5, with a large DCSB-component, is confirmed in Ref. [194]. 
Only a particular kinematic arrangement was readily accessible in that lattice simulation 
but this is enough to learn that, at the current-quark mass considered: tz, is roughly two 
orders-of-magnitude larger than the perturbative form (see Fig. 6.4); and 



Vp'>0: ir5(p,-p;2p)l>!A3(p,p;0)| 



(6.19) 



The magnitude of the lattice result is consistent with instanton-liquid model estimates 
[195,196]. 



50 Craig D. Roberts 

This large chromomagnetic moment is likely to have a broad impact on the properties 
of light-quark systems [196,197]. In particular, as will be illustrated in Sec. 6.3, it can 
explain the longstanding puzzle of the mass splitting between the a±- and p-mesons 
in the hadron spectrum [184]. Here a different novelty will be elucidated; viz., the 
manner in which the quark's chromomagnetic moment generates a quark anomalous 
e/eciromagnetic moment. This demonstration is only possible now that the method of 
Ref. [152] is available. It was accomplished [182] using the same simplification of the 
effective interaction in Ref. [127] that produced Figs. 6.2. 

In order to understand the vertex Ansatz used in Ref. [182], it is necessary to return 
to perturbation theory. As mentioned above Eq. (6.17), one can determine from Ref. [193] 
that at leading-order in the coupling, Os, the three-gluon vertex does not contribute to 
the QCD analogue of Eq. (6.11) and the Abelian-like diagram produces the finite and 
negative correction (— as/[127r]). The complete cancellation of ultraviolet divergences 
occurs only because of the dynamical generation of another term in the transverse part 
of the quark-gluon vertex; namely, 

rM""*(p/,pO = [elj ■ k + i'if,(Jyplukp\TA(p;,pi) , (6.20) 

with T^i/ — 5fj,i, — k^k^/k^ , aj :— T^i/O^. (N.B. The tensor denominated F* here is 
labelled T* in Refs. [186,187].) 

Cognisant of this, one may use a simple Ansatz to express the dynamical generation 
of an anomalous chromomagnetic moment via the dressed-quark gluon vertex; viz., 

f^(p/,pO = Ff (p;,pO + rr"(p/,pO, (6.21) 

V^;'^[pf,p,) = Fr^^(p^,pO+Fr"^(p/,pO, (6.22) 



with T5{pf,Pi) = (-7/4)As(pj,Pi) and 



TiipfyPi) =J^{z) 



Sj? a ^„2 „2 



AB{pj,pt)-AA{pj,pi] 



(6.23) 



Me 

The damping factor J"(z) = (l-exp(-2))/2, z = (pf+p^ -2M|)/A3r, Ayr = 1 GeV, sim- 
plifies numerical analysis but is otherwise irrelevant; and Me = {p|p > 0,p^ = M^{p^)} 
is the Euclidean constituent-quark mass. A confined quark does not possess a mass-shell 
(Sec. 2). Hence, one cannot unambiguously assign a single value to its anomalous mag- 
netic moment. One can nonetheless compute a magnetic moment distribution. At each 
value of p^, spinors can be defined to satisfy the free-particle Euclidean Dirac equation 
with mass m — >■ M{p^) =: ?, so that 

«(p/;?)Fp(p/,p.;fc)u(p,;<;)=^t(p/)[Fi(fc')7^ + ^a^,fc.F2(fc')]«(pO (6.24) 



and then, from Eqs. (6.21) - (6.23), 



-2^VS[ 



^^""(0 = ^. .:' f .... . (6.25) 



2<;^5'. + 2^S' 



<; ' 



where ct^ — Ea (?,?), S\ = AA(<;,'i), etc. The numerator's simplicity owes to a pre- 
meditated cancellation between T4 and rs, which replicates the one at leading-order in 
perturbation theory. Where a comparison of terms is possible, this vertex Ansatz is 
semi-quantitatively in agreement with Refs. [194,188]. However, the presence and under- 
standing of the role of Y^f^"- is a novel contribution by Ref. [182]. N.B. It is apparent 
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Figure 6.5. Left panel - f-y (GeV"^) in Eq. (6.29) cf. (-7/4)As(p^p2), both computed 
using Eq. (6.21) and the same simphfication of the interaction in Ref. [127]. Right panel 
- Anomalous chromo- and electro-magnetic moment distributions for a dressed-quark, 
computed using Eq. (6.30). The dashed-curve in both panels is the rainbow-ladder (RL) 
truncation result. 



from Eq. (6.25) that k^™ oc m^ in the absence of DCSB, so that K^'^^/pm] ^ in the 
chiral limit. 

The BSE for the quark-photon vertex can be written following the method of Ref. [152] . 
Since the method guarantees preservation of the Ward-Takahashi identities, the general 
form of the solution is 



^l{Pf,Pi 



rM^(p/,pO + r^(p/,pO. 



(6.26) 



^tiiPf^Pi) = 'yfj, Fi +a^,^kvF2 +T^p(Tp,J„e- kFs + [e^'y ■ k + i'yfj_a„pe„kp]F4 



-iCh +C'~i-ke.-kh-tli-tFT + ela^pg^kpFs ., 



(6.27) 



where {Fi\i = 1, . . . , 8} are scalar functions of Lorentz-invariants constructed from pf, 
Pi, k. The Ward-Takahashi identity is plainly satisfied; viz., 



kpiTf,{pf,Pi) ^ kpirf,'^{pf,Pi) = S ^(p/)-S ^(p,). 



(6.28) 



Figure 6.5 depicts the results obtained for the quark's anomalous electromagnetic 
moment form factor 



/t,(p) := lim — — 



trcr^,fc^r2:(p/,p) = F2 + -p^Fs . 



(6.29) 



The result is evidently sizable. It is worth reiterating that f-y is completely nonperturba- 
tive: in the chiral limit, at any finite order in perturbation theory, f~, = 0. For contrast 
the figure also displays the result computed in the rainbow-ladder truncation of QCD's 
DSEs. As the leading-order in a systematic but stepwise symmetry-preserving scheme 
[97], this truncation only partially expresses DCSB: it is exhibited by the dressed-quark 
propagator but not present in the quark-gluon vertex. In this case f-y is nonzero but 
small. These are artefacts of the truncation that cannot be remedied at any finite order 
of the procedure in Ref. [97] or a kindred scheme. 
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Employing Eq. (6.24), in connection with the dressed-quark-photon vertex, one can 
write an expression for the quark's anomalous electromagnetic moment distribution 

2,F. + 2,^A + A40^ (6.30) 

where: A„(?) = 2<;^(5^ — 2';S''g — i^Fs — ';^Ft, and the Fi are evaluated at p} = pi — 
M{p'j)^ =: <j^, k^ = 0. Plainly, k(<;) = in the chiral limit when chiral symmetry is 
not dynamically broken. Moreover, as a consequence of asymptotic freedom, n{q) — >■ 
rapidly with increasing momentum. The computed distribution is depicted in Fig. 6.5. 
It yields Euclidean mass-shell values: 

Mf^„ = 0.44GeV, <T = -0.22 , <T = 0.45 

cf . M|l = 0.35 GeV, k^^ = , k^^ = 0.048. 

It is thus apparent that DCSB produces a dressed light-quark with a momentum- 
dependent anomalous chromomagnetic moment, which is large at infrared momenta. 
Significant amongst the consequences is the generation of an anomalous electromagnetic 
moment for the dressed light-quark with commensurate size but opposite sign. (N.B. 
This result was anticipated in Ref. [198], which argued that DCSB usually triggers the 
generation of a measurable anomalous magnetic moment for light-quarks.) The infrared 
dimension of both moments is determined by the Euclidean constituent-quark mass. This 
is two orders-of-magnitude greater than the physical light-quark current-mass, which sets 
the scale of the perturbative result for both these quantities. 

There are two more notable features; namely, the rainbow-ladder truncation, and 
low-order stepwise improvements thereof, underestimate these effects by an order of mag- 
nitude; and both the T4, and rs terms in the dressed-quark-gluon vertex are indispensable 
for a realistic description of hadron phenomena. Whilst a simple interaction was used to 
illustrate these outcomes, they are robust. 

These results are stimulating a reanalysis of hadron elastic and transition electromag- 
netic form factors [185,199], and the hadron spectrum, results of which will be described 
below. Furthermore, given the magnitude of the muon "(j^ — 2 anomaly" and its assumed 
importance as an harbinger of physics beyond the Standard Model [200], it might also 
be worthwhile to make a quantitative estimate of the contribution to (?^ — 2 from the 
quark's DCSB- induced anomalous moments following, e.g., the computational pattern 
for the hadronic light-by-light scattering component of the photon polarization tensor 
indicated in Ref. [201]. 



6.3 a^-p mass splitting 

The analysis in Ref. [152] enables one to construct a symmetry-preserving kernel for the 
BSE given any form for F^. Owing to the importance of symmetries in forming the 
spectrum of a quantum field theory, this is a pivotal advance. One may now use all 
information available, from any reliable source, to construct the best possible vertex 
Ansatz. The last section illustrated that this enables one to incorporate crucial nonper- 
turbative effects, which any finite sum of contributions is incapable of capturing, and 
thereby prove that DCSB generates material, momentum-dependent anomalous chromo- 
and electro-magnetic moments for dressed light-quarks. 

The vertex described in Sec. 6.2 contains a great deal of information about DCSB. It 
is the best motivated Ansatz to date, has stimulated a detailed reanalysis of the quark- 
photon coupling [187], and may be used in the calculation of the masses of ground-state 
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Figure 6.6. Illustration of the procedure used to determine meson masses, which is fully 
described in Ref. [202] and analogous to the method used in lattice-QCD. Solid curve - 
ai-meson, nonperturbative kernel; dot- dash- dash - ai, kernel derived from Eq. (6.8) only 
(Ball-Chiu, BC); and dash - a\, kernel derived from just the first term in Eq. (6.8) (IBC, 
a minimal renormalization improvement [203] of the leading-order - RL, rainbow-ladder 
- kernel [97]). Dot-dash curve - p-meson, nonperturbative kernel; Dot-dash-dot - p, 
BC-kernel; and dotted - p, IBC-kernel. Points - values of l/r(fc — 0; P^) in the given 
channel computed with the kernel described. Fade approximants are constructed in each 
case; and the location of the zero is identified with (— m^eson)- 



spin-zero and -one light-quark mesons in order to illuminate the impact of DCSB on the 
hadron spectrum. This analysis expands significantly on the discussion of scalar and 
pseudoscalar mesons in Sec. 6.1. 

A prediction for the spectrum therefore follows once the gap equation's kernel is 
specified and the Ward-Identity solved for Ag^„. In the pseudoscalar and axial- vector 
channels the Ward-Takahashi identity for the Bethe-Salpeter kernel is solved by 



2A5,8(M) = [r;3(?+,fc+)+75r,3(g-,fc-)75] 



rr5(M)(fc;^) 



s-i(-g+) + 5-1(9-) 



s-i(fc+) + s-i(-fc-) 

[75f^(g+,fc+)75-hf;3(g_,fc_)], (6.32) 



where F is the chosen Ansatz for the quark-gluon vertex. Kernels for other channels are 
readily constructed. 

Reference [184] computes ground-state masses using the interaction described in 
Ref. [94], which produced Fig. 2.8; the vertex model explicated in Sec. 6.2; and the method 
for solving the inhomogeneous Bethe-Salpeter equation that is detailed in Ref. [202], 
which ensures one need only solve the gap and Bethe-Salpeter equations at spacelike 
momenta. This simplifies the numerical problem. To explain, the inhomogeneous BSE 
is solved for the complete Bethe-Salpeter amplitude in a particular channel on a domain 
of spacelike total-momenta, P^ > 0. Any bound-state in that channel appears as a pole 



in the solution at P — 



Denoting the leading Chebyshev moment of the am- 



plitude's dominant Dirac structure by F(fc;P), where k is the relative momentum, then 
l/r(fc = 0; P^) exhibits a zero at (— m^eson)- The location of that zero is determined via 
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Table 6.1. Col. 1: Spectrum obtained with the full nonperturbative Bethe-Salpeter 
kernels described herein, which express effects of DCSB. The method of Ref. [202] was 
used: the error reveals the sensitivity to varying the order of Fade approximant. Col. 2 
- Experimental values; computed, except rria, from isospin mass-squared averages [55]. 
Col. 3 - Masses determined from the inhomogeneous BSE at leading-order in the DSE 
truncation scheme of Ref. [97] using the interaction in Ref. [152] (with this simple kernel, 
the Fade error is negligible); and Col. 4 - results in Ref. [204], obtained directly from the 
homogeneous BSE at the same order of truncation. 





Ref. [184] 


Expt. 


RL-Pade 


RL- direct 


TOtt 


0.138 


0.138 


0.138 


0.137 


nip 


0.817 ±0.016 


0.777 


0.754 


0.758 


ma- 


0.90 ±0.05 


0.4- 1.2 


0.645 


0.645 


rriai 


1.30 ±0.11 


1.24 ±0.04 


0.938 


0.927 


mbi 


1.15 ±0.07 


1.21 ±0.02 


0.904 


0.912 


rua, - mp 


0.48 ±0.12 


0.46 ±0.04 


0.18 


0.17 


rUb, - mp 


0.33 ±0.09 


0.43 ±0.02 


0.15 


0.15 



extrapolation of a Fade approximant to the spacelike-behavior of l/r(fc = 0; P^). This 
is illustrated for the p- and ai-channels in Fig. 6.6. 

A full set of results is listed in Table 6.1, wherein the level of agreement between 
Cols. 3 and 4 illustrates the efficacy of the method used to compute masses: no differ- 
ence is greater than 1%. Next consider mtr and compare Cols. 1-3. It is an algebraic 
result that in the RL-truncation of QCD's DSEs, rricr ~ 2M, where M is a constituent- 
like quark mass [205]. On the other hand, incorporating the quark mass function into 
the Bethe-Salpeter kernel via V^'^ generates a strong spin-orbit interaction, which signifi- 
cantly boosts JTio- [152]. This feature is evidently unaffected by the inclusion of VJf^^; i.e., 
those terms associated with a dressed-quark anomalous chromomagnetic moment. Since 
terms associated with pion final-state interactions were deliberately omitted from the 
nonperturbative kernel derived in Ref. [184], it is noteworthy that rricr in Col. 1 matches 
estimates for the mass of the dressed-quark-core component of the cr-meson obtained us- 
ing unitarised chiral perturbation theory [206,207]. N.B. In addition to providing a width, 
such final-state interactions necessarily reduce the real part of the mass [103,208,209]. 

Now compare the entries in Rows 2, 4-6. The p- and ai-mesons have been known for 
more than thirty years and are typically judged to be parity-partners; i.e., they would 
be degenerate if chiral symmetry were manifest in QCD. Flainly, they are not, being 
split by roughly 450 MeV (i.e., > mp/2). It is suspected that this large splitting owes to 
DCSB. Hitherto, however, no symmetry-preserving bound-state treatment could explain 
it. This is illustrated by Cols. 3, 4, which show that whilst a good estimate of rUp is 
readily obtained at leading-order in the systematic DSE truncation scheme of Ref. [97], 
the axial-vector masses are much underestimated. The fiaw persists at next-to-leading- 
order [170,173]. 

The analysis in Ref. [184] points to a remedy for this longstanding failure. Using the 
Foincare-covariant, symmetry preserving formulation of the meson bound-state prob- 
lem enabled by Ref. [152], with nonperturbative kernels for the gap and Bethe-Salpeter 
equations, which incorporate effects of DCSB that are impossible to capture in any 
step-by-step procedure for improving upon the rainbow-ladder truncation, it provides 
realistic estimates of axial-vector meson masses. In obtaining these results, Ref. [184] 
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showed that the vertex Ansatz used most widely in studies of DCSB, F^^' , is inade- 
quate as a tool in hadron physics. Used alone, it increases both lUp and mai but yields 
rriai — rUp = 0.21 GeV, qualitatively unchanged from the rainbow-ladder-like result (see 
Fig. 6.6). A good description of axial- vector mesons is only achieved by including interac- 
tions derived from T^'^; i.e., connected with the dressed-quark anomalous chromomag- 
netic moment [182]. Moreover, used alone, neither term in Fj^,'^™, Eq. (6.22), can produce 
a satisfactory result. The full vertex Ansatz and the associated gap and Bethe-Salpeter 
kernels described in Sec. 6.2 are the minimum required. 

Row 5 contains additional information. The leading-covariant in the 6i-meson channel 
is 75fc^. The appearance of k^ suggests that dressed-quark orbital angular momentum 
will play a significant role in this meson's structure, even more so than in the ai-channel 
for which the dominant covariant is 757^,. (N.B. In a simple quark-model, constituent 
spins are parallel within the ai but antiparallel within the 61 . Constituents of the 61 may 
therefore become closer, so that spin-orbit repulsion can exert a greater influence.) This 
expectation is borne out by the following: with the full kernel, rrib^ is far more sensitive 
to the interaction's momentum-space range parameter than any other state, decreasing 
rapidly as the interaction's spatial- variation is increasingly suppressed. 

The results reviewed in this subsection rest on an Ansatz for the quark-gluon vertex 
and whilst the best available information was used in its construction, improvement 
is nonetheless possible. That will involve elucidating the role of Dirac covariants in 
the quark-gluon vertex which have not yet been considered, as in Ref. [187], and of 
resonant contributions; viz., meson loop effects that give widths to some of the states 
considered. In cases for which empirical width-to-mass ratios are < 25%, one might judge 
that such contributions can reliably be obtained via bound-state perturbation theory 
[210]. Contemporary studies indicate that these effects reduce bound-state masses but 
the reduction can uniformly be compensated by a modest inflation of the interaction's 
mass-scale [205,211], so that the masses in Table 6.1 are semiquantitatively unchanged. 
The case of the cr-meson is more complicated. However, the prediction of a large mass for 
this meson's dressed-quark core leaves sufficient room for a strong reduction by resonant 
contributions [206,207]. 

This section reviewed a continuum framework for computing and explaining the me- 
son spectrum, which combines a veracious description of pion properties with estimates 
for masses of light-quark mesons heavier than nip. (A contemporary lattice-QCD per- 
spective on this problem may be drawn from Refs. [212,213].) The method therefore offers 
the promise of a first reliable Poincare-invariant, symmetry-preserving computation of 
the spectrum of light-quark hybrids and exotics; i.e., those putative states which are 
impossible to construct in a quantum mechanics based upon constituent-quark degrees- 
of-freedom. So long as the promise is promptly fulfilled, the approach will provide pre- 
dictions to guide the forthcoming generation of facilities. 



7 Probing the Hadronic Interior 
7.1 Elastic form factors 

Form factors have long been recognised as a basic tool for elucidating bound-state prop- 
erties and are of particular value in hadron physics because they provide information 
on hadron structure as a function of Q^, the squared momentum-transfer: small-Q^ is 
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Figure 7.1. Apart from the direct contribution (left image), electron scattering from a 
nucleon may also proceed via meson- nucleon final state interactions (right image) , which 
have a measurable effect. 



the nonperturbative domain; and large-Q^ is the perturbative domain. Nonperturbative 
methods in hadron physics must explain their behaviour from Q^ = into and through 
the transition domain, whereupon the behaviour is currently being measured. Experi- 
mental and theoretical studies of hadron electromagnetic form factors have made rapid 
and significant progress during the last several years, including new data in the timelike 
region, and material gains have been made in studying the pion form factor. 

Despite this, many urgent questions remain unanswered. For example; how can 
one use experiment to chart the long-range behaviour of the /3-function in QCD; given 
the low mass of the pion and its strong coupling to protons and neutrons, how can 
one discern features produced by final-state interactions and disentangle them from the 
intrinsic properties of the target, Fig. 7.1; and at which momentum-transfer does the 
transition from nonperturbative-QCD to perturbative-QCD take place? Overviews of 
the recent status of experiment and theory are provided in Refs. [214,215]. However, 
with the intense interest and investment in this field, the status changes rapidly. 

7.1.1 Pion Form Factor. As an illustration, I will consider the electromagnetic 
charged-pion form factor. Measurement of this form factor is not simple because the pion 
is an unstable meson with a short lifetime - roughly 30ns - and hence one cannot construct 
a stable pion target. On the other hand, since the pion is often described as a two body 
system, many theorists have pretended that a reliable computation of the form factor is 
easy. However, owing to the intimate connection between the pion and DCSB, that is 
not the case and instead one must employ a tool that is capable of veraciously expressing 
that relationship. The DSEs fulfill that requirement and a quantitative prediction was 
obtained [120] by combining the: dressed- rainbow gap equation, Eq. (3.1) with Vl — >■ 7,^; 
dressed-ladder Bethe-Salpeter equation, Eq. (3.10) with the kernel in the top row of 
Fig. 3.1; and dressed-impulse approximation to the form factor. For an incoming pion 
with momentum pi = K — Q/2 and an outgoing pion with momentum p2 — K + Q/2, 
the latter is provided by 



iK^FriQ') = ^-/(£jltrD 



r^(t;-P2) 



xS{t + P2)Tf,{t + P2,t + pi)S{t + Pi) r^(f ; Pi) s{t) 



(7.1) 



where the trace is over spinor indices alone. This expression is represented by the diagram 
in Fig. 7.2. 

The impulse approximation to the form factor is completely determined once the 
interaction in the gap equation is specified; viz., a choice made for Z\g^ D^j,v{k)T^{q,p) 
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' ^Jq^Q) - Dressed-quark-photon vertex: 
■ Computed via irhomogeneous 

\ Bethe-Salpeter equation 



■5fpj - dressed-quark propagator 
computed via the Gap Equation 




\. 



r^(ti;P) - Pion Bethe-Salpeter amplitude: 
computed via tine Bethe-Saipeter equation 

Figure 7.2. Generalised impulse approximation to the electromagnetic charged-pion form 
factor [179]. As depicted, the relevant Ward-Takahashi identities are satisfied so long as 
each element in the computation is obtained from the appropriate DSE at leading-order 
in the truncation scheme of Ref. [97]. 



in Eq. (3.1). That employed in Ref. [120] grew from the ground laid by a series of studies 
[80,118,127,216-219]. It combines a one-parameter model for the infrared behaviour of 
the interaction with renormalisation-group-improved one-gluon exchange at ultraviolet 
momenta so that the results of perturbative QCD are recovered. In connection with 
ground-state vector and flavoured pseudoscalar mesons, this interaction correlates more 
than forty observables with a root-mean-square relative error of 15% [5] and has hitherto 
provided the most efficacious tool for JLab physics. It is corrected and updated in 
Ref. [94]. The pion for factor obtained as a parameter- free prediction in Ref. [120] was 
confirmed by subsequent experiments at ,ILab, as illustrated by the dashed curve in 
Fig. 7.3. 

Another of the reasons I focus on F^"^{Q^) is the existence of a prediction that 
Q^FTr{Q^) w constant for Q^ ^ m^ in a theory whose interaction is mediated by massless 
vector-bosons. To be explicit [223-225]: 

Q^F^iQ^) "^ -'"" 16^/^0(0^), (7.2) 

which takes the value 0.13 GeV^ at Q^ = 10 GeV^ if one uses the one-loop result for the 
strong running-coupling a{Q^ = lOGeV^) ~ 0.3. The verification of this prediction is 
a strong motivation for modern experiment [220-222], which can also be viewed as an 
attempt to constrain and map experimentally the pointwise behaviour of the exchange 
interaction that binds the pion. 

Section 4.1 details some extraordinary consequences of DCSB, amongst them the 
Goldberger-Treiman relations of Eqs. (4.5) - (4.8). Of these, Eqs. (4.6) and (4.7) entail 
that the pion possesses components of pseudovector origin which alter the asymptotic 
form of -F™(Q'^) by a multiplicative factor of Q^ cf. the result obtained in their absence 
[116]. Here, in connection with Eq. (7.2), Fig. 2.9 and Eqs. (4.5) - (4.8) are critical. 

In the electromagnetic elastic scattering process, the momentum transfer, Q, is pri- 
marily shared equally between the pion's constituents because the bound-state Bethe- 
Salpeter amplitude is peaked at zero relative momentum. Thus, one can consider k ~ 
Q/2. The Goldberger-Trieman-like relations express a mapping between the relative 
momentum of the pion's constituents and the one-body momentum of dressed-quark; 
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Figure 7.3. Q^F°™(Q^) - Dashed curve : DSE prediction [120], which employed a 
momentum-dependent renormalisation-group-improved gluon exchange interaction and 
produces r-n = 0.67 fm cf. experiment [55] 0.672 ± 0.0081m. Owing to the interaction's 
momentum-dependence, this computation involves a dressed-quark mass-lunction which 
evolves in the chiral limit as ~ l/p'^ for p^ ^ ^qcd- Solid curve : Q^F^'^{Q^) ob- 
tained with a symmetry-preserving DSE treatment of a momentum-independent one- 
gluon exchange interaction, Sec. 2.2 and Refs. [77,117]. In this computation, the com- 
puted dressed-quark mass is momentum independent; viz., M(p^) = M. For Q^ > 
0.17 GeV^ « M^, marked by the vertical dotted line, the contact interaction result for 
F^'-'^iQ^) differs from that in Ref. [120] by more than 20%. The data are from Refs. [220- 
222]. 



and the momentum dependence of the dressed-quark mass function is well-described by 
perturbation theory when fc^ > 2 GeV^ . Hence, one should expect a perturbative-QCD 



analysis of the pion form factor to be valid for k = Q /4 > 2 GeV ; i.e. 



Fr (Q') ^ f: 



icm pQCD 



(Q^)forQ^ >8GeV^ 



(7.3) 



This explains the result in Ref. [116], which is reproduced in Fig. 7.4. A similar argument 
for baryons suggests that the nucleon form factors should be perturbative for Q^ > 
ISGeV^ 



7.2 Deep inelastic scattering 

As mentioned in Sec. 2, quarks were discovered in deep inelastic scattering (DIS) experi- 
ments at SLAG, performed during the period 1966-1978. DIS is completely different to 
elastic scattering. In this process, one disintegrates the target instead of keeping only 
those events in which it remains intact. On a well-defined kinematic domain; namely, 
the Bjorken limit: 



OO , XBjoriicii . — 



2P- 



= fixed, 



(7.4) 
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Figure 7.4. Influence of tfie pion's pseudovector components on tfie asymptotic befiaviour 
of Q^F^^{Q^). The existence of such components is a necessary consequence of DCSB 
[114,116]. Dashed curve - Result obtained if one ignores the pseudovector components, 
Ftt{Q^) ~ 1/Q*; and solid curves - result produced via inclusion of the pseudovector 
parts of the pion's Bethe-Salpeter amplitude, Q^F-^{Q^) ~ constant. In the latter case 
the two curves result from slightly different assumptions about the evolution of the strong 
running-coupling from the infrared to the ultraviolet, something which again exposes the 
sensitivity of this observable to the momentum-dependence of the strong-interaction's 
/3-function. 



where P is the target's four-momentum and q is the momentum transfer, the cross-section 
can rigorously be interpreted as a measurement of the momentum-fraction probability 
distribution for quarks and gluons within the target hadron: q(x), g{x). These quantities 
describe the probability that a quark/gluon within the target will carry a fraction x of 
the bound-state's momentum, as defined in the infinite-momentum or light-front frame. 
(The light-front formulation of quantum field theory is built upon Dirac's front form of 
relativistic quantum dynamics [226].) 

The past forty years have seen a tremendous effort to deduce the parton distribution 
functions (PDFs) of the most accessible hadrons - the proton, neutron and pion. There 
are many reasons for this long sustained and thriving interest [11] but in large part it is 
motivated by the suspected process-independence of the usual parton distribution func- 
tions and hence an ability to unify many hadronic processes through their computation. 
In connection with uncovering the essence of the strong interaction, the behaviour of the 
valence-quark distribution functions at large Bjorken-i is most relevant. Furthermore, 
an accurate determination of the behavior of distribution functions in the valence region 
is also important to high-energy physics. Particle discovery experiments and Standard 
Model tests with colliders are only possible if the QCD background is completely under- 
stood. QCD evolution, apparent in the so-called scaling violations by parton distribution 
functions,'' entails that with increasing center-of-mass energy, the support at large- a; in 
the distributions evolves to small-x and thereby contributes materially to the collider 



^DGLAP evolution is described in Sec.IID of Ref. [11], The evolution equations are derived 
in perturbative QCD and determine the rate of change of parton densities when the energy-scale 
chosen for their definition is varied. 



60 Craig D. Roberts 




Pion ^-^ ^vProton 



Figure 7.5. nN Drell-Yan process, in which, e.g., a valence antiquark from the pion 
annihilates with a valence quark in the nucleon to produce a fi^ fi~ pair. 



background. N.B. The nucleon PDFs are now fairly well determined for x < 0.8 but the 
pion and kaon PDFs remain poorly known on the entire domain of x. 

7.2.1 Pion and kaon valence-quark distributions. Owing to the dichotomous 
nature of Goldstone bosons, understanding the valence-quark distribution functions in 
the pion and kaon is of great importance. Moreover, given the large value of the ratio of 
s-to-u current-quark masses, a comparison between the pion and kaon structure functions 
offers the chance to chart effects of explicit chiral symmetry breaking on the structure 
of would-be Goldstone modes. There is also the prediction [227,228] that a theory in 
which the quarks interact via 1/fc^-vector-boson exchange will produce valence-quark 
distribution functions for which 

q^{x)(x{l-xf+^ , x>0.S5, (7.5) 

where 7 > is an anomalous dimension that grows with increasing momentum transfer. 
(See Sec. VLB. 3 of Ref. [11] for a detailed discussion.) 

Owing to the absence of pseudoscalar meson targets, experimental knowledge of the 
parton structure of the pion and kaon arises primarily from pionic or kaonic Drell-Yan 
processes, illustrated in Fig. 7.5, involving nucleons in heavy nuclei [229-233]. Theoreti- 
cally, given that DCSB plays a crucial role in connection with pseudoscalar mesons, one 
must employ an approach that realistically expresses this phenomenon. The DSEs there- 
fore provide a natural framework: studies of the pion and kaon exist and will be reviewed 
here. The first [41] computed pion PDFs, using efficacious parametrisations of both the 
Bethe-Salpeter amplitude and dressed-quark propagators [179,234,235]. The second [43] 
employed direct, numerical DSE solutions in the computation of the pion and kaon PDFs, 
adapting the approach employed in successful predictions of electromagnetic form factors 
[108,119,120,236]; and also studied the ratio UK{x)/uTr{x) in order to elucidate aspects 
of the influence of an hadronic environment. 

fn rainbow-ladder truncation, one obtains the pion's valence-quark distribution from 

"-W^-i f 7^^^-cd ir^{e,^p) S4£)r"{i-x) S4£)r^{£,p) SdH -^ p)] , (7.6) 

2 J [Itv)* 

wherein the Bethe-Salpeter amplitudes and dressed-quark propagators are discussed 
above and F" {£; x) is a generalization of the dressed-quark-photon vertex, describing 
a dressed-quark scattering from a zero momentum photon. It satisfies a BSE (here with 
a rainbow-ladder kernel) with the inhomogeneous term i'y ■ n 5{£ ■ n — xP ■ n) . In Eq. (7.6), 
n^ is a light-like vector satisfying n^ = 0. In choosing rainbow-ladder truncation one im- 
plements a precise parallel to the symmetry-preserving treatment of the pion charge form 
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Figure 7.6. Pion valence quark distribution function evolved to (5.2 GeV)^. Solid curve 

- full DSE calculation [43]; dot-dashed curve - semi-phenomenological DSE-based calcu- 
lation in Ref. [41]; filled circles - experimental data from Ref. [232], at scale (4.05 GeV)^; 
dashed curve - NLO reanalysis of the experimental data [233]; and dot- dot- dashed curve 

- NLO reanalysis of experimental data with inclusion of soft-gluon resummation [42] . 



factor at Q^ — 0, wherein the vector current is conserved by use of ladder dynamics at all 
three vertices and rainbow dynamics for all three quark propagators [116,119,120,179]. 
Equation (7.6) ensures automatically that 



(Xf) :— / dxqj{x) — 1 for f = u,d, 



(7.7) 



since J dxT"{£; x) gives the Ward-identity vertex and the Bethe-Salpeter amplitudes are 
canonically normalised. 

Figure 7.6 displays the DSE results for the pion's valence u-quark distribution, evolved 
from a resolving scale Qo = (0.6 GeV)^ to Q^ = (5.2 GeV)'^ using leading-order DGLAP 
evolution (Sec.IID of Ref. [11]), and a comparison with ttN Drell-Yan data [232] at a scale 
Q^ ~ (4.05 GeV)^, inferred via a leading-order analysis. The computation's resolving 
scale, Qo, was fixed by matching the {x")^ moments for n = 1, 2,3 to an experimental 
analysis at (2GeV)^ [237]. 

It is notable that at Qo the DSE results yield 



2(x)5„=0.7, 2(x)^„=0.8. 



(7.8) 



(For comparison, the parametrised valence-like pion parton distributions of Ref. [238] 
yield a gluon momentum fraction of {xg)Q g^^ = 0.3.) In each case the remainder of the 
hadron's momentum is carried by gluons, which effect binding of the meson bound state 
and are invisible to the electromagnetic probe. Some fraction of the hadron's momentum 
is carried by gluons at all resolving scales unless the hadron is a point particle [11]. 
Indeed, it is a simple algebraic exercise to demonstrate that the only non-increasing, 
convex function which can produce (x^) = 1 and {x) = i, is the distribution u{x) = 1, 
which is uniquely connected with a pointlike meson; viz., a meson whose Bethe-Salpeter 
amplitude is momentum-independent. Thus Eqs. (7.8) are an essential consequence of 
momentum conservation. 
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Whilst the DSE results in Fig. 7.6 are both consistent with Eq. (7.5); i.e., they produce 
algebraically the precise behaviour predicted by perturbative QCD, on the valence-quark 
domain - which is uniquely sensitive to the behaviour of the dressed-quark mass-function, 
M{p^) " it is evident that they disagree markedly with the Drell-Yan data reported in 
Ref. [232]. This tension was long seen as a crucial mystery for a QCD description of 
the lightest and subtlest hadron [11]. Its re-emergence with Ref [41] motivated a NLO 
reanalysis of the Drell-Yan data [233], the result of which is also displayed in Fig. 7.6. 
At NLO the extracted PDF is softer at high- a; but the discrepancy nevertheless remains. 
To be precise, Ref. [233] determined a high-a: exponent of /? ~ 1.5 whereas the exponents 
produced by the DSE studies [41,43] are, respectively, 2.1 and 2.4 at the common model 
scale. They do not allow much room for a harder PDF at high- a;. ^ 

Following the highlighting of this discrepancy in Ref. [11], a resolution of the conflict 
between data and well-constrained theory was proposed. In Ref. [42] a long-overlooked 
effect was incorporated; namely, "soft-gluon resummation." With the inclusion of this 
next-to-leading-logarithmic threshold resummation effect in the calculation of the Drell- 
Yan cross section, a considerably softer valence-quark distribution was obtained at high- 
X. This is readily understood. The Drell-Yan cross-section factorises into two pieces; 
one hard and the other soft. The soft piece involves the PDF and the hard piece is 
calculable in perturbation theory. Adding additional interactions to the latter, which are 
important at large-x; viz., soft gluons, provides greater strength in the hard piece on the 
valence-quark domain. Hence a description of the data is obtained with a softer PDF. 
Indeed, the distribution obtained thereby matches precisely the expectations based on 
perturbative-QCD and obtained using DSEs. This is evident in a comparison between 
the dash-dot and dash-dot-dot curves in Fig. 7.6. This outcome again emphasises the 
predictive power and strength of using a single internally-consistent, well-constrained 
framework to correlate and unify the description of hadron observables. 

The ratio uk /u-^ measures the effect of the local hadronic environment. In the kaon, 
the u-quark is bound with a heavier partner than in the pion (rris ~ 25md) and this 
should cause uk{x) to peak at lower-a; than u-n^x). In fact, one finds that the s-quark 
distribution peaks at a value of x which is just 15% larger than that of the u-quark. 
Hence, even though DIS is a high-Q'^ process, constituent-quark-like mass-scales explain 
this shift: Ms/Mu « 1.25. The DSE prediction of uk/u^ [11,43] is shown in Fig. 7.7 
along with available Drell-Yan data [229,230]. The parameter-free DSE result agrees well 
with the data. I note that 

UKJQ) DGLAPiQ^^oo ., ,„^, 

«7r(0) 

viz, the ratio approaches one under evolution to larger resolving scales owing to the 
increasingly large population of sea-quarks produced thereby [239] . On the other hand, 
the value at a:: = 1 is a fixed-point under evolution: 



VQ? > Ql ^ 



DGLAP:Q2-^Qf UKiX) 



r(l) 



uk{1) 



U,r(l) 



(7.10) 



Ql 



i.e., it is the same at every value of the resolving scale Q^, and is therefore a persistent 
probe of nonperturbative dynamics [11]. 

With Ref. [43] a significant milestone was achieved; viz., unification of the computa- 
tion of distribution functions that arise in analyses of deep inelastic scattering with that 



^Here, "hard" means pointlike but can also identify a process with large momentum transfer. 
On the other hand, "soft" is associated with the structure of diffuse composite objects or processes 
with small momentum transfer. 
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Figure 7.7. DSE prediction for the ratio of M-quark distributions in the kaon and pion 
[11,43]. The full Bethe-Salpeter amplitude produces the solid curve; a reduced BSE ver- 
tex produces the dashed curve. The reduced amplitude retains only the invariants and 
amplitudes involving pseudoscalar and axial vector Dirac matrices, and ignores depen- 
dence on the variable q ■ P. These are part of the reductions that occur in a pointlike 
treatment of pseudoscalar mesons. The experimental data is from [229,230]. 



of numerous other properties of pseudoscalar mesons, including meson-meson scattering 
[240,241] and the successful prediction of electromagnetic elastic and transition form fac- 
tors. The results confirm the large- a; behavior of distribution functions predicted by the 
QCD parton model; provide a good account of the tt-N Drell-Yan data for UTr{x); and a 
parameter- free prediction for the ratio uk{x)/u-k{x) that agrees with extant data, show- 
ing a strong environment-dependence of the u-quark distribution. The new Drell-Yan 
experiment running at FNAL is capable of validating this comparison, as is the COM- 
PASS II experiment at CERN. Such an experiment should be done so that complete 
understanding of QCD's Goldstone modes can be claimed. 



8 Charting the Interaction between Light-Quarks 

Let us now return to some remarks from Sec. 2.4. Namely, confinement can be related to 
the analytic properties of QCD's Schwinger functions so that the question of light-quark 
confinement may be translated into the challenge of charting the infrared behavior of 
QCD's universal /^-function. This is a well-posed problem whose solution is an elemen- 
tal goal of modern hadron physics [94,242,243]. The answer provides QCD's running 
coupling. Naturally, the behaviour of the /3-function on the perturbative domain is well 
known; indeed, to fourth-loop order [55]. 

The DSEs provide a mapping between experimental observables and the pointwise 
behaviour of QCD's /3-function, including the infrared domain. Hence, comparison be- 
tween computations and observations - of, e.g., the hadron mass spectrum, and elastic 
and transition form factors - can be used to chart the /^-function's long-range behaviour, 
wherefrom arises the pattern of chiral symmetry breaking. Whilst extant studies show 
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Table 8.1. Row 1 - Contact-interaction meson-related results obtained with [117] 
oir/tt = 0.93 in Eq. (2.10), and (in GeV): m = 0.007 in Eq. (2.11), and Ai^ = 0.24, 
Auv = 0.905 in Eq. (2.14). k-„ is the in-pion condensate [129,130]; and f-„^p are the 
mesons' leptonic decay constants. Row 2 ~ Same quantities computed using an inter- 
action based on renormalisation-group-improved one-gluon exchange [118]. Empirical 
values are k^ « (0.22 GeV)^ and [55] f-n = 0.092 GeV, fp = 0.153 GeV. 



interaction 


M 


^TT 


m^ 


■m.p 
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/. 


contact 
QCD-like 


0.37 
0.34 


0.24 
0.24 


0.14 
0.14 


0.93 
0.74 


0.101 
0.093 


0.13 
0.15 



the static properties of hadron excited states to be more sensitive to the long-range be- 
haviour of the /3-function than those of ground states [94,95], dynamical properties of 
ground states possess quite some discriminating power. 

I will now illustrate these remarks by explicating the impact of differing assumptions 
about the behaviour of the Bethe-Salpeter kernel on a few selected hadron properties. In 
this way we will see that it is possible to build a stock of material which can be used to 
identify unambiguous signals in experiment for the pointwise behaviour of: the interac- 
tion between light-quarks; the light-quarks' mass-function; and other similar quantities. 
Whilst these are particular qualities, taken together they will enable a characterisation of 
the nonperturbative behaviour of the theory underlying strong interaction phenomena. 

Consider therefore the interaction kernel specified by Eq. (2.10), and the collection 
of DSEs derived therefrom via the confining, symmetry-preserving regularisation scheme 
indicated in Sec. 2.2. A raft of studies - static properties of vr- and p-meson [77,117,205], 
and those of the neutron and proton [185,205] - have shown that contact-interaction 
results are not realistically distinguishable from those obtained with renormalisation- 
group-improved one-gluon exchange for processes involving momentum transfers Q^ < 
A/^, where M is the dressed-quark mass in Eq. (2.12). This is exemplified in Table 8.1 
and Fig. 8.1 

The picture changes dramatically, however, for processes involving larger momentum 
transfers. This may again be illustrated via the charged pion electromagnetic form factor. 
With a symmetry preserving regularisation of the interaction in Eq. (2.10), no meson's 
Bethe-Salpeter amplitude can depend on relative momentum. Hence Eq. (4.4) becomes 



r,(P) = 75 [iE.[P) + ^7 • PMP)\ 



.1) 



the dressed-quark propagator is given in Eq. (2.12), and the explicit form of the model's 
ladder BSE is 



E^{P) 
F.{P) 



4qir 



ICee ICef 
K,FE ICef 



E„{P) 
F.iP) 



(8.2) 



where, with m = = P^, anticipating the Goldstone character of the pion, anticipating 
the Goldstone character of the pion. 



K.EE = C{M \Tiy:,T^^) , ICeE=0, 2K.EE = Ci{M ;''"ir,Tm,), ICfE 



-2K.I. 



(8.3) 



and Ci{z) — —zC'{z), where I have suppressed the dependence on n^^uv The solution of 
Eq. (8.2) gives the pion's chiral-limit Bethe-Salpeter amplitude, which, for the computa- 
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Figure 8.1. Solid curve - Contact interaction result for the charged pion's electromagnetic 
form factor [117]; and dotted curve, result obtained with renormalisation-group-improved 
one-gluon exchange [120]. There is no meaningful difference for Q^ < K'P . 



tion of observables, should be normalised canonically; viz., 

^M = ^ctr|^r.(-P)^S(g + P)r.(P)S(g). 



Hence, in the chiral limit, 



1 = i^ A^ Ci(^'^^ ^i?' ■^l^)E.[E. - 2F^], 



and the pion's leptonic decay constant is 

iVc 1 



/° = ^^jCi{M'-rt„Tt^)[E. - 2F.\ 



(8.4) 



(8.5) 



(8.6) 



If one has preserved Eq. (3.9), then, for m = in the neighbourhood of P^ = 0, the 
solution of the axial- vector BSE has the form: 



r5M(fc+,fc) = -^2/°r,(P) + 757M^fl(^) 

and the following subset of Eqs. (4.5) - (4.8) will hold: 

/X = M, 2^+Fr = 1. 



(8.7) 



(8.8) 



Hence i<V {P) is necessarily nonzero in a vector exchange theory, irrespective of the point- 
wise behaviour of the interaction. In fact, for m = 0, E-jriP) = 3.568, FTr(P) = 0.459, so 
that F-n (P) has a measurable impact on the value of f-^ . 

More importantly, it has a big impact on the form factor, which may be computed 
using Eq. (7.1). This yields the solid curve in Fig. 7.3, which exposes that a contact 
interaction, Eq. (2.10), generates 



F^'^{Q — >■ <x) = constant. 



(8.9) 
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Equation (8.9) should not come as a surprise: with a symmetry-preserving regularisa- 
tion of the interaction in Eq. (2.10), the pion's Bethe-Salpeter amplitude cannot depend 
on the constituent's relative momentum. This is characteristic of a pointlike particle, 
which must have a hard form factor. Notwithstanding this, for more than twenty-years 
numerous practitioners held to the notion that the pion form factor could not distinguish 
between a contact interaction and QCD. Indeed, this misapprehension still prevails in 
other contexts, as explained in Ref. [44]. 

This illustration emphasises that when a momentum-independent vector-exchange 
interaction is regularised in a symmetry-preserving manner, the results are directly 
comparable with experiment, computations based on well-defined and systematically- 
improvable truncations of QCD's DSEs [120], and perturbative QCD. In this context it 
will be apparent that a contact interaction, whilst capable of describing pion static prop- 
erties well. Table 8.1, generates a momentum-independent dressed-quark mass, which 
entails a form factor whose evolution with Q^ deviates markedly from experiment for 
Q^ > 0.17 GeV^ « M^ and produces asymptotic power-law behaviour, Eq. (8.9), in seri- 
ous conflict with QCD [223-225]. Plainly, experiment is very sensitive to the momentum- 
dependence of the interaction between light-quarks, and the raft of computations reported 
in Refs. [44,77,117,185] therefore provide important benchmarks for use in experimen- 
tally determining QCD's dressed-quark mass function and therefrom constraining the 
/3-function's infrared behaviour. 

It is important now to strive toward the reliable computation of properties of meson 
excited states, and the spectrum of hybrid and exotic mesons, states which are hypo- 
thetically possible in QCD but have not yet been seen experimentally. As mentioned 
above, hybrids are defined to be states that possess normal constituent-quark model me- 
son quantum numbers but non-quark-model decay patterns whilst exotics are described 
by quantum numbers not possible for quantum mechanical quark-antiquark systems. If 
such states exist, then it is suspected that they will possess "constituent gluon" content; 
i.e., they are expected to have a large overlap with interpolating fields that explicitly 
contain gluon fields. This being the case, their existence would completely eliminate the 
distinction between matter- and force-fields that has held since Maxwell's time. Given 
such a composition, however, it is possible that a Bethe-Salpeter equation treatment is 
inadequate to understand their properties. One might instead have to work with a three- 
body Faddeev equation that describes the interactions between two dressed- quarks and 
a dressed-gluon.^ This is a daunting prospect. If one is going to work with the Faddeev 
equation, then it is certainly best to first study systems about which a good deal is known 
experimentally; viz, baryons in QCD. 



9 Describing Baryons and Mesons Simultaneously 

While a symmetry-preserving description of mesons is essential, it is only part of the 
physics that nonperturbative QCD must describe because Nature also presents us with 
baryons: light-quarks in three-particle composites. An explanation of the spectrum of 
baryons and the nature of interactions between them is basic to understanding the Stan- 



^Glueballs, states with no valence quark content at all, are found in numerical simulations of 
lattice-regularised QCD [244]. However, the light-quark current-mass is four-times larger than 
the physical value. Such predictions will not be firm until they are produced with more realistic 
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dard Model. The present and planned experimental programmes at JLab, and other 
facilities worldwide, are critical elements in this effort. 

No approach to QCD is comprehensive if it cannot provide a unified explanation of 
both mesons and baryons. We have explained that DCSB is a keystone of the Standard 
Model, which is evident in the momentum-dependence of the dressed-quark mass function 
- Fig. 2.9: it is just as important to baryons as it is to mesons. Since constituent-quark- 
like models cannot incorporate the momentum-dependent dressed-quark mass-function, 
they are not a viable tool for use in this programme. The DSEs furnish the only extant 
continuum framework that can simultaneously connect both meson and baryon observ- 
ables with this basic feature of QCD, having provided, e.g., a direct correlation of meson 
and baryon properties via a single interaction kernel, which preserves QCD's one- loop 
renormalisation group behaviour and can systematically be improved. This is evident in 
the preceding sections and their combination with Refs. [180,183,211,245]. 

In order to illustrate this programme, I will focus on nucleon electromagnetic form 
factors, which were introduced and described in connection with Eqs. (2.3) - (2.4). As 
noted thereabouts, for a structureless or simply-structured fermion Fi(Q^) = 1 and 
F2{Q^) = 0, so that Ge{Q^) = Gm(Q^) and the distribution of charge and magnetisation 
is identical. This was believed to be the case for the proton until 1999. In that year, 
enabled by the high luminosity of the accelerator at JLab, a new method was employed 
to measure the ratio G^(Q'^)/G'^j^((5^) [246]. The result astonished the community: 
whilst G%{Q'^)/GIj{Q^) ^ 1 for Q2 < lGeV^ G%{Q'^)/Glj{Q^) is a rapidly decreasing 
function for Q^ > I GeV^ (see Fig. 9.6). How is this to be understood? 



9.1 Faddeev equation 

In quantum field theory a baryon appears as a pole in a six-point quark Green func- 
tion. The residue is proportional to the baryon's Faddeev amplitude, which is obtained 
from a Poincare covariant Faddeev equation that sums all possible exchanges and inter- 
actions that can take place between three dressed-quarks. A tractable Faddeev equa- 
tion for baryons [247] is founded on the observation that an interaction which describes 
colour-singlet mesons also generates nonpointlike quark-quark (diquark) correlations in 
the colour-3 (antitriplet) channel [248]. The dominant correlations for ground state octet 
and decuplet baryons are scalar (0^) and axial- vector (1^) diquarks because, for ex- 
ample, the associated mass-scales are smaller than the baryons' masses [169,249] and 
their parity matches that of these baryons. It follows that only they need be retained in 
approximating the quark-quark scattering matrix which appears as part of the Faddeev 
equation [180,181,205]. On the other hand, pseudoscalar (0^) and vector (1~) diquarks 
dominate in the parity-partners of ground state octet and decuplet baryons [205]. The 
DSE approach treats mesons and baryons on the same footing and, in particular, enables 
the impact of DCSB to be expressed in the prediction of baryon properties. 

It is important to appreciate that diquarks do not appear in the strong interaction 
spectrum [97,167,168]. However, the attraction between quarks in this channel justifies 
a picture of baryons in which two quarks within a baryon are always correlated as a 
colour-3 diquark pseudoparticle, and binding is effected by the iterated exchange of roles 
between the bystander and diquark-participant quarks. Here it is important to emphasise 
strongly that QCD supports nonpointlike diquark correlations [117,250]. Hence models 
that employ pointlike diquark degrees of freedom have little connection with QCD. 

The Faddeev equation, illustrated in Fig. 9.1, is a linear homogeneous matrix equa- 
tion, in many respects similar to a Bethe-Salpeter equation. Its solution is the nucleon's 
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Figure 9.1. Diagrammatic representation of a Poincare covariant Faddeev equation for a 
baryon. ^ is the Faddeev amplitude for a baryon of total momentum P = Pq + Pd- It 
expresses the relative momentum correlation between the dressed-quark and -diquarks 
within the baryon. The shaded region demarcates the kernel of the Faddeev equation, 
in which: the single line denotes the dressed-quark propagator, F is the diquark Bethe- 
Salpeter amplitude; and the double line is the diquark propagator. 



Poincare-covariant Faddeev amplitude, which describes quark-diquark relative motion 
within the nucleon. The composite nature of the diquark correlations, and their dy- 
namical breakup and reformation through dressed-quark exchange, pictured in Fig. 9.1, 
is crucial to maintaining fermion statistics for the nucleon bound-state. Furthermore, 
owing to the critical importance of both scalar- and axial-vector-diquark correlations, 
the nucleon's rest-frame wave- function possesses S-, P- and D-watve correlations; i.e., a 
nucleon should a prion be expected to contain significant dressed-quark orbital angular 
momentum. This is verified in Ref. [172], which shows that in the nucleon's rest frame 
just 37% of the total spin of the nucleon is contained within components of the Faddeev 
amplitude which possess zero quark orbital angular momentum. 

It is worthwhile expressing one Faddeev equation concretely and I choose that for the 
A- resonance, a J = 3/2 bound-state of three valence light-quarks. It is possible to obtain 
a simple yet extremely informative equation if one employs the interaction presented 
in Eq. (2.10) and an additional simplification; i.e., representing the quark exchanged 
between the diquarks as 



S^ik)^§: 



(9.1) 



where qa = 1.56 [205]. This is a variant of the so-called "static approximation," which it- 
self was introduced in Ref. [251] and has subsequently been used in studying a range of nu- 
cleon properties [252]. In combination with diquark correlations generated by Eq. (2.10), 
whose Bethe-Salpeter amplitudes are momentum-independent, Eq. (9.1) generates Fad- 
deev equation kernels which themselves are independent of the external relative momen- 
tum variable. The dramatic simplifications which this produces are the merit of Eq. (9.1). 
Following this through, one derives [205] 



1 = 



gl Eqg^^ 
M 



'191 + 



Pa ^iii+ 1 



M 



1? 



X I da {in 
Jo 



WJo 



_dH' 

(2 



da 



("T-991+ + (1 - a)^mi)(amA + M) 



(9.2) 



'99,+ 



+ (1 



mA){amA + M)C'i{aA{a, M, niqq^_^_ , mA)),(9.3) 



Strong QCD and Dyson-Schwinger Equations 



69 



1 1 

• 
■ 


1 1 1 1 1 

■ 


- 


- 


• 

1 1 1 1 1 1 1 1 1 



1 P U, <P a, I* a,* p* f . 




N N'N-N 



A A* A- A*- 



fi fi* fl- n*- 



Figure 9.2. Comparison between DSE-computed hadron masses {filled circles) and: bare 
baryon masses from Ref. [253] {filled diamonds) and Ref. [254] {filled triangles); and ex- 
periment [55], filled-squares. For the coupled- channels models a symbol at the lower 
extremity indicates that no associated state is found in the analysis, whilst a symbol 
at the upper extremity indicates that the analysis reports a dynamically-generated res- 
onance with no corresponding bare-baryon state. In connection with Sl-baryons the 
open-circles represent a shift downward in the computed results by 100 MeV. This is an 
estimate of the effect produced by pseudoscalar-meson loop corrections in A-like systems 
at a s-quark current-mass. 



where Eqq ^ is the canonically-normalised axial-vector diquark Bethe-Salpeter ampli- 
tude, mqq , is the computed mass for this correlation. 
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(9.4) 



and Ci (z) — Ci^{z)/z [see text around Eqs. (2.15), (8.3)]. Equation (9.3) is an eigenvalue 
problem whose solution yields the mass for the dressed-quark-core of the A-resonance, 
rriA. It is one dimensional because only the axial- vector diquark correlation contributes 
to the structure of the A. (The nucleon, on the other hand, is constituted from scalar- 
and axial- vector-diquarks and presents a five-dimensional eigenvalue problem.) 

With experience, one can look at Eq. (9.3) and see that increasing the current-quark 
mass will boost the mass of the bound-state. This is just one of the Faddeev equations 
in Ref. [205]. In fact, building on lessons from meson studies [14], a unified spectrum 
of u, d-quark hadrons was obtained therein using the symmetry-preserving regulariza- 
tion of a vector x vector contact interaction that I have briefiy described herein. Refer- 
ence [205] reports a study that simultaneously correlates the masses of meson and baryon 
ground- and excited-states within a single framework. In comparison with relevant quan- 
tities, the computation produces fnTs=13%, where mis is the root-mean-square-relative- 
error/degree-of freedom. As evident in Fig. 9.2, the prediction uniformly overestimates 
the PDG values of meson and baryon masses [55]. Given that the employed truncation 
deliberately omitted meson-cloud effects in the Faddeev kernel - analogues of those in 
the right image of Fig. 7.1 - this is a good outcome, since inclusion of such contributions 
acts to reduce the computed masses. 

Following this line of reasoning, a striking result is agreement between the DSE- 
computed baryon masses [205] and the bare masses employed in modern coupled-channels 
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Figure 9.3. The Excited Baryon Analysis Center (EBAC) examined the Pn-channel 
and found that the two poles associated with the Roper resonance and the next higher 
resonance were all associated with the same seed dressed-quark state. Coupling to the 
continuum of meson-baryon final states induces multiple observed resonances from the 
same bare state. In EBAC's analysis, all PDG-identified resonances were found to consist 
of a core state plus meson-baryon components. (Adapted from Ref. [253].) 



models of pion-nucleon reactions [253,254], see Fig. 9.2 and also Ref. [185]. The Roper 
resonance is very interesting. The DSE study [205] produces an excitation of the nucleon 
at 1.82 ± 0.07 GeV. This state is predominantly a radial excitation of the quark-diquark 
system, with both the scalar- and axial- vector diquark correlations in their ground state. 
Its predicted mass lies precisely at the value determined in the analysis of Ref. [253] . This 
is significant because for almost 50 years the "Roper resonance" has defied understanding. 
Discovered in 1963, it appears to be an exact copy of the proton except that its mass 
is 50% greater. The mass was the problem: hitherto it could not be explained by any 
symmetry-preserving QCD-based tool. That has now changed. Combined, see Fig. 9.3, 
Refs. [205,253] demonstrate that the Roper resonance is indeed the proton's first radial 
excitation, and that its mass is far lighter than normal for such an excitation because 
the Roper obscures its dressed-quark-core with a dense cloud of pions and other mesons. 
Such feedback between QCD-based theory and reaction models is critical now and for the 
foreseeable future, especially since analyses of experimental data on nucleon-resonance 
electrocouplings suggest strongly that this structure is typical; i.e., most low-lying A^*- 
states can best be understood as an internal quark-core dressed additionally by a meson 
cloud [255,256]. 

Additional analysis [185] suggests a fascinating new feature of the Roper. To elu- 
cidate, I focus first on the nucleon, whose Faddeev amplitude describes a ground-state 
that is dominated by its scalar diquark component (78%). The axial- vector component 
is significantly smaller but nevertheless important. This heavy weighting of the scalar 
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diquark component persists in solutions obtained with more sophisticated Faddeev equa- 
tion kernels (see, e.g., Table 2 in Ref. [181]). From a perspective provided by the nucleon's 
parity partner and the radial excitation of that state, in which the scalar and axial- vector 
diquark probabilities are [257] 51%-49% and 43%-57%, respectively, the scalar diquark 
component of the ground-state nucleon actually appears to be unnaturally large. 

One can nevertheless understand the structure of the nucleon. As with so much else, 
the composition of the nucleon is intimately connected with dynamical chiral symmetry 
breaking. In a two-color version of QCD, the scalar diquark is a Goldstone mode, just 
like the pion [165]. (This is a long-known result of Pauli-Giirsey symmetry [258,259].) A 
memory of this persists in the three-color theory and is evident in many ways. Amongst 
them, through a large value of the canonically normalized Bethe-Salpeter amplitude 
and hence a strong quark-f-quark— diquark coupling within the nucleon. (A qualitatively 
identical effect explains the large value of the nN coupling constant.) There is no such 
enhancement mechanism associated with the axial-vector diquark. Therefore the scalar 
diquark dominates the nucleon. 

With the Faddeev equation treatment described herein, the effect on the Roper is 
dramatic: orthogonality of the ground- and excited-states forces the Roper to be consti- 
tuted almost entirely (81%) from the axial-vector diquark correlation. It is important 
to check whether this outcome survives with a Faddeev equation kernel built from a 
momentum-dependent interaction. 



9.2 Nucleon form factors 

With masses and Faddeev amplitudes in hand, it is possible to compute baryon electro- 
magnetic form factors. For the nucleon, studies of the Faddeev equation exist that are 
based on the one-loop renormalisation-group-improved interaction that was used effica- 
ciously in the study of mesons [180,181]. These studies retain the scalar and axial- vector 
diquark correlations, for the reasons explained in Sec. 9.1. 

A nonpointlike composite nucleon must interact with the photon via a sophisticated 
current, whose form is constrained by vector Ward-Takahashi identities. This continues a 
thread that pervades these notes. For the bound-state described by the Faddeev equation 
in Fig. 9.1 that current is described in Ref. [260] and depicted in Fig. 9.4: Diagrams 4-6 
represent eight-dimensional integrals, which can be evaluated using Monte-Carlo tech- 
niques. 

The pattern of the computation should now be clear. In principle, one specifies an 
interaction kernel for the gap equation and solves for the dressed-quark propagator. This 
completes the specification of the Bethe-Salpeter kernel, so one may compute the masses 
and amplitudes for the diquark correlations. With all these quantities determined, the 
Faddeev equation is defined and can be solved for a baryon's mass and amplitude. In 
combination with the current. Fig. 9.4, it is a straightforward numerical task to compute 
the form factors. 

Some of the steps indicated here can be simplified by employing algebraic parametri- 
sations of one or more elements (propagators and/or amplitudes). Such parametrisa- 
tions, based upon numerical solutions of the gap- and Bethe-Salpeter equations, have 
long been employed efhcaciously [179,234,262]. Following that path, Ref. [181] produced 
a comprehensive survey of nucleon electromagnetic form factors, the quality of which is 
characterised by the results depicted in Fig. 9.5. With these results, Ref. [181] unified the 
computation of meson and nucleon form factors, and also their valence-quark distribution 
functions (see Sec. 9.3). 




Figure 9.4. Vertex which ensures a conserved current for on-shell nucleons described 
by the Faddeev amplitudes, ^Pij, obtained from the equation depicted in Fig. 9.1. The 
photon probing the nucleon introduces momentum Q and is represented by the wiggly 
line; the single line represents S(p), the dressed-quark propagator, and the double line, 
the diquark propagator; and F is the diquark Bethe-Salpeter amplitude. The vertex in 
Diagram 4 (second row, right column) represents an electromagnetically induced transi- 
tion between the scalar and axial-vector diquarks; and the vertices in row 3 are "seagull 
terms," which appear as partners to Diagram 3 and arise because binding in the nucle- 
ons' Faddeev equations is, in general, effected by the exchange of a momentum-carrying 
dressed-quark between nonpointlike diquark correlations [260]. 



Figure 9.5 also depicts nucleon form factors obtained [185] with a symmetry-preserving 
DSE-treatment of the contact interaction in Eq. (2.10). These form factors charac- 
terise a nucleon that is constructed from diquarks whose Bethe-Salpeter amplitudes 
are momentum-independent and dressed-quarks with a momentum-independent mass- 
function, which inputs to the Faddeev equation yield a bound-state described by a 
momentum-independent Faddeev amplitude. This last is the hallmark of a pointlike 
composite particle and explains the hardness of the computed form factors, which is evi- 
dent in Figs. 9.5. The hardness contrasts starkly with results obtained from a momentum- 
dependent Faddeev amplitude produced by dressed-quark propagators and diquark Bethe- 
Salpeter amplitudes with QCD-like momentum-dependence; and with experiment. As 
I have remarked above, evidence for a connection between the momentum-dependence 
of each of these elements and the behaviour of QCD's /3-function is accumulating; e.g., 
Refs. [44,77,117,120,180,183,189]. The comparisons in Figs. 9.5 add to this evidence, in 
connection here with readily accessible observables, and support the view, presented 
herein, that experiment is a sensitive probe of the running of the /3-function to infrared 
momenta. 
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Figure 9.5. Nucleon Dirac (upper panel) and Pauli (lower panel) form factors, as a func- 



tion of a; = Q /r 



proton on left and neutron on right, dashed curve - result obtained 



in Ref. [181], which employed QCD-like momentum-dependence for the dressed-quark 
propagators and diquark Bethe-Salpeter amplitudes in solving the Faddeev equation; 
Solid curve - result obtained using a symmetry-preserving regularisation of a contact- 
interaction [185] and hence a dressed-quark mass-function and diquark Bethe-Salpeter 
amplitudes that are momentum-independent; and dot- dashed curve - a parametrisation 
of experimental data [261]. 



This returns us to the ratio G^((5^)/G^,j^(Q^) described in the introduction to this 
section and depicted in Fig. 9.6. One observes that, independent of the nature of the 
interaction, computations of this ratio exhibit a zero. Its location, however, is very 
sensitive to the interaction, although not to the electromagnetic size of the diquark 
correlations [181]. (I note that the DSE prediction of the same ratio for the neutron has 
been confirmed by recent experiments [268], as also is the trend of the nucleon's flavour 
separated Dirac and Pauli form factors [185,268,269].) 

In order to assist in explaining the origin and location of a zero in the Sachs form 
factor ratio, in the top panel of Fig. 9.6 I depict the ratio of Pauli and Dirac form factors: 
both the actual contact-interaction result and that obtained when the Pauli form factor 
is artificially "softened;" viz., 



F'zpiQ^ 



F2p(Q^ 



l-f 



r/(4m%) 



(9.5) 



As observed in Ref. [270] , a softening of the proton's Pauli form factor has the effect of 
shifting the zero to larger values of Q^. In fact, if F2p becomes soft quickly enough, then 
the zero disappears completely. 

The Pauli form factor is a gauge of the distribution of magnetisation within the 
proton. Ultimately, this magnetisation is carried by the dressed-quarks and influenced 
by correlations amongst them, which are expressed in the Faddeev wave-function. If 
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Figure 9.6. Upper panel: Normalised ratio of proton Pauli and Dirac form factors. 
Solid curve - contact interaction [185]; long-dashed curve - result from Ref. [199], which 
employed QCD-like momentum-dependence for the dressed-quark propagators and di- 
quark Bethe-Salpeter amplitudes; long-dash-dotted curve - drawn from parametrisation 
of experimental data in Ref. [261]; and dotted curve - softened contact-interaction result, 
described in connection with Eq. (9.5). Lower panel: Normalised ratio of proton Sachs 
electric and magnetic form factors. Solid curve and long-dashed curve, as above; dot- 
dashed curve - linear fit to data in Refs. [246,263-266], constrained to one at Q^ = 0; 
short-dashed curve - [1, 1]-Pade fit to that data; and dotted curve - softened contact- 
interaction result, described in connection with Eq. (9.5). In addition, we have repre- 
sented a selection of data explicitly; filled-squares [263]; circles [265]; up-triangles [266]; 
and open-squares [267]. 



the dressed-quarks are described by a momentum-independent mass-function, then they 
behave as Dirac particles with constant Dirac values for their magnetic moments and 
produce a hard Pauli form factor. Alternatively, suppose that the dressed-quarks possess 
a momentum-dependent mass-function, which is large at infrared momenta but vanishes 
as their momentum increases. At small momenta they will then behave as constituent- 
like particles with a large magnetic moment, but their mass and magnetic moment will 
drop toward zero as the probe momentum grows. (N.B. As described in Sec. 6.2, massless 
fermions do not possess a measurable magnetic moment.) Such dressed-quarks will pro- 
duce a proton Pauli form factor that is large for Q^ ~ but drops rapidly on the domain 
of transition between nonperturbative and perturbative QCD, to give a very small result 
at large-Q^ . The precise form of the Q^-dependence will depend on the evolving nature 
of the angular momentum correlations between the dressed-quarks. From this perspec- 
tive, existence, and location if so, of the zero in fJ-pGEp{Q^)/GMp{Q^) are a fairly direct 



dv{x) 


, where 


dy{x) 
u^{x) 




Uvix) 


. F2(-) 






^ F^M 



Strong QCD and Dyson-Schwinger Equations 75 

measure of the location and width of the transition region between the nonperturbative 
and perturbative domains of QCD as expressed in the momentum-dependence of the 
dressed-quark mass- function. 

One may expect that a mass-function which rapidly becomes partonic - namely, is 
very soft - will not produce a zero; has seen that a constant mass-function produces a zero 
at a small value of Q^, and knows that a mass- function which resembles that obtained in 
the best available DSE studies [94,99] and via lattice-QCD simulations [101], produces a 
zero at a location that is consistent with extant data. There is plainly an opportunity here 
for very constructive feedback between future experiments and theory. It is anticipated 
that experiments at JLab in the 12 GeV era will establish conclusively whether or not 
this ratio possesses a zero. The result will assist greatly in refining understanding of the 
dressed-quark mass function and therefrom QCD's /3-function. 



9.3 Valence-quark distributions at a; = 1 

Before closing this section I would like to exploit a connection between the Q^ = 
values of elastic form factors and the Bjorken-a; = 1 values of the dimensionless structure 
functions of deep inelastic scattering, F2'^{x). First recall that the x = 1 value of a 
structure function is invariant under the evolution equations. Sec. 7.2.1. Hence the value 
of 

. ..^ 4SM - 1 

(9.6) 

is a scale-invariant feature of QCD and a discriminator between models. Next, when 
Bjorken-i is unity, then Q^ + 2P ■ Q = 0; i.e., one is dealing with elastic scattering. 
Therefore, in the neighbourhood of a; = 1 the structure functions are determined by the 
target's elastic form factors. The ratio in Eq. (9.6) expresses the relative probability of 
finding a d-quark carrying all the proton's light-front momentum compared with that of 
a it-quark doing the same or, equally, owing to invariance under evolution, the relative 
probability that a Q^ = probe either scatters from a d-quark or a it-quark; viz., 

^^ =^ (97) 

In constituent-quark models with 5f/(6)-symmetric spin-flavour wave-functions the 
right-hand-side of Eq. (9.7) is 1/2 because there is nothing to distinguish between the 
wave-functions of u- and d-quarks, and the proton is constituted from it-quarks and one 
d-quark. On the other hand, when a Poincare-covariant Faddeev equation is employed 
to describe the nucleon, 

pp,d 2 pp,a I 1 pp,m 

f_l ^ 3-^1 """3-^1 /q c,-. 

pp,u pP,s I 1 r>P,a I 2 r>P,m i ^^-"1 

^1 ^1 + 3-^1 + 3-^1 

where I have used the notation of Ref. [181]. Namely, P['" — F(p{Q^ — 0) is the con- 
tribution to the proton's charge arising from diagrams with a scalar diquark component 
in both the initial and final state: u[ud\ (g) 7 (g) M[ud]. The diquark-photon interaction is 
far softer than the quark-photon interaction and hence this diagram contributes solely 
to Ui, at x = 1. P*"'" — Fip{Q^ = 0), is the kindred axial-vector diquark contribution; 
viz., 2d{uu} (g) 7 (g) d{uu} + u{ud} (gi 7 (g) u{ud}. At a; = 1 this contributes twice as much 
to dv as it does to Uv P['"^ = Fi^{Q^ = 0), is the contribution to the proton's charge 
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arising from diagrams with a different diquark component in the initial and final state. 
The existence of this contribution relies on the exchange of a quark between the diquark 
correlations and hence it contributes twice as much to u„ as it does to d„. If one uses 
the "static approximation" to the nucleon form factor, Eq. (9.1), as with the treatment 
of the contact-interaction in Ref [185], then P^'"^ = 0. It is plain from Eq. (9.8) that 
dv/uv = in the absence of axial- vector diquark correlations; i.e., in scalar-diquark-only 
models of the nucleon, which were once common and, despite their weaknesses, still too 
often employed. 

Using the probabilities presented in Refs. [181,185], respectively, one obtains: 

pP,s pP,a pP.m d^ F^ 

1 1 1 u^ F^ 

0.60 0.25 0.15 0.28 0.49 . (9-9) 



"M(7) 
M=constant 



0.78 0.22 0.18 0.41 



Both rows in Eq. (9.9) are consistent with dv/uv — 0.23 ± 0.09 (90% confidence level, 
F2 I ^2 ~ 0-45 i 0.08) inferred recently via consideration of electron-nucleus scattering 
at a; > 1 [271]. On the other hand, this is also true of the result obtained through a 
naive consideration of the isospin and helicity structure of a proton's light-front quark 
wave function at x ~ 1, which leads one to expect that rf-quarks are five-times less likely 
than u-quarks to possess the same helicity as the proton they comprise; viz., d^/u^ = 0.2 
[228]. Plainly, contemporary experiment-based analyses do not provide a particularly 
discriminating constraint. Future experiments with a tritium target should help [272], 
emphasising again the critical interplay between experiment and theory in elucidating 
the nature of the strong interaction. 



10 Epilogue 

QCD is the most interesting part of the Standard Model and Nature's only example 
of an essentially nonperturbative fundamental theory. Whilst confinement remains a 
puzzle, it is recognised that dynamical chiral symmetry breaking (DCSB) is a fact in 
QCD. It is manifest in dressed-propagators and vertices, and, amongst other things, it is 
responsible for: the transformation of the light current-quarks in QCD's Lagrangian into 
heavy constituent-like quarks, in terms of which order was first brought to the hadron 
spectrum; the unnaturally small values of the masses of light-quark pseudoscalar mesons 
and the 77-7;' splitting; the unnaturally strong coupling of pseudoscalar mesons to light- 
quarks - Qjtqq ~ 4.3; and the unnaturally strong coupling of pseudoscalar mesons to the 
lightest baryons - fl.^jviv ~ 12.8 ^ Sgnqq- 

Herein I have illustrated the dramatic impact that DCSB has upon observables in 
hadron physics. A "smoking gun" for DCSB is the behaviour of the dressed-quark mass 
function. The momentum dependence manifest in Fig. 2.9 is an essentially quantum 
field theoretical effect. Exposing and elucidating its consequences therefore requires a 
nonperturbative and symmetry-preserving approach, where the latter means preserving 
Poincare covariance, chiral and electromagnetic current-conservation, etc. The Dyson- 
Schwinger equations (DSEs) provide such a framework. I have explained the nature of 
some of the experimental and theoretical studies that are underway which can potentially 
identify observable signals of M{p'^) and thereby confirm and explain the mechanism 
responsible for the vast bulk of visible mass in the Universe. 
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Along the way I have described a number of exact results proved in QCD using the 
DSEs, amongst them: 

• Light-quark confinement is a dynamical phenomenon, which cannot in principle be 
expressed via a potential; 

• Goldstone's theorem is fundamentally an expression of equivalence between the 
one-body problem and the two-body problem in the pseudoscalar channel; 

• quarks are not Dirac particles - they possess anomalous chromo- and electro- 
magnetic moments which are large at infrared momenta; 

• and gluons are nonperturbatively massive, being described by a mass-function 
which is large in the infrared but diminishes with power-law behaviour in the 
ultraviolet. 

Numerous items could be added to this list, some of which are described above. 

There are many reasons why this is an exciting time. As per my brief, I have fo- 
cused on one. Namely, through the DSEs, one is unifying phenomena as apparently 
diverse as: the hadron spectrum; hadron elastic and transition form factors, from small- 
to large-Q^; and parton distribution functions. The key is an understanding of both 
the fundamental origin of visible mass and the far-reaching consequences of the mecha- 
nism responsible; namely, DCSB. Through continuing feedback between experiment and 
theory, these studies should lead us to an explanation of confinement, the phenomenon 
that makes nonperturbative QCD the most interesting piece of the Standard Model. 
They might also provide an understanding of nonperturbative physics that enables the 
formulation of a realistic extension of that model. 
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